AVISION 


. 
7 


 MOLUME 87 NO. 
— 
iil & AMERICAN SOCIETY 


BASIC REQUIREMENTS FOR "MANUSCRIPTS 
riginal papers discussions of current papers should ‘ie submitted to 
_ Manager of Technical Publications, ASCE. . Authors should indicate the technical _ 
_ division to which the paper is referred. The final date on which a discussi on should 


teach the Society is given as a footnote with each paper. Those who are ae e 
i yl ‘to submit material will expedite the review and publication 


1. Titles must have a length not esol 50 characters and spaces. 


3 it. 
on one side of 8%4-inch by 11-inch ‘paper. 


all illustrations, tables, etc., must be included. 
The author's full name, Society membership grade, 


ing present em loyment must appear on the first age of the | aper. 


Mathematics are recompose ed from m the Be 
of this, it is necessary that letters be drawn carefully, and that special symbols be — 
properly identified. The letter used. should be defined where 


in the illustrations or 


‘copita), ca one side of ‘inch by 11-inch paper. Speci ill 
must be aide ‘in the text each tabl 


must wi ha in black ink on one side of 8 54-inch by 1 1l- inch 


paper Because illustrations will be reproduced with a width of between 3- width, 


ext for each 


The desi is 18,000 As an each full page of 


text, or illustration is the equivalent of 300 words. 
Re; Technical papers ded ablicati on must be written in the third 


d a* 
this may be made on condition that the full title, name 


e, and da e of publication by the ‘ 


‘published by the America Society of Civi 
r Engineers. Publication office is at 2500 South State Street, Ann Arbor, Michigan. 
Editorial and General | Offices are at 33 West 39 Street, New York 18, New York. 


‘The index for 1959 was published as ASCE Publication 1960- 10 (list pri 
$2. 00); ; indexes for previous years are also available. ; : 


— 
__2. Asummary of approximately 50 words must accompany the paper, a 300- be 
— 
— 
— 
— 
— 
— 
— 
— 
otograpns snou. pe su Xpianations ana Gescr 1p 
— 
— 
— 
— 
— «26 
— 
— 
— 
ae Thee Society 1s responsible TOF aly stavement maue 
— 
— 
— 
— 
— 


Further Studies of the Strength of Beam- -Columns 
re Robert L. Ketter .... 135 


87 No. ST 


Journal o of the 


STRUCTURAL. DIVISION | 


eedings of the American Society of Civil Engineers — 


EXECUTIVE COMMITTEE 


a _ Emerson J. Ruble, Chairman; Nathan D. Whitman, Jr., Vice Chairman; 


- Robert D. Dewell; Theodore R. ‘Higgins; John D. Haltiwanger, Secretary 
aan _ Elmer K. Timby, Board Contact Member ; 


ON N PUBLICATIONS 


CONTENTS 


96 
August, , 1961 


Torsional Behavior of of Suspension Bridge Towers aan 
= Frank Baron and Anthony G. Arioto 


Dead- Load Stress in Model Dams by Me Method of of Integration — at on 


‘General Solution of Space Frameworks 
Ignacio Martin and Jose E. 


Mournir Badir 


Copyright 1961 by the American. Society of Civil Engineers. 

_ __ The three preceding issues of this Journal are dated March 1961, pril 


a 

| 
| 


if 
| 
Gerald Abbott Frank; 4 
Wil 
4 
| 
1 
| 


Modes of the Fixed Base Column, by Donald 
‘Sawyer. (July, 1960. Prior discussion: October, 1960, March, 1961 


a Donald A. A. Sawyer | (closure) . 


Response of Elasto- Plastic Frames, by Joseph Penzien. 
Guly, 1960. Prior discussion: June, 1960, March, 


_ by Joseph Penzien (closure) 


; Freezing and Thawing Effects on Prestressed ‘Concrete, by M. J. 


Gutzwiller and F. E. Musleh. (October, 1960. Prior discussion: 
March, 1961. Discussionclosed.) 
by M. J. Gutzwiller and F. E. Musleh (closure) . wid? 


in Structural , by Kurt H. Gerstle. 


(October, 1960. Prior discussion: ‘March, April, 1961. 


Concepts of Structural Safety, by B. Brown. (December, 
1960. Prior discussion: April, (1961. Discussion closed. 
by George B. Begs, Jr. 


"Lateral Inelastic ‘Buckling of Tied Arches 5, by Chil Chin Fung Kee. 
 Canuary, 1961. Prior discussion: None. Discussion | closed. 


"Safety, Reliability and ‘Structural De Design, by A. M. Freudenthal. 
(March, 1961. Prior discussion: None. Discussion closed.) an 


by Jack R. and Theodore C. ‘Zsutty. 


Strength of PlateGirders in Bending 
by Konrad Basler and Bruno Thirlimann ........ 
— 
yf 
: 
(185 — 
= 
-~_ 
— 
d 
= 
= 
| 


STRUCTURAL 


By F Frank Baron, 1 Fr. ASCE and Anthony G. Arioto,2 A 


= Pes 
re 


oo 
SYNOPSIS 
- Sa An interpretive study is presented of the influences of different vari bles on a 
the torsional behaviors of suspension . bridge towers. General and specific 
classes of towers are considered. Among the variables are the numbers of _ 
struts, proportions of panels, and relative sizes and struts. -Dimen- 
sionless ratios are used to relate the behaviors. : ; 


INTRODUCTION 
The torsional characteristics of suspension bridge towers are of a 


in studies dealing with the effects of (a) traffic loads on one side of the road- 
way, and (b) aerodynamic forces which tend to twist the roadway of a bridge. 
_ Each of these loading conditions produces a torque at the top of a tower about 
_ the longitudinal axis of the tower. The torque produces an angle of twist =~ 
= it is the result of the unbalanced horizontal components of cable tension — 
at the top of the tower. Its magnitude is dependent on the characteristics of . 
the loads and of the bridge. With reference to those of the bridge, , the ~ll 
4 on a tower is dependent on the torsional rigidity of the tower and on the re- | 
a lative resistances 0 of the side and center span cables to to Caner & at - 
~ é _ The torsional rigidities of the towers can influence the deformations of a 
the towers, the roadway, the stiffening trusses, and the bracing elements = 
_ Note.—Discussion open until January 1, 1962. To extend the closing date one month, i 
a written request must be filed with the Executive Secretary, ASCE. This paper is part | 
of the copyrighted Journal of the Structural Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No. ST 6, August, 1961. 
a Prof. of Civ. Engrg., Univ. of California, Berkeley, Calif. 
2 Assoc. Civil Engr., County of Santa Clara, Calif. 
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TORSIONAL BEHAVIOR OF SUSPENSION BRIDGE TOWERS | 
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the roadway system of a The influence of these factors pre- 
a sented in| a) a report of the Chief Engineer of the Golden Gate Bridge, and 
(ob) in the —— of the Boards of Engineers4 appointed to study the causes 
was yiven ina previous paper by the writers. 
The torsional rigidities towers and with the distributions o of bending 
moments, twisting moments, and shears in the legs and struts of the towers — 
are analyzed. ‘Two groups of towers are considered, One group c consists of i 
actual towers, specifically those of the Golden Gate Bridge, the first Tacoma 
Narrows Bridge, and the Walt Whitman Bridge. The other group consists of 
: hypothetical towers. Among the variables are the height to width ratio of a 
tower, the number of struts inthe tower, and the flexural and torsional rigidi- 7 
ties of the legs and struts. - 
Results obtained from an analytical procedure, presented elsewhere, 5 for 7m 
determining the torsional characteristics of towers are given, The ccc 


> consists of beginning at the base of a tower and obtaining successive statical 
and geometrical relationships at each section of the tower until the top of the 
a : tower is reached, For the actual towers, the procedure was adjusted to con- _ 8 


sider the effects of the vertical saddle reactions at the tower tops, variations - 
a in I and J of the legs, and variations in the canes distances between the * 


_ The notations are the same as those listed a pan for a few 
modifications where necessary. The notations for the actual towers and the 
‘hypothetical towers are illustrated in ‘Figs. 3. 1(a) and 1(b), respectively. The 


horizontal distance between center r lines of ath or oa 
ates of struts in a tower; 7 


Young’s modulus = 30 x 103 k in.-2; 


‘modulus, which is Poisson’s ratio (u= 0.3); 
= moment of inertia (Subscripts designate axis, location 


section, or kind of member); 


= torsional rigidity y of a member, in which J is taken shout 


the appropriate axis of the member; 


height | of a panel at the base of a tower; ; 


_ 3 The Golden Gate Bridge,” by Joseph B. Strauss, Bd. of Directors, Golden | Gate 
and Highway Dist., California, September, 1937, pp. 90-107. 
4 *The Failure of the Tacoma Narrows Bridge,” Texas Engrg. Experiment Sta., a 
 letin No. 78, 1944. Part (a): Report to John M. Carmody by Bd. of Engrs., Othmar H. _ 
_ Amman, Theodore von Karman, and Glenn B. Woodruff, March 28, 1941, pp. 61, 62, 125, 
and 133. Part (b): Reports to Paul Carew by Bd. of Narrows Bridge Loss Committee, 
_ Clifford E. Paine, Chmn., Hardy Cross, ‘Shortridge Hardesty, Holton D. Robinson, and — 
Wilbur M. Wilson, June 2, 1941, pp. 42-48. Part (c): _ Report to James A. Davis by —_ 
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“Torsional Analysis of Suspension Bridge Towers,” by Frank Baron and AnthonyG. 
Bi a Arioto, Proceedings, ASCE, Vol. 86, No. ST 1, January, 1960, pp. 143-169. ae &g 
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FIG. 1.—DIMENSIONS | 
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= flexural eaninns of a member, in whichI is taken about the 7 
appropriate axis of the mer member; 


* torsional rigidity of of the entire re tower; 


=transverse shear along the z- axis (Additional eal 
designate the section | or member); 


gnate the section 
= bending moments or about x and y axes, 
_ respectively, depending on whether the member isa leg or -. 

strut ‘subscripts designate the section or ‘mem-— 
ber); 


= 
= torque about the longitudinal axis of the tower; = 
= vertical reaction at each saddle of f the tower; a 


ee ‘reaction ; at each saddle, c caused by the dead load 
of the ‘suspended system; 


= an Euler load for each leg 
z of a tower (The buckling load is for a leg fixed at the base 


and free at the is about the — 


direction; 


= angle of twist of a at the 
relative to the base and about the y-axis. | 
In the section of the paper which deals with the hypothetical towers, addi- 
tional notations are used, giving dimensionless ratios as follows: = | 


| 


the s subscripts and L refer | to the legs of the towers, 
* respectively. Note that h is the same for each panel of an hypothetical tower, 


he is the height- to-width ratio of the panel, an an h/s = rh for an hypothetical 


wl sign convention is illustrated in Fig. 2. The shears are represented - 


| 
= 
— 


and y-axes, positive angle of twist dy ‘is associated with a positive torque 


TORSIONAL CHARACTERISTICS OF ACTUAL TOWERS 
The towers are those of the Golden Gate ate Bridge, the first Tacoma Narrows 
‘Bridge, and the Walt Whitman Bridge. These towers were selected for 7 
because of the wide range which they represent in tower parameters. The only 7 
comparisons that are made are those the torsional characteristics 


(0) GOLDEN GATE 

~MACKINAC “STRAITS 

F FRANCISCO - OAKLAND 

) FIRST TACOMA NARROWS” 


WALT 
= DEALT WITH IN TEXT 
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Of the three towers. Obviously, the design of of ‘each tower is is based on on many 
considerations besides those of torsion. _ 
Description of the Towers.—The outlines and over-all dimensions of the | 
three towers are shown in Fig. 3. If additional details are desired, the reader — 


is referred to the plates and figures which | are given elsewhere. 6,7,8 — 


_ 6 The Golden Gate Bridge,” by Joseph B. Strauss, Bd. of Directors, Golden Gate 
“Bridge and Highway Dist., California, September, 1937, pp. 90-107, Plates V to IX. 

‘s . “The Failure of the Tacoma Narrows Bridge,” Texas Engrg. Experiment Sta., Bul- 
letin No. 78, 1944, Part (c), Figs. 2, 5, 5(a); Part (b), Fig. 2. 
_ 8 “Anchorages and Superstructure of Walt Whitman Bridge, ge,” by Homer R. Seely, 
Civil Engineering, February,1956, 52-57 
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tower is composed “legs which 2 are and two or more 

- struts. The legs of each tower are cellular in form and vary in cross-section | a 
from the base to the top of the tower. The cross-sections at the bases of the 
towers are illustrated in Fig. 4. The variation in cross- -section along tt the length 


ia, plates and in the sizes of the scien second, (a) by a change in the number of 
cells for the Golden Gate towers, (b) by a change in the transverse | dimensions © 
a of the legs and a change in the number of cells for the first Tacoma Narrows 


FIRST TACOMA. 
NARROWS 


a GOLDEN GATE "WHITMAN 


FIG. 5. _cRoss- SECTIONS OF STRUTS 


a, The numbers and cross- -sections of the struts are different for the three 
_e towers and are illustratedin Figs. 5to 8. The struts of the Golden Gate towers 

a are composed of trusses, whereas those of other towers are composed of 

= thin-walled cells. The trusses of the three top struts of the Golden Gate Bridge _ 
n towers are laced together in pairs. It is reported3 that the lacing of the web_ 
_ members was increased as a result of an mveatigation of the effects of a 
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“BRIDGE TOWERS 
wl Dimensions and Properties of Towers.—A summary summary is given in Table 1 of 
the panel _ proportions and of the horizontal and vertical dimensions of the — 


Other factors which influence the torsional characteristics of the towers are 
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FIG. 7.—TACOMA NARROWS TOWERS 
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TOWERS 


and for a strut ne 06% 25 


and 


In Eqs. the to the axes that are shown in Fig. 
Values of < * and J were eee for the top and bottom sections ct 
every tier in the elegs (Figs. 6 to 8), and for the ends and center cross- sections 


1.. —SUMMARY OF TOWER DIMENSIONS 


Tecoma 
‘Narrows 


TABLE 2, VALUES OF I AND J 


wa 


Sec, 01 ‘Bec. rq Strut 
__ Values of | Values of Valuesof |  Valuesof 


0 in’ in4 (108in 4) (108in 


First Tacoma 


7 - of the struts. The | distributions: of the values of I and J along the legs of the - 
towers are indicated in Figs. 6 to 8; and representative values of I and J are 

7 _ ated in Table 2. The representative values are for the midsections of the top 

4 ~ and bottom struts, and for the top andbottom sections of the legs. Figs. 6 to 8 
& _ show that the values of I;, and Jy, are greatest for the legs of the Golden Gate 
“oa towers , and ‘smallest for a. of the first Tacoma Narrows towers. Table 
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are greatest for the Walt Whitman towers, and smallest for the Golden Gate 
<a The computations of I are based on the gross areas of the cross- -sections. \- 
The computations of J are based on the gross areas of the sections minus one- 
‘half of the areas of the stiffener angles. In Table 2, equivalent values of I and 


| are given for the trussed struts of the Golden Gate towers. _ ma 
er The computations for J are basedon the assumptions of the ordinary theory 


a of torsion, and on the application of Bredt’s formula9,10 to each section of the 


, towers. For example, consider the section that is shown in Fig. 9. The section 
is subjected to a total torque T, and is composed of walls of equal thickness © 


and of five equal cells. The walls are so thin in proportion to the transverse _ : 
dimensions of the cells that the stresses t/t canbe assumed 


| 


t is smal 
to 
and d ‘a 


4 


“4 OF CROSS- SECTION IS shown 


“FIG. $ ear 10.- SINGLE CELL 


distributed across the of the of the walls is ont 
_ included in this study. A distribution of forces f per unit of lengthas shown in > 
iH the figure satisfies the requirements of static equilibrium because the quantity - 


of force f flowing into a junction equals the quantity of force f tee 7 amy 


“Theory of by S. . Timoshenko, Book Co., New York, 1934, 

_ 10 “Torsion of Multiconnected Thin-Walled Cylinders,” by Frank Baron, Journal of 

‘Applied Mechanics, June 1942, pp. A72-A74, 
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iB Bredt’ > formula is ‘based on Hooke’ 8 Law and on n the requirement | of ¢ con- 


in which A, is the area within a closed circuit and a is a measure of the angle 


TABLE 3 _-RELATIVE VALUES FIANDJ 


Sec. 01 


Values of Values of Values of _ 


Towers 


e's 


+5 


> > 


First 
Narrows 
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13 


I 73 Values of Cy Values of a 


‘First Tacoma, . 
2 


_ Application of Bredt’s formula to each cell of the given section (Fig. 9) 
but two independent equations, namely 4 


t...... 


— 

Because of the symmetry of the section, T; = Tg and — 
2 
y Of delormations along aclosedpatnsuchas tne peripnery Olacell, 
4 
| 
| | @ | ® | © | mM |] & (9) 
— | 7.99 | 2.58 | 6.71 | 0.427 | 1.11 | 0.429 | 111 
73 “9.102 183 
34= 20.2 = = = J = =0, 
34SC, S0,263| 1.35S*S2.40| 0.475Scy $0,485] 0.66624 
Walt Whitman|0.331 $C, $0.347| 1.22331.31 | 0.345Sc, $0,358 0.53452 S0.654 
— 
3 


| Values of 


‘Golden -0.309 | 0.484 | 0.258 | 0. 7 
9,006 70.080 | =-0,040 | 0,864 
Narrows 


0,245 | 0,368 | 0.316 | 0.633 


Solving these ¢ equations for and T,, 


= 
(14) 


The determination of J because ecient complex as the number of | 
t cells in a section increases, because of the Senne increase inthe | 
= of simultaneous equations that must be solved. 
Estimating the Properties of Other Towers.- were obtained 
_ between A, I, and J for the various sections of the towers. The relationships | - 
can be useful in estimating the properties of other towers and are summarized ys » —_— 
in Tables 3 and 4. In Table 3, the representative sections shown are the same > 


For the various sections, the values of I can be estimated by means of the 


— 
| 
| 
| 
Sad 
in which the angle of twist per unit 
Consequently, the value ¢ of J for the section is given by 
1 
a 


"BRIDGE TOWERS ? 


= 0) (BASED On ACTUAL ‘MOMENTS OF OF INERTIAS) et 


| 


0,945 


0. 


(0.428 


' ‘in which c isa coefficient and ddenotes the depth of the section — to the 
axis of 1 The _ coefficient c depends on the distribution of the material in the y 
section. Its value is between 0 and 1/2 and it is practically an invariant for a 
om any of section.!1 In Table 4, ranges of values are given for cy of the - 
legs an Cy of the cellular struts. Note that the | values of Cx are about the 
strong axes of the legs, 
a - The values of J for the cellular sections can be estimated ~— means of the 
= ratios between ss (or ly ) and J. The ranges of values for 1, /J of the legs rte 
I y/J of the cellular struts are givenin Table4. 
a = An additional aid is given in Fig. 10 for estimating values of c and VJ of 
- multi- -cellular thin-walled sections. The cross-section shown in Fig. 10 re 
rectangular. It is composed of a single hollow cell, and has a constant wall 
thickness t which is small relative to the lateral dimensions b and d, In esti- 


SS Can n be used by substituting an average rectangular c cross- s-section for the multi- ‘ 


a. Torsional Behaviors of the Towers.—Four cases of each tower are 4 


ss the condition of a torque about the y-axis with the vertical. saddle reactions | 
_ 11 “Structural Design in Steel,” by T. C. Shedd, John Wiley and Sons, Inc., New — 


 &§ 0,968 | 0.016 [ 2.05 | 0,234 | 0,383 
| 
141 
— 
720 | 
0,214 | +#0.473 | 0.756 | 0.878 | 1... 

7 

| be stated in the same way as for sections that are ordinarily used for columns; 

4 
4 
hy! 
— 

ag 
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‘TABLE 6. —SUMMARY OF DISPLACE) EMENTS AND TORSIONAL- RIGIDITIES 


Golden 
| 1.98 | 459.0 


Tacoma 


= 


equal zero, pan cases are designated by A, B, € and 'D. Case A is the ye actual 
tower (Figs. 6, 7, and 8). In Cases B, C, and D, the legs are identical to those = 


of Case A but the struts are different from those of the actual tower. In Case 


he B, the struts are the same in number as in Case A but the values of I, and . : 
are infinite. In Case C, the tower is braced with only the top strut of the eee 7 
tower. In Case D, the tower is braced with no strut. APRA BOP a. — 
ies: The results of the analyses are summarized in Figs. 11 to 13 and Tables © 

| 5 and 6. . The distributions of m;/My, my/My, and b vz/M. , along the legs of 

the towers are indicated in Figs. 11 to 13. Representative values of 


7 same ratios are listed in Table 5. The representative values of my /My, 


—mty/My, and b vz/M, are for sections at the ends of the top and bottom © . 
- struts, and for sections at the : tops and bottoms of the legs (Fig. 1 for the lo- ; 
— cations of sections rq and 01). The values of the moments and shears in _ 7 


Z tower are dependent o on the relative magnitudes, not on | the absolute — _ 


an depend on the over-all dimensions of the towers and on the absolute values of 
E I andGJ ofthe various members of the towers. The of the deflections 
: and chord rotations in Table 6 are for values of M, = 108 by 10° in. Ib, 
E = 30 x 10° Ib in.-2 , and E/G = 2.6. ‘Further, the values of the torsional 


canta are independent of the values of My, because the torsional rigidity 


* mements of inertia in bending and torsion for each strut and tower leg seg- 


would eumeed that aad is assumption is valid for the Golden Gate towers, but 


— 

towers L | Az, in| in 1012. A,. in 
L | | Rad. | Ibin.? | in, | Rad 

4 0.104 | 0,042 “0.078 11,7000 

0255 | 1.45 | 6.21 | 92.2 

(Walt 3,33 | 0.774 | 1.44 | 318.0 0.477 | 0,323 | 0,600 160 

— 
ag 
_—— =8§}=©in Table 6, a summary is given of the deflections and chord rotations at 
the tane af the tnawer in rigidities of + towers factor 

— 
— 


‘(BASED ON ACTUAL |! MOMENTS ¢ 


Red... 
(13) 5 


381.0 | 0,086 8.32 | 109.0 


0.4 426 
0,426 


"questionable for the first Tacoma Narrows and Walt Whitman towers. Further-_ 
_ more, the effect of the saddle reaction at the tower top on the longitudinal bend- 
= of the tower legs is neglected inthe simplified analysis given elsewhere.? | 
4 The analysis ‘was modified to include these factors and a study was made to 
- determine the degree of approximation that is inherent in the assumptions of 
the simplified analysis. Tables” 7 and 8 show the relative values for several — 
ame compared to those given in Tables 5 and 6. The relative values of Table 
7 are for Case | A of each tower when the saddle reactions, P, are equal to 
"zero, These values are based on the assumption of average but uniform mo- 7 
- ments of inertia in torsion and bending for each strut and tower leg | segment be 2 
‘The relative values in Table 8 are for Cases C and D of each tower with 
the saddle reactions P = 0. 4 Pe, and for Case A of the Walt Whitman tower 
with P = = Ppt. The actual moments of inertias are considered in Tables 5, 6, ; 
and 8. _ Table 8 shows (a) the influence of the vertical saddle reactions is ap- 
‘preciable for ‘Case D (when only the longitudinal bending of the tower legs is” 


4 


on the torsional behaviors of the actual towers is practically negligible. a 

> Although | the results of Table 8 8 are based on a deflection theory, superposition © 
=. provided that the values of P are held constant. 120 

_ Tables 5, 6, and 8 also show the differences that exist in the torsional and 
| rigidities of each tower A and Note & 


== 


FOR A ACTUAL TOWERS» 


> 


= 
4 —— - 
— | . | Wy in| Stow,] ay, | | ¥y,in |Gitow, | 
SS in} 10-3 | 1012” | in in 1012 
0025 
| 1,60 7,69] 593 | 008 
ug — 
— 
g 
and Pe for these towers. The flexural rigidities of these Sowers 
with P = Ppy are almost nil. 
7 a, the hypothetical towers. The reader may wish tO compare at a later Stage the © e 
“Theory of Elastic Stability,” by S. Timoshenko, McGraw-Hill Book Co.,New York, 


August, 19681 
results” of both groups of towers, selecting in each instance an pypothetical 
tower that is considered to be the most representative of an actual tower. 2 
ae _ Statical Relationships. —For Case A, the values of mx01/My are greatest 
for- the Golden Gate towers (1.32) and smallest for the first Paseo Narrows 
_ towers (0.836). These values are principally dependent onthe values of a; and > a 


are in agreement with those obtained for the hypothetical t towers, To ot 


which have ‘two or more struts and which have values of aj and ag that are 
approximately the same as those of the actual towers. The values of aj are 
* 99 and 3.90 for the base sections of the Golden Gate and first Tacoma Se 7 
rows towers , respectively. The average values of ag are approximately 1. 7 
and 1.21 for the same towers, respectively. Now obtain the values of myg;/My 
for” the corresponding hypothetical towers. The latter values should check 
ithin reason the values for the actual towers. fe 
7.—RESULTS ON UNIFORM 2M MOMENTS OF INERTIA; 


"Relative Values, Values In Tables 5 5 and 6 As 1,00 100 


af 


Ts illustrates the procedure to be followed in correlating the results of 
the hypothetical and the actual towers. It does not imply that complete agree-_ 
7 ment will be Obtained in each instance between the two sets of results, 
‘“ In Figs. 11 to 13, it will be observed that the distributions of the bending 
re ‘moments m;/ My along the legs of the towers are uniform for the Golden Gate 
uniform for the first Tacoma and Walt Whitman — 
towers. The principal causes for these differences are the number and kinds» e 
of struts that are used in each tower. For example, there are six flexible 4 
struts for the Golden Gate towers, whereas, there are two stiff struts for he 
Walt Whitman . Also note the influences of the and the values: 


of GJ on n the bending moments: along the legs « of the towers. The latter in- = 


"fluences a are for the first Tacoma Narrows towers. 
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BRIDGE TOWERS» 
Z For the struts of the various towers, the values of the bending moments ~ 


 my/My) and the twisting moments (m,/My)s are biggest for the Walt 
Whitman towers and smallest for the Golden Gate towers. This is principally 7 
_ because of the differences in the relative stiffnesses of the struts and the legs 
2 of the towers. For each tower the biggest values of (my/My)s and (mx/My )s 
occur in the top strut of the tower. This is most apparent for the top al of 
the ‘Walt Whitman and the first Tacoma Narrows tor wers. (it is in 


TABLE 8. —RELATIVE EFFECTS OF VERTICAL SADDLE 


(BASED ON ACTUAL MOMENTS OF INERTIAS) 


Golden Gate First Tacoma Narrows | Walt Whitman 
. 

Seen 33,86 x 10° lb .. x x 10 lb 


Pe 


Relative Values, Assuming Values in Tables S5and6As 1.00 


this analysis that the bottom strut adequate t to resist the t bending» 
and torsional deformations at the base of the tower 
= gl _ The values of b Vz/My and 2 my/My are of particular interest for the top " 
ok and bottom tiers | Of the towers. Because the values of 2 my/M are related to 
those of b vz/My, the succeeding study is restricted to the values of b vz/My. — 7 
‘This value indicates for a panel the proportion of the total torque | which is 
ele by the transverse shears in the legs of the panel, For the top panels 


of the three towers, values of b are for the 


— 
«| 
— 
| 
— 
| 
04P, | 55.2 2,96x10°lb | 25,5x 106 Ib 
| q 024500 P9640 0,288 
Cor 
~ P= | P= | P= | P=| | P=] P a 
O44 PL] 04 PL] | 0.4 PL | 04 PL 
— |... | aoe]... [ 203 aoe 
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‘towers (48.4%) and smallest for the first Tacoma Narrows towers (almost 0%). 
_ These values are in agreement with those of a set of hypothetical towers 
which have a single strut and values of ay through h ag that that are about the same 
as those of the top panels of the towers. 
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Per 


—MOMENTS AND SHEARS IN GOLDEN 1 GATE BRIDGE TOWERS — 


and 0.184 for the Golden Gate towers. These values are decreased by 3 an in 

_ crease in the number of struts and, to a lesser extent, by a decrease in the | 

sizes of struts. The ee are e obtained for the — 
towers. 
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FIG. 12.~MOMENTS AND SHEARS IN FIRST TACOMA 
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ngle of twist a at t the: top ofa 


i 


2 in which A zis the displacement at the 2 top of a saddle and b is the width between | 
= cables. “The values ofAz, Vy y, and G Jtow for the actual towers are summarized 
-. Table 6. This table owe that fora value of My = 108 by 106 in. lb the dis- — 

placement Az of the Golden Gate tower is 1.07 in., and the corresponding dis- 
_ placement for the Walt Whitman and first Tacoma Narrows towers are 72% and 
477% of this value, respectively. This comparison does” not show all of the 
characteristics of the towers subjected to a torque, because of the es differ- 
ences in the heights, widths, and makeups of the towers. tae 
‘Further comparisons ‘of the three towers are given in Table 6, in which 
values of G Jtow and G dow/(E Iz,) are listed. The value of E Iy, at the base 
of each tower is selected as a reference. The values of G ~. /(E iL) for the 


‘tow 


at Case B, the towers are e the same as ; in Case A, except the struts of Case 
are infinite. The results of Case B are at complete variance with those of 7 
‘The statical relationships for the towers of Case c Ceili differ- 
ent from those of Case A. The torsional deformations, however, are about the 
same for the corresponding towers of Cases Aand C. Therefore, the torsional 
rigidity of a tower can be approximated as if it were braced with its top strut 
~ only. This approximation may be surprising to some readers. It is possible | 
sy because of the similarities in the distributions of "the a along the legs of the : 
s 


F towers for Cases A and C, For - example, compare the distributions of ‘m,/My 


Nua for Case C closely fo follows the corresponding line for Case . 


a more e detailed study is given of the influence of different 
factors on the torsional behaviors of suspension bridge towers. The analysis _ 
fora group of hypothetical towers in which the and 


q ie Description of the Towers.—An outline of an hypothetical tower is shown Pas 
in Fig. 1(b). The tower is comprised of two parallel legs, r struts, and is sub- ‘a b 
jected to a torque The legs are prismatic and the values of E I; and G Jy, + 
are constant. Allstruts are identical, equally ‘spaced, and have values of E I, At 
and G Jg which are constant. The struts divide the tower into r panels, each > 


— 
and the torsional rigigity of a towe 

| 
the 

hip is obtained for a similar group of hypothetical towers. However, 4 
for the two groups are not the same because 
absolute values of G for 
— 

| 

J 
— 
7 


| 


and by by the e absolute \ values of hg and E ly. i In this study, | dimensiones ratios 
are used. These ratios are useful in design, because they define (1) the function 

of the lateral bracing in resisting the torsional deformations of a tower, and Ea 
(2) the vennee of changes i inthe number of struts in oe tower, the values of 
me The following sets of values are considered forthe dimensionless ratios of a 


= 0.1; 
‘a, = 1.0; 2.0; 3.0; 4.0 
40... 


‘Mmiting value in mind, the values | of the r and a terms in Eqs. 22 to 26 at 


the characteristics of 1,360 different kinds of hypothetical towers. Sire ie 
_ Torsional Behavior of the Towers .—The torsional behavior of the hypothet- 
ical towers was determined by means of the analytical procedure given else- 
5 
where. A computer program was written for and com- 
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influences of the. ditterent variables on 1 the values of 


‘summarized in Figs” 14 to 21. In these figures, values of 
and a; and aq are considered. Furthermore, the values of 2 my ili and — 
i, = b vz: :01/My fora iret tower are shown in the same figure, because 


a requirement of statics for each tower. 
. The results of the computations are grouped in the figures to demonstrate © 


the influence of a change ineach variable on the torsional behavior | of a tower. 
_ In each figure, three of the ave variables (r_. and aj to aq) are held constant 4 


_ Figs. 14 and 15 show the effect of eiabioaiieta of ay and r, for two groups e 
_ towers, E and F. In Group E, the towers have a height-to- width ratio of 4:1; 
and in Group F a height-to-width ratio of 8:1. The values of ag and ag — 
the same for ¢ each group and are constant. (This does not mean that the values 
_ of the bending and torsional stiffnesses of the struts are a constant when a, is 


i the values of aj; whereas, the values of, G Jtow/(E I) are decreased; 
_ (2) The values of mx, -01/ My and bvz;01 /My are increased by an increase in 
the value of r; whereas, the values of G Jtow/(E are not 


(3) The values of mx;91/ /M and G Jtoy,/ (Elz) ar are influenced 
 Ttis interesting to compare the torsional behaviors of towers which have r | 
struts” to those of Similar towers which have no struts. For a tower without 
struts, the values of the moment and shear at of the tower, and the 
torsional rigidity, are given by 


— 
— 
— 
= 
— 
| 
oa’ ‘e 


4 4, 
wh) 
no struts 


with those 

0 Figs. 14 to 21, a major change occurs in the behavior of a tower as soon as 
as strut of een ae is introduced at the top of the tower. For ex- 
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ample, the torsional rigidity of a a tower of Group E in Fig. 15 is from 4 to 35 
times greater than that of a similar tower which has no struts. For a tower _ —_ 


> _ of Group F, the torsional rigidity is from 12 to 100 times greater than that of 
a similar tower which has NO struts. Inthese and other examples, the increases: 


i of aj. The values of mx;01/My and b vz; ;01/ My are appreciably reduced in a 
_ tower as soonasthetops strut is introduced. A reduction in the value of b vz;01 | 
is accompanied by an increase in the value of 2 m -01/ oo — 
Briefly, the major change in he mgngereen behavior of a tower isa a result 


in which tx and A¢yare angle- changes about the x- y-axes 
i“ The angle- changes about the two axes are coupled by means of Eq. 35, and an 
. « interaction occurs between the torsional and flexural deformations of the legs 


sonata oscillations ‘of the b and of Further 
should be given to the influence of the torsional rigidity of a tower on the re- 
_ gulting deflections of a bridge and on the response to torsional oscilliations, 

It appears (Part a)4 that the fundamental torsional mode of a bridge is largely 
prevented by the resistance at the towers to a twist about their vertical axes. 
a. The following quotation (Part a)4 is of interest: “The towers, however, offer 

"much greater resistance to twisting than to bending parallel to the bridge axis 
5 and contribute therefore a material resistance to torsional oscilliations of the 

In Figs. 16 and the values” of a2 and ag are varied. The influences ot 

_ these variations on the behavior of a tower are illustrated for two groups of 

towers designated as G and H. These groups are similar in certain respects 

to the two which are considered in ‘Figs. 14 and 15. In Gre p G, the towers 
fo two struts each and a height-to-width ratio of 4:1, and in Group H ~~ 
ie towers have four struts each and a height-to-width ratio of 8:1. For both oth 


Ih Figs. 16 ae hd the effects of increases in the values of ag and ag are h 


The values of 1 


b vz- -01/My and GJtow/ (E Ir, are e slightly increased. The latter increases ; are 


biggest for towers with small values of aj, 
5A (2) The values of the statical relationships a and the torsional rigidities of i) 
the towers are almost the same for both groups of 


The « comparisons which are dealt with in Figs. 18 and 19 are 2 of a \ different - 


category from those dealt with in Figs. 14to 17. In mel 18 and 19, the towers 7 


¥ | 
= 
- 
| 
— 

9 
= 
and of the strut. Similar coupling or interaction effects occur in other kinds 
of structural systems in which bracing is introduced. 
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FIG. 16. VARIATIONS. IN IN STIFFNESSES OF STRUTS 
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FIG. . 19.—TORSIONAL RIGIDITIES OF GROUPS J AND K 


fluences of these on ‘the torsional behaviors of the towers are 
_ trated for two groups of towers, namely, Groups J and K, In | Group J, J, final value ‘alue 7 
of ag) is 1.0, and in Group K the value of agis2.0. — —_—s 


To 
are radially different in m and 
of h,g/b are varied as the values Or a9, 4g, = a4 are stant, ine in- 
q 
— 
- 


oles Figs. 18 and 19, ‘the e effect of changes in follows for of 


if The values of mx. 01/My of are increased by an increase 


= For a given» value of the values of mx:01/ M b Vz; ;01/My, 
2 my :01/My, and G Stow/(E Iz,) are generally different for the tower with one 
4 strut than they are for the towers with twoor more struts. 

The relationships in Figs. 18 and 19 can be used in comparing the values" 

of the bending moments, twisting moments and shears at interior - sections of a 
given tower. _ For example, consider the values of mx.jx /My, b vz k/My, 
a and 2 m ik / My at an arbitrary section jk of the general hypothetical tower — 
e (Fig. 1(b)). Asgume that the values of a, through aq are held constant. Figs. 
18 and 19 indicate that the values of my, jk /M >» D Vg jk /My, and 2 my. i, /My 
practically the same for all sections whic h are or more paneis re- 


moved from the top of the tower. . The results of other computations support - 


Figs. 18 and 19 indicate the need to examine more closely the torsional — 
behaviors of towers with single struts 
_ The torsional behaviors of towers: with single struts are illustrated in 
: ‘Figs. 20 and 21. They are illustrated for a group of towers (Group L) in 
which the values of a; and ag are varied and the values of ag and ag are — | 
constant. The same kind of information is givenin these figures as in the pre-— 
_ ceding figures. Additional information, however can be directly deduced anil 


for the shears in the top struts. The moments —= the ends a the top struts | a 


the same as the corresponding moments at the tops of the legs. _ satiate. 
- Only the ‘most important observations will be stated for towers — are 


similar those of Group L. They are as follows: 


ef (2) In general, the values of G Stow/(E Ip) a are dnenennet by an n increase — 
in the values" of ay. ‘For a constant value of a4, th the values of G Stow/(E 


+ 
i 
| 
J @ | 
q (1) Practically every variable (a; through aq) has an important influence 
— 


hey are dependent on 1 the value of ag. Finally, the values of ao ‘and. ag oo 


m 


VALUES OF bv 


2 


little influence on the value of G Stow/(E I;,), excep — for small values of a; ay. 
The reader is reminded that all values of aj through wh ag are assumed to be 
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ry 
—y 1 as a 
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Statement (1) also holds for the ye panelof a a multi-strutted tower. In fact, 


fixed at its base. 


The preceding study illustrates the of different factors, 
and the approximations which can be made in the torsional be- 


of vertical on towers. 
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‘By Jerome M. Raph ael 


fraction of the total stress in the Existing experimental | 


ods that have been used to determine dead-load stresses in structural models 

are reviewed and evaluated. A new experimental technique termed the “meth- 

od of integration” is described. The theoretical background of the test is de- | 
scribed, the similitude factors are derived, description is given of a typical — 
experiment, and types of results that may be expected are given. The method — vl 
can be used to study alternative tive construction procedures to a in desired ae 


of dead-load stress. 


on weight ons of the most important loadings on the structure. 
Measurements of massive structures in service2 have shown that stresses due 
to dead weight establish the predominating stress pattern, and that live- ¥ 


_ stresses, due to water loading and temperature changes, modify this stress 
_ pattern, Furthermore, observations of actual structures have shown that con- 
- struction procedures | have : a considerable effect on establishing the pattern of 
a. Note,—Discussion open until February 1, 1962, To extend the closing date one month, | 
a written request must be filed with the Executive Secretary, ASCE, This paper is part 
of the copyrighted Journal of the Structural Division, Proceedings of | the American So- | 
of Civil ‘Vol. 87, No. ST 6, August, 1961. 
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DEAD-LOAD STRESS IN MODEL DAMS BY METHOD OF INTEGRATION a 7 

7 

— 

asta Dam,” by Jerome M. Raphael, Trans-— 


Hence, when it was decided to perform an stress analysis of 
_a large and unusually shaped concrete dam, one of the prime objectives of the ~ 
investigation was the ; development of a procedure for the determination of the a ” 


dead-load stresses during construction and on completion of the structure. Al 


of dead-load Stresses in completed structures, none showed the development 
of stresses in the dam taking into account ‘construction procedures. it is be- 
lieved that the method of integration developed to show the gradual accumulation 


atte] CM 


| } 


ELEVATION 


nee FIG. 1.— -MODEL DAM SHOWING POSITIONS OF TUBES / AND LINES FOR 


“4 “of stress has many points of superiority over other experimental methods, Be- 

fore describing the method of integration, a brief historical sketch will be given | 


of the development of experimental aeenaueed for the determination of dead- — 
stresses. 


— 


METHOD OF APPLIED LOADS. 


— 
x= 
= 
‘The method o rimental 
is stress analyses of dams. In the classic J. S. Wilson and William Gore exper- * 


ST 6 -DEAD-LOAD STRESS 


g iments3 in Great Britain, , reported in 1908, a rubber slab representing a a typi- 


; cal section of a gravity dam together with part of its foundation was tested, as a 
e shown in Fig. 1. After dividing the cross section into a number of incremental 


areas, pins driven through the model at the center of mass of each subdivision — 
were loaded with dead weights in proportion to the dead weight of eachincre- | 
ment. Strains measured on the model | before and after loading were used t a 

A variation of this method is weed at the Istituto Sperimentale | Strutture e a 
Modelli4 (ISMES) at Bergamo, Italy. Here, as shown in Fig. 2, the load is ap- 
_ plied internally through rods anchored in the mass of the structure. The rods, 7 

_ in turn, are loaded either with dead weights, or tied into a mechanical system 
= with hydraulic jacks, shown in Fig. 3. If the profile of the dam model i 


thus loaded is not too irregular, this internal loading system has the advantage 
of leaving the surface free for the application of strain gages. However, some 
_ local stress concentrations are bound to exist in the immediate vicinity of the | 
_ anchoring devices and — must be taken to avoid strain measurements at this 


METHOD OF 


— 


The use of the method of inversion for the experimental determination of : 
-dead- load stresses in essentially two-dimensional models, has been described — 
by the Laboratorio Nacional Engenharia Civil (LNEC). 5 In this method, a 
model of the structure is fabricated and fitted with strain gages. The difference © 
is strains between the model in its upright and inverted position is measured. : ; 


__ By inverting the model, a loading effect is produced which is equal to twice be 
that of gravity. _ The sensitivity of this test can be increased by increasing the 
density and decreasing the modulus of elasticity of the model material. . Fig. 40 
_ shows the test setup used at LNEC on a model of abuttress dam. Considerable 
_ difficulty is experienced in attempting to utilize this method in a three- -dimen- — é 
sional model of any size or complexity. Furthermore, unless the foundation is 
s _ made so rigid that deformations are minimized on inverting the model, strains — 
will be induced in a three-dimensional model, from the deformed repeal 


which are unrelated to the normal gravity stresses. This might be a special 
weakness for utilizing this method in the testing of ar arch dams, because arches 


are quite sensitive to abutment 


-METHOD OF IMMERSION 
A. Biot has the analogy between the two- ‘state of 
stress” resulting from body- -force loadings, such as the earth’s gravitational 


~ 


“Stresses in An Experimental Investigation by of India-Rubber Mod-— 


els,” by J. S. Wilson and William Gore, Proceedings, I.C.E., Paper No. 3705, Vol. 172, 7 


Materiaux pour Mod@les Reduit et Installations de Charge,” by E. Fumegeitt, 
ISMES Publication No. 13, Report presented at the Internatl. Colloquium on Structural 7 


5 “Analysis of Concrete Dams by Model Tests,” by M. Rocha and J. L. Serafin, »Com- 


26> = 


munication No. C36, 5th Congress on Large Dams, Paris, 1955. = 
6 “Distributed Gravity and Temperature Loading in Two- Dimensional Elasticity Re-— “4 
placed by Pressures nad Dislocations,” ” Iby M. AL Biot, Journal of 
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3.—JACKS AN AND LINKAGE FOR APPLYING 
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~ BEING TESTED AT LNEC _— 
— 
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ESS 
DEAD- -LOAD STRESS 


.-* is inverted and lowered into a 1 pool of mercury ‘whose upper surface just 
touches the ground surface, as indicated in Fig. 5. Stresses at various points 
in the slice are then related to corresponding stresses in the prototype dam by 
ie use « of Biot’s similitude equations, which take into account the scale factor and _ 

the relative densities of the immersing fluid and the materialof the dam. This | 
7 _ method ha has been weed in LNEC | - Lisbon for finding the stress distribution in 


— 


FOUNDATION 


= 


2 


5,-MODEL DAM IMMERSED IN MERCURY. 


a section of a gravity or of a buttress dam. ‘Similar — has 


for the three- dimensional state of stres 
_ The principle of the centrifugal method for the cupetateh determination 
_ of dead-load stresses is the substitution of centrifugal forces for gravitational - 
- forces by rotating the model in a centrifuge.? This method has been used with © 
 ——_ success in the Engineering Materials Laboratory of the University 
of California by H. D. Eberhart in two-dimensional photoelastic models of the 
cross sections of several gravity dams.8 The model, made of a sheet of photo- 
; gas material, was rotated at high speeds at elevated temperatures (200° F) 
several hours un until it reached thermal equilibrium. While maintaining con- 
4s The Stress Function and Photo-elasticity Applied to Dams,” by J. H. A. alan, 
7 8 *Analys sis of Gravity Load Stresses by. Photoelastic Methods,” * by J. J. Polivka, 
Proceedings, 16th Semi-Annual Eastern Photoelastic Conf., 1942. AY 


— 


op ary loads. Inthis meth- _ 
| 
not been made 


196 


room temperature was reached. "When ‘the model was removed, it had frozen 
7 within it a stress pattern resulting from the equivalent of an increased gravi- — 
 —— field which could then be analyzed by photoelastic methods. pes => ae ; 
In another. experiment, 9 A. J. Durelli, went a little further toward three- 
= analysis. His experiment was concerned with the distribution “~ ur. 
dead-load stresses ina single buttress of a massive head buttress =§ 
model was cast using a photoelastic material of a single buttress together v with 
part of its foundation. As soon asthis casting had solidified, it was transferred 
to a centrifuge and rotated rapidly while curing progressed. At the end of the — 7 


3 model was s then cut into two-dimensional slices, polished, and the ‘stresses at 


various planes determined photoelastically. 
 Bothof these test methods are excellent solutions for essentially two-dimen- — 
problems. However, certain limitations must be recognized when 
“tending this method to three- ‘dimensional problems. The most serious aifti- a 
wide disparity means that ‘three- dimensional models pany of. plastic do 
not deform in all respects like the prototype, and that, as a consequence, the _ 
_ stresses are not distributed in a plastic model like those in the prototype. » 
second and less stringent limitation is that unlike gravity, centrifugal force is 
> not constant throughout the model, but varies as the square of the radius from 7 
_ the center of rotation. Only by rotating a fairly small model on along arm can | 
. the variation inthe » centrifugal field be held within feasible limits. A third lim- 4 
_itation is that of rma _ Both centrifuges mentioned were designed to handle _ 
fairly small models. ‘Difficulties can be anticipated in testing a model of a 
complex shape in sizes that can be handled by centrifuges. aviation and 
: missile industry has developed large centrifuges that make it possible to test 
_ fairly large plaster models having the proper Possion’s ratio for twree- ~dimen- 
sional analysis. However, ‘the cost of equipment for such a test makes it 
feasible for the average experimental stress analysis. 
wee. 


HOD OF AGE: 
Per 
A significant number of large concrete dams are constructed in stages, ‘o 
balancing the capacity of the reservoir with increasing hydrologic require- 


ments. In this way construction costs are spread over a longer period of in- 


dams are constructed by stages, : since in each case definite controllable vol-. J 
umes of concrete are added as increments during construction. 
- _ The method of construction by stages was devised at the LNEC in studying s* 
_ how to increase the height of existing dams where it was required toknow with 
-_ precision the effect of the weight of the dam taking into account the construction © 


ae 9 “Experimental Stress Analysis of the Bekhme Dam,” by A. J. Durelli, Armour Re- 7 
search Final Report, Issued for Harza Engineering Co., December, 


3 creasing demand rather than — concentrated at the beginning of the period, — 
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DEAD-LOAD STRESS 


1 sequence. J. J. L. Serafim and J. P. da Costa described10 the philosophy of this — 

“The various stages of of a dam correspond to an equiva- 


al lent number of differently shaped soils and consequently the final state 
- of stress in the solid will be the sum of all the states of stress resulting ices 
_ from each stage. This final state of stress is not equal to that which - Bote. 


= with its own n weight in o in one 2 operation; the greater the number of ¢ connec- an 
- tions in the dam during the construction, the greater will the difference - 


Proceeding from this on an experimental method was devised 


reversing the order of construction of a dam. _A model was constructed of a 
plastic and fitted with strain ‘gages. ‘The model was then systematically cut 


down from the top by stages in a reverse order of the proposed construction 
sequence. - Readings of the strain gages reflect the removal of the dead weight, - a 


and, by reversing this record, the gradual development of stress as construc-_ -_ 


tion proceeds can be followed easily. It can be seen that even with fairly high 
density model materials, the sensitivity of thistype of test is not great because ~ 
" of the low general dead-load stress level ina model. It was to increase the al 


4 sensitivity of — -load tests that the method of integration eee herein 


id 


Four of the five me methods for dead- load testing described previously — 
plied loads, , inversion, immersion, and centrifuging - had an - important com- 
mon defect. Any test in which a completely constructed model is suddenly en- 

_ dowed in its entirety with weight does not in the slightest degree reflect the 
_ incremental manner in which weight is added to the prototype structure, and» 


the: manner in which the structure gains its strength, and leads to ot 


struction proceeds. The fifth method, the method of construction by stages, 

lacks the requisite sensitivity for accurate analysis. 

_ The importance of the construction program in determining the distribution | 7 

= é stresses in a concrete dam has been demonstrated. 2 In studying a proposed . 

_ design for a dam, question ha had been raised as to the behavior of the structure — - 

because of particularly large overhangs. Hence, a new experimental method 


| 
| 
» 


was devised by which the stresses could be determined not only inthe complete 
structure, but also at various stages of partial completion, following any given 
construction program. This method has been termed the “method of integra- 
As finally developed, in the method of integration for the determination of 
_ dead-load stresses, the model is cut down by increments leaving a cet 
surface representing the top of a lift. On this ourtace , a uniform load is 7 
plied, 
“Methods and Materials for the Study of the Weight Stresses in Dams by 


of Models,” by J. L. Serafin and J. P. da Costa, RILEM mnenee. deteniieastenin on Models 
of Structures, Madrid, June, 1959. 


conclusions as to the development of dead-load stresses in the dam as con- q 
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a es » ietened by adding 1 ft of concrete at s successive > elevations. This is ac con- 
a tinuous function which can be integrated graphically or numerically to indicate 
_ the effect of successive stages of construction, including the dead load stresses” 7 
4 in the complete dam. Thus, critical stages in the construction when dead-load f 
_ stresses area ‘maximum which exceed the stress state in the complete dam ° 
at be indicated by a of a ae of stress versus height o! of | dam. Py 


The general ‘considerations that apply to. the ‘derivation on of the similitude 
- factors used inthe method of integration for determining | dead stresses indams 
af 


STRESS, _ 
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Yo 


om 


toa final elevation H. At an intermediate stage of construction, it will have — a 
reached elevation h. In testing, the model dam is cut down by stages to some a 
elevation y. The similitude problem és to find the effect of the increment of 
height dy on the strains and stresses at location Xo YO | 
For a geometrically similar model, the ratio of model length 
type length Lp is defined as 


‘The ratio of the density of of model it loading Ym 

ing Yp is defined as 


wee 


Gages on the model! factors, this is then related to 
See aes, ¢& the effect on the prototype dam of adding 1 ft of concrete at that elevation. : 
— 
— 

| 
— 
— 
— 
— prototype load- 


‘Thus, for a dead- load test of a sae dam amde ot a plaster weighing 40 
4 ‘pef representing a prototype dam made of mass concrete weighing | 155 pef, the 
Expressed in these units, force then ‘becomes L3 and givens which 
is force per unit area, becomes 


T Thus, the relationship of stresses between model and can be 


 simititude factors. for strain € can be derived starting with the elastic ; 
stress-strain 


Consider a load per unit area of wym to be applied at elevation 
dy, then this will apply a stress dy a at y. 


_ If modeland prototype are geometrically similar, and both have linear elas- 


- tie properties, the ratio of the ‘stress or to the model tothe stress applied 
to the prototype will be 


itt ites that the increment of stress in | the prototype due to this incremental 


— 
q 
From Eqs. 3 and 5, one can write we 
3 
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arly for dead- load stress the dam: 


“It will be noted that the model- — factor: wym i is taken within the integral 
because it may vary with each elevation tested, and is thus a function of y. 
ee: - Ina typical test, a model was loaded uniformly with 1,500 lb of lead bricks 
7 applied to an area of 250 sq in. Thus, the increment of stress in the prototype 
to an increment of 1 ft of concrete to the dam was 


Em Y 240,000 x 158 x1 


(13) 


044 €ym (13) 


where | strains are measured in micro- inches perinch | oe 7 

: me In testing, a model is cut down to definite elevations y convenient for anal- 

ysis, and loaded heavily enough to give a well- defined response in the on 


Using the similitude factors contained in Eq. 9 curves of 


cannot be determined experimentally, be- 
, and it is necessary to compute the stresses at - 
the gage | elevation yg due to the 1e effect of loading at that elevation. By Hooke’s S- 
law 


wee 


(14) 
1-p2- 
But 90 is the equal to the density of the concrete of 158 p pef o or 
Pe 10 = per ft of concrete. ff uu is equal to 0.24, Eq. 14 reduces to 
Ep dego = = -1. 04 - 0. 24Epdeg . 


_ Since the direction of this stress is known, 


dego + mere: Ep d €45 ° oe 


— — he prot emental loa 
structed to elevation h is 
— 
— 
do = 
— lc ku 
i 
i 
ag 


| 
For any strain gage ronstie, Bes. 15 and 16 must be satisfied, and all ‘three 
strains must be consistent with the strains developed at adjacent points y. 
_ With the incremental stress function completely defined, the stress at 70. 
; due to construction to height h can be found by graphical or numerical integra-— Z 
tion of the curves representing Eq. 10. 
The principal stresses for any stage of construction can then be calculated 
inthe usual fashion from the stresses along tl the individual directions shown on : 


“dures used totest the 1:200 model shown in Fig. 1. This is the central portion | 
of one design proposed, but not used, for the 750- -ft high, 5,000 ft long Oroville — 
Dam, 11 after conditions at the dam site led to final decision to construct an 
embankment dam instead of the concrete dam on which this study is based. — 
The model was constructed by casting a plaster- -celite mixture in a rigid fiber - 
glass mold on an extensive foundation assembled from prefabricated plaster- 
4 celite slabs. The central portion of the dam consists of an inclinedarch which 


carries on it a pair of riding buttresses. | The riding buttresses not only « carry” 


load into the flanks of the arch, but they stiffen the arch in an indeterminate — 
fashion. It can be seen that in so indeterminate and complex a structure, ex- 7 
"perimental s stress analysis has many advantages over a numerical procedure. 
Ih the gravity loading test, a uniformly distributed load was applied toa 
u - series of horizontal surfaces, after all live-loading tests had been completed. — 
a This uniform load was approximated with lead bricks. At each elevation to be — 
‘4 tested, the area of the exposed surface was measured, and divided by a plani- : 
_ meter into a number of equal smaller areas. The center of gravity of each of 
these small areas was then found, and the position for each lead brick carefully — 
drawn, so that the weights were equally distributed. . Between 50 and 70 lead» 
« bricks, each weighing 26 lb, were used for each load test. Thus the model sup- 
ported nearly a ton of dead weight at various stages of the testing. The cut — 
_ model with its dead weight in place is shown in Fig. 8. When a number of ere 
4 of lead weights had to be used, sheets of paper were interposed between layers | 


to increase the stability of the mass 


sections, cut at even n 100- ft elevations throughout the height ‘of the 
; Three identical tests were conducted at each elevation. Since all the lead bricks 
had to be placed by hand, it was found most expedient to load the model during 
the day, let it stand overnight, and then early the next morning, when all tem-— 
pos had been equalized and the air in the laboratory was still, a reading 
was taken of all strain gages with the dam in the loaded condition, the lead 
removed, and a second series of readings taken of the unloaded model. 
4 Since the model progressively reduced in height during testing, the number of © 
gages r read for any test elevation was likewise reduced. Thus the time neces- - 
sary to read the gages differed at each test elevation. However, using the test — rt 
Ben up shown in Fig. 8, which was designed around the Fairchild RBI (Resis- * 


11 “Structural Model Investigations for Oroville Dam,” by J. M. Inst 
-Engreg. Research, Univ. of California, Series 100, Issue 6, February, 1960. 


— 

| As an example of the use 0 e me Of integration for determining the Jf 

a) 
7 

| 

| 

| 
15: 

— 
— 

| 

— 
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FIG. 8.— DEAD- LOAD TEST BY METHOD OF INTEGRATION AT 


= 
| 


DEAD-LOAD STRESS 
tance Bridge Indicator), 150 strain gages could be read ina period of 10 min. 
_ without difficulty. It was found that model strains were sensitive to the exact 
placement of the uniform load, but an average of three tests tended to eliminate | 
& error. _ Likewise, the use of the unloading cycle for strain measurements 
hd enabled more care tobe taken with the exact eed of each lead brick, and 


‘Fig. 9 shows the and of the principal dead- load stresses 


i the completed structure. It can be seen for that for the intrados, the prin- 
cipal compressive stresses increase gradually from top to bottom and that the 
directions are . parallel to the boundaries upstream and downstream, varying 
7 gradually from face to face. For the extrados, the structural discontinuity in- _ 


troduced by the riding buttress is evidenced in higher compressive § stresses © 
near its lower boundary. Along the outstanding edge of the arch, the maximum ~ 
‘compressive stresses increase fromtop tobottom. The considerable variation 
; of the direction of these maximum principal stresses from the center line of 
- the outstanding leg can be explained in part by twisting, shearing action at this 
io boundary. Fairly high tensile stresses are found near the boundary of the — 
' and plug for both the extrados and intrados. Since these two structures were 
4 bonded together in the model at the arch- plug interface, a structural discon- 
-timuity, proposed for the construction of the prototype would reduce or largely 
eliminate these tensile stresses, 
Considerable speculation was raised as to the behavior of the main arch 
during construction. . A glance at the downstream view of the dam shows the ; 
; convergence of the sides of the arch as construction proceeds upwards, to the q 
- os -— where they overhang the base. In Fig. 10 some light is shed on the grad- ap 
ual development of the principal stresses for the intrados and extrados at the - 
base as construction proceeds upwards. In Fig. 10, the intrados of the dam is 
shown at the left and the extrados at the right. Near the bottom foundation line, 4 
_ circles show the locations for which principal stresses were computed. _Di- 
rectly above each circle, and at elevations corresponding to the top elevation — ; 
of the concrete at any particular time, the drawing shows the magnitude and 
x direction of the principal stresses when concrete had reached that — 
For example, for the farthest downstream location at the intrados, the maxi-— 
< mum compressive principle stress was 135 psi whenthe dam was concreted to 
‘EL. 400, 520 psi when construction had reached El. 700, and 600 psi when the | 
dam was ‘completed. It can be noted that the directions of the principal stress 4 
for the intrados change very little during the construction of the dam, and that 
the increase gradually a as construction ‘progresses. However, 


‘The gradual ¢ development of the vertical stress at the intrados and extrados . 
is given in Fig. 11. Here, superimposed on an outline drawing of the intrados _ 
and extrados, are curves indicating the magnitude and sense of the vertical — 
eee at any stage of construction, for various locations at El. 450 and along» 
- the base. Thus, the extreme left-hand curve at the bottom of the figure indi- 
ten the gradual increase of compressive vertical stress at the downstream 
end . the intrados from 0 to 400 psi as construction of the dam proceeds a 
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_ stress: is compressive. . Fig. 11 also indicates that the compression increases — 

fairly rapidly until construction reaches about El. 700, and after that the in- ie 

crease is very slight. However, for the extrados, the curves indicate that the 
vertical stress is compressive at the upstream locations, tensile at the ceo 
downstream locations, and that it changes from tension to compression nea 
the center of the extrados as construction passes El. 70. 


FIG. —DEVELOPMENT OF V VERTICAL AL STRESS DUE TO 
A slightly different picture is presented for the itiiaiis of | vertical 
stress at El. 450 due to construction of the dam. This information is shown 
on a set of curves at the top of Fig. 11. The effect of the weight of the vertical 
load from the buttress is re readily apparent at El. 450. For this elevation, the 


_ vertical stresses are everywhere compressive for all stages of construction. ’ 
4 The increase in stress is fairly uniform throughout the construction. Looking 


at Fig. . 10 one can see that the vertical load coming down th the a buttress 


—. 

| 
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causes heavy twisting in t the arch below the junction of the ‘riding buttress and 
aoe help make this stress picture more uniform below El. 450. weed 
_. With so much tension shown at the base for the extrados, question might be 
and as to the stability of the arch during construction. In ‘evaluating the sta- 
bility of the base, it was assumed that wherever tensile stresses existed, con- | 
_ erete could crack all the way back to zero stress. It was further assumed that 
_ if the concrete did crack, that the moment and thrust that caused the stresses 
would remain without redistribution to the remainder of the structure, which 
is a conservative assumption for the base. Having made these assumptions, it 
was found that for every cross-section for which stresses were measured, the 
- resultant fell within the base, and it ‘could thus be concluded that the arch was 


; everywhere stable. This was not construed to mean that the method of con- 


struction in horizontal layers was best for this type of dam - another method 
_ involving sloped pours might very well result in ey lessened ananes PE 


CONCLUSIONS 


om A new method of experimental stress analysis, called the method of inte- 7 
A gration, is available for the determination of the dead-load stresses in dams. ) 
, _ This method can be used to find the state of stress ina partially completed — 
dam, and to evaluate the effect of different construction procedures on the dead : 
_ load stress pattern remaining within the structure. This method should have — 
_ considerable utility in predicting the behavior during and after construction of 


_the overhanging portions of f doubly curved, arched 
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GENERAL SOLUTE OF SPACE FRAMEWORKS 


>i The solution of space frameworks with members in any rdieection, , under er any 
‘system of loads, in any direction, is established by applying the slope-defection- _ 

- gyration method. 3 It provides aset of simultaneous equations that can bes ‘solved 

_ The procedure presented herein can be extended to the approximate solu-- 

“tion of plane or warpedsurface structure, concentrating the — of sections 

of the amare in an arbitrary system of generatrices. 


INTRODUCTION 
TRODUCTION 


matical problems, and for. that reason it is possible to solve sets of simultane- 

- ous equations that are established in structuralanalysis. This could be achieved - 
without having to resort to methods of successive approximations that contrib- 

uted so much to the progress of the structural analysis of continuous stru 


tures when the necessary elements to solve large sets of simultaneous equa- 
tions were notavailable; 
“ _ Note.—Discussion open until January 1, 1962. To extend the closing date one month, 
: a written request must be filed with the Executive Secretary, ASCE. This paper is part | 
_ of the copyrighted Journal of the Structural Division, Proceedings of the American aia 
ciety of Civil Engineers, Vol. 87, No. ST 6, August, 1961. A sivicna | Mas 
_ 1 Partner, SACMAG of Puerto Rico, Inc., Cons. Engrs., ais Juan, P cae | 
2 Struct. Asst., Saenz- -Cancio-Martin-Gutierrez, (SACMAG) -arquitectos, 


3 “Orthogonal Gridworks Loaded | Normally t to their Planes,” by I. Martin and J. E. 
_Herngndez, Proceedings. ASCE, Vol. 86, No. ST 1, 
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rk with members having any direction in 
~ under any system of loads can be attacked d by establishing the senmeal al 


179 


t 


ment, torsional moment, and ‘shear in the bar anc the 
equations of forces and moments in each joint. + 
red _ Any warped surface under any system of loads, can be approximately solved 


— 


SPACE 
the tributary surface is concentrated, m: making the e surface structure a frame-_ 


DERIVATION OF THE SLOPE DEF LECTION- GYRATION- 

‘Notation. .—The letter symbols adopted for use in t this paper are defined dand- 

for convenience of reference, in Appendix I. 


BE coincide, as shown in Fig. will be An arbitrary cartesian 
cartesian 
system m of coordinates will be used. 
The direction of each bar is defined by its cosines a, b, andc with : 
veapest to the system of cartesian coordinates selected. = | 
m Each joint is subjected to two types of deformations such as angular defor- ; 


_ mation and displacement. Both deformations can be considered as vectors, 


-- 
= 


tm ox 
a ‘Tt a member AB shown in Fig. 2 is considered, three axes can a he tattnet 


ie 


“the procedure, the direction of the torsionalaxes of the can be 


by following the sequence of the letters assigned tothe joints. Having. done this, { 
_ the direction of the principal (1) and secondary (2) flexural axes can be ‘ 
_ mined for each bar, sothat the three axes T, 1, and 2 form a direction trihedron. 
+The direction cosines of each axis of the bar AB, 


TAB 
| 


— 
“a 
| 
| 
le 
; 


fap + 


pS 

» 


Convention.—Consider bending and , shear forces, 
and also deflections and rotations as vectors. — Hence they all shall be — i 
or negative as they coincide with the — or — direction of axes T, 
2,or X,Y, Z. 
- When | expressing moments in terms of deflections and rotations, these will 
7 be — considered ow and and the moment will show the sign that corre- a 
wi AB can be decomposed into loads ees the three bar axes, and covatennd 
- moments along axis T-T and flexural moments along axes 1-1 and 2-2 can be 


_ The basic slope-deflection-gyration, Eqs. 3, can be applied to | the torsional 
‘moments along axis s T- and the flexural axes: 1-1 and 2- 


MPLAB ‘has 


By 


— 


+ Kopa a 
TAB * 


— 
| 
direction cosines comply with the following two conditions 
> 
Moap = Mpoas * Koap (1 + 82 “2AB Lap) — 
— ) — 
— 


SPACE FRAMEWORES 
- The moments ‘corresponding to any other bar r r may be. obtained by changing 


ie 


Ato Band BtoC inthe preceding equations. 


The following equations can be established f for joints A and B, acco 


PAZ 


+ biaB Pay * C1AB Paz (10a) 
= PAX + Say + C2AB 
°BX + 
41a = 4ax + biaB 4ay 


Soa = Dax + bap 4 “AZ 


sit ‘should be observed that the displacement component 6 7 TB along axis T-T 
- not have any influence inthe torsional or flexural moments of member AB. 
Similar equations can be established for all the other bars that coincide in 


_ If the values obtained in Eqs. 9(a) to 13(b), inclusive, a: are substituted in 
Ee. 6(a) to 8(b), inclusive, the following equations are: 
B Y C2AB 4AZ 
Mipa MF1BA - +B ) 


@2BA 4px + bopa Apy + ABZ 


‘ipa (@1aB + biaB * az) 
+ + bipa + C1BA 


Mropa kopa (1 + 


koBA ax + boap Pay + + CoaB 


* Kopa (apa x* PBA * 


se, 


co" 


| 
= 
cto 
— 9a) 
; 
— 
gq — 
| 
{ 
lm | 
: 
& 
| 


TBA (2TAB + bTAB Say + CTAB (16) 

moments BA> and MTBA@ are projectedalong the three 
axes, the following equations can be obtained: 
= + 82AB M2BA + aTAB - 


MiBA * Mapa * eT, TAB Mrpa 
ae ~ values obtained in Eqs. 14(a), 15(a), and 16(a) are substituted in Eqs. 


17(b), and 17(c), the following ‘equations 


M = +a + 
“BAX ~ [FIBA MF2BA aTBA MPTBA 


rrr 
2 


“TBA TAB “TBA 


“TBA “TBA | 


te 


1BA 2AB 1BA 1 2BA | 


Sax 


> 


7 
a 


ipa ax * 9px) 


* Popa Kopa ‘AX 
* Prpa “ra “TBA TBA 
3X 


fee 
T 


ee 


* 


* Cons rap (1+ 
PoaB * a Pal 
| 1BA “2AB — 


°rpa "ras “TBA 


— 


= 

— 

“BX 

 =séCéddg 

| 

— 

¢ 4 4 

=; 


pa “rea ( 


4 If joint B is in equilibrium with eaten to rotation, it is wegryd that the 
sum of the moments along each axis should be equal to zero. If the effect of 7 


the moments ‘produced by the relative displacements of the bars and the axial 


is equations of the moments are 


= 


7 


- 


th the va 

-19(b), and 19(c) and the procedure is s repeated for all the other bars that coin- 

— in joint B, three > equations as a function of the linear and angular defor- 


tablished, whieh for bar AB 


sel 


(2 


7 


“1BA ~ Tap] 


ugust, 1961 sTé § 
7 B 
— 
Inthe usual case of framedstructures, the moments due to the relative dis- — 
placements of the ba es q 
» 
e) 4 
— 
directions 
— 
» 


_- it the observer faces the plane in which the shears ot the bars develop so 
that the axis perpendicular to that plane points towards him , the positive mo- 
that th produce positive shears at the left end and negative shears at the right 


Moments can be ibe expressed as functions of their — along the three > 


Mpa 1BA Mpax * Mpay * “1p 1BA “BAZ 


* 


= Mpax * Popa Mpay * Copa MBAZ 
i I the values of Eqs. 22(a) to 23(b), inclusive, are substituted in Eqs. 2010) ~) 


and 21(b), asa function of the direction cosines at end B, the following ¢ — 


(24) 


Mpax - + beBA MBAY sana 


a The components of the shear force along the three cartesian axes at end B, j 
BAX ~ "18a * 


“in Eqs. , and the following equations are 


Ba 


\ 
+ 


— 

— 7 

=) 

ln 

and 

— 
: 

— ) are substituted 
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Popa | Von 


bapa 
BA (My, 


_ 


— (M 


22BA - C2BA 


BA 


negative, and the Lap inthe equations for VABX, VaBy> and VBAZ is posi 
Mapy: Mapz, MBAx; MBAY; MBAZ; can be substituted by 
_ values as functions of the components of the linear and angular deformations : 
along the three cartesian axes, found in equations of the type of Eqs. , 11(a) a 


7% —-17(b), and 17(c). The components of shear force along the three cartesian 
= will be expressedas a function of the components of the linear: and angular = 
if joint B is in equilibrium with respect to displacements, it ts necessary 
that the sum of all forces along each axis should be equal to zero, so that the 
components along the three coordinate axes of the shears and the axial load 


along the T-axis of each bar concurring at oe B — os external loads ap-— 


at the joint should be to zero. 


The values” York, and may be obtained changing Ato 0B and 
to C in VABx,; VABY, and VABZ, which are obtainedas 
s.27(a),27(b), and27(c). 


of the Geometrical Equations .— 


— 

| 

| 

— 
7 

With these equations, itis nec- _ q 
mene rometric conditions of the struc- 
— 


FRAMEWORKS 


ture whieh consist in making the difference of the deformations along the T- 
axis of each bar equal to the axial deformation of the same bar. Presse an 


— 


be that the equations are reduced 


Computation of Moments and —Eqs. 19(c), 19(d), 19(e), 
— 28(a 1), 28(b), 28(c), 29(c), and 29(d) are expressed as a function of the com- | 
_ ponents along the three cartesian axes of the linear and angular deformations 2 
- at joint B, and of the axial loads of the bars. Similar equations can be estab- 
_ lished for all the other joints. These equations can be solved simultaneously 
in order to obtain the values of the linear and angular deformations at each | 
4 joint, and of the axial loads of the bars. As soon asthe deformations and axial 
4 forces of the bars ateach jointare obtained, the values of the flexural and tor- 


= ‘moments and shear forces can be obtained for each 


» gi putations tobe made before forming the equations, as this tends to avoid eae de 
After the data have been entered in the tables, the filling- -in can be done me- 

- chanically by pontetants who need not be acquainted with the basic theory of 

. The tables have been used to solve Examples Iand I. In Example Ia sim- 

re zx for undisplaceable joint is used. In Example II all the tables are 


SIMPLIFICATION OF TABLES 


7 
— 
sbarscan 
w 8 
| 
“s se of systematizing the solution of the slope-deflection-gyration Jf 
| 
a 
Undisplaceabl ars that meet at the 
joint are so arranged that no movement of the joint is possible under load. It 


- then only necessary to fill in Tables ay 2- 1 and 5 for that ae 
case of Example I, solvedlater, 


a Similar simplifications are possible for the cases of fixed joints in seveten, 


flexure, or supported ones, when undisplaceable. 


a 2. Concentrated Loads at Joints Only.—In- this case it is not necessary to 


. aie Symmetrical Bar.—Here the matrices of the bar under consideration and 
the inverted bar in Tables 1 and 2 are identical. 
se 4. Symmetrical Bar With Symmetrical Loads. it symmetrical loads act 
ona symmetrical bar, the matrices of the ane caeameanes and the inverted bar 


. Symmetrical Structure with or Loads.— -In this 


case itis proven that the deformations are either equal, or equal 
Triorthogonal .—It is convenient in this case to make the 
ate axes parallel, respectively, to the torsion axes of the bars. By substituting — 
es values of the direction cosines in Table 1 and filling inthe other tables, the 7 
matrix of the system of of slope-deflection- gyration equations is obtained, whose — 
elements are explicitly expressed as functions of the rigidities and the trans- 
_ mission factors of each bar. The matrixis shownin Table5. = 
ow : Biorthogonal Frame Loaded Normally to Its Plane.—To obtain the ma- 
trix in this -case3 it is sufficient to consider the equilibrium of the moments | 
7 with respect to the X and Y axes and the equilibrium of the shear with respect | 
tothe Z axis in Table 5, since the unknowns @Z, AX, and AY are equal to zero. 
8. Biorthogonal Frame With Loads Within Its Plane. .—In this case, only the 
moment equilibrium equation with respect to axis Z and the shear equilibrium 
2 equation with respect to axes X and Z are taken from Table 5. This simplifi- 
a _ cation explicitly furnishes the matrix of the system of Maney’s slope-deflection 
equations for the structure. 
Simplifications 6, 7, and 8 refer to cases frequently met in practice. They 
can also be used to study the continuous systems of two or three dimensions. 7 
“i For that purpose these would be substituted by approximately equivalent or- 
thogonal grids, the of continuous inthe axes 


CONCLUsONE 
slope- -deflection- -gyration provides a meanto obtain the exact 
tion of space frameworks with members in any direction under any system of 
 * _ It is possibleto prepare a program for an electronic computer following 
Other space structures, as warpedsurfaces can be 


| 
- 
— 
— 
— 


4 “C4lculo de Estructuras Reticulares,” 
torial Dossat, Madrid, 1948. 


B, having fixed ends opposite to B 3). 
_ The ‘Position of the | ends of the bars are indicated in Fig. + 7 


‘ 


by Ferndndez Casado, Quinta Edi- 


Example E-1.—Undisplaceable Joint.—Considerthe case presented byFér- }§ 
se A moment is appliedat joint B which can be expressed by its cartesian com-_ a 
a ‘Table shows the values of the direction cosines, the relative rigidities 
and establish the following stmultancous equations: 


+ 
262 0, 458 oy 4, 054. 


* 


Since the joint is undisplaceable, components do not 
_ By substituting the values previously obtained inthe equations of the emeeiedl 7 
a d torsional moments and shears, the following results are obtained a ; 


641 in. Ib * 215 in. 490 in. 


341 = - 189 in. Ib 


ie 


A 


Fac = 4.211 


‘unequal 1 bennches Soom! ing a 90° horizontal engie between them, , without support 7 


_ at the common joint and fixed ends. The figure also indicates the plan, eleva- _ 


_ tion, and cross section as well as the T-axis of he branches of the stairway. 
_ The origin of coordinates coincides with joint B; the X-axis with the hori- 
zontal projection of AB; the Y-axis with the horizontal projection of BC, and 
_ the Z-axis withthe vertical projection of BC. The coordinates of the ‘Stairway 
joints are: 
BC 


000 = 0 
— 
— 
-. 
| 

Vaan =- 2.730 = 2440 Ib 
=- 11.791 l BD 12.525 lb 
The flexural momen site to joint B have 
the value of the flexural moments of these bars at joint B. The torsional 
x moments at the ends of the bars opposite joint B are equal in value to torsional 
The results obtained are equal to those found by Fernandez Casado* by the 
-moment distribution method. 
| 
ree 

i 
Uy. 


os Tables 1 to4 inclusive show ihe computation of the elements of Table 5 which 74 
to the slope-deflection-gyration equations. = = j= 
In orderto simplify the numerical computations, all the 
in terms of the flexural principal rigidity of bar AB. i. oe i 


The fo following values are obtained d by so solving the the slope- -deflection- iii: 


4 


4 = - 13,569 
7 iia the preceding amounts in Table 6 furnishes the he values tor the 


flexural and torsional moments (in ke. m): 


0 0 


| 
6,769.5 


671.0 M 
BA 


> absolute values for rotations and deflections are desired, it is 1 necessary 
to divide the relative values (obtained by solving the equations of Table 5) by 
the absolute value of flexural principal rigidity of | s AB, which is 3' hes calle 


Kgm. This way the following values are obtained: = 
radians 


00134 


0000542 radians 


000184 radians =-. 408 


= = The final moments, however, were obtained substituting directly in Table 6 
"the relative values of rotations and gyrations obtained by solving the — 


“APPENDIX I. —NOTATION 


bar AB with respect to 


29 AB’ 

cosine a with respect to Y- axe: 
‘ww 
direction cosines on i, 1-1, 2-2 2and T- T, 


bpaB 
of bar . with to Y- axis; 
= direction cosine of a bar with to Z- axis; 


Z-axis components of the axial forces of the bars that 
— 


concur at joint By Te 


direction cosine of a bar with respect to 
<4 
clively, 
on axes 1-1, 2-Zand T-T, respe — 
X-axis components of the axial forces o 
BA components of the axial forces of the bar tha 


| 


= flexural righty of of bar AB at A, with | respect to ants 1 


= flexural rigidity of bar AB at B, with respect 1 to axis 2-2; 7 
= = length of bar AB, always considered positive; — ss 


= length bar AB, with h plus si sign; 


= = length of b bar AB, ‘with minus s signs 
‘resultant of bending moments of bar AB, at at ends A and 
with h respect to axis 1-1; 
oe resultant of bending - moments of bar r AB, at ends A and 
B, respectively, with respect to axis 2-2; 
of moments of bar AB, at ends A and B, 
respectively, with respect to axis ‘T-T; 
components of of bending moments of each 
_ bar concurrent at B, in a plane perpendicular to axis XX; 
’ 


7 Marz components of the resultant of bending moments ot each 
a bar concurrent at B, in a ‘plane perpendicular to axis ZZ; 


Ma. fixed end ‘bending m moments with respect. to axis 1 of 


ria’ Mripa, bar AB, at ends A and B, respectively; 


end bending moments with respect to axis 2- of 


‘F2BA bar AB, at ends A and B, respectively; — a 
Tee , = fixed end torsion moments with respect. to axis T-T of | 


FTAB’ AB, at ends A and B, respectively; 
; = bending moments at B, in a plane per adidas to the 
pL pe 
- X-axis, due to the relative displacements of the ends of 


the bars and the axial forces of the same ‘bars, corre- 


= defined similarly to the previous case, and following ithe 
symb lisn hich cor d 
olism whi fe) is in each 


Mappy’ 


: 
a SAB = flexural rigidity of bar AB at A, with respect to axis 2-2; —_— 
BA 
“ABO 
M....M...,. 
— | 
— TAB’ TBA 
| 
BCX” 
Mpa: 
| 
framed structures these values can be disregarded; 


external moment B, ina plane perpendic= 


external moment applied at joint B, a plane 


ulartoaxis YY; 


= components along axes XX, ‘YY,a and ZZ, of 


load a lied at at B B; 
AB, on axis 1-1; 


"1B BA shear at B of bar AB, AB, on axis 1-1; 


ViaB’ af 


May 


= isostatic shears on axis 1-1, due oe the loads acting on on 
bar AB, at A and 


= isostatic ites at A and B, respectively, on axis | 2. -2, due 
the loads acting on bar AB; 
of at B of bar. axis 


. = nw 
= components along axis XX cof the resultant shear at 


— the concurrent at 


sy ; = component of resultant shear at B of bar AB, on axis YY; 
Bay’ 


Vv = components along axis YY of the resultant shear at B of 


shear at me of warAB, on axis ZZ; 


4 


V. 


VBaz’ 


| 


pn transmission fac factor with respect toaxis 1- 1, from 
A to B, on bar AB; 


& flexure transmission factor with respect toaxis 1-1, from 
 BtoA, onbar AB; 
flexure transmission factor with to 2,from 


A to B, on bar AB; eee — irom 


i 
a 


i 

i shear at B of bar AB, on axis 2-2; 
| 
| 

— — 

a 


SPACE 
il 


= displacement | of joint A, along axis 


= displacement of joint B, , along axis 1- 1; A im 
= displacement of joint A, » along a axis 2-2 


= ‘displacement of joint B, along axis 2- 2: 


deformations along the the T-axis is of the bars S concurring at 
joint B, due to axial loads Fpa and Fee: . For framed 


"structures these values are negligible; 


= rotation of joint Ay with respect to axis 1. 1; 


= rotation of joint A, with ‘Fespect to axis 
‘rotation of joint at A, with respect to axis 
rotation of B, with to axis axis 


rotatic pect to ar T- T; 
a th res ctt a 
rotath on of e pect o axis 


: = angular deformation ¢ of joint B ina pene perpendicular 


angular deformation of joint Bina plane p perpendicular 


angular deformation of joint a plane perpendicular | 
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— TABLE 2,—EXAMPLE EL 
— ako? | | € kac 655 | o | 1620; | | 868 | -244 | 262 
2 kab €kb | | 855 | O o «6669 323 | —244 | 1235 | 135 
| | 4275) © | Oo | B10 o | o | 370 188 | 2586 
2 2] 2 | | 427.5 | o | 205.1] 3339 | 155 | | 129 
&g 3 €Bkoc <€Bkbe| o | o | ° 333.9 | 10.4 | 2586 129 | -62 
&§ — 
3 0 ° | | 6869) 492.4] 5206 264 330 
Re | @fa 3 2 s7e | 724] | : 
i; 


ABLE 2. —EXAMPLE E2—CONTINUED 


| Bers _8C 


j +40720 | -28625 


| + 051862! 


6 3012 ~21.086 


| 70.356 


=f 


€ koe | | € ke? -29.669/ | 101.88 —13.939 | 46.933 Le 
— teste | | | | 
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TABLE. 2. 2,—EXAMPLE E2 E2—CONTINUED 


Bere is 


as 


ta 
a 


— 1 &4, -3.7086 | o | o | 4 
2 2 | €W,R) | €W,R) |-46.778| © | 5941 |-19889) 


1961 


| 48.031 | -.025931 
40720. 086648 


J nN), \-7 014884 | +.049737 


+ 22966 —.78264 a | 


3 


-.78264 


° - 014884 
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5. EXAMPLE E2, EQUATIONS (JOIN 


Common Jomt=— 8 


INDE PENDENT 


13148 (086648) 


5.52862 13939 


-29.669 134 148.8! 


| 


26.475 0084848 ase |-45.031 


44362 
‘| 175.77) 0 587.33 


o 
508.55) 


= v:0 5 | 22673 te 959) 4 
q 
T 1 2 T 2 T 2 
6 4 Ors | © ° ° 348.29 4436.2 535.70 
m. 
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E2—CONTINUED 


4058.! -i4162 


<| .?72 |=772 
46.716 |}-46.716|110.25 [5512 


«Bar inverted 


W9.9739-09739 SiO 25.5 
Consideres | #£Bar 


~14178 
348.29 4436.2 53570 | 


Bor inverted Inverted 
14178 -3980.! 
eae? 
4 Ore 203 4436.2 535.70 
| 


1281.0 ~39697 
2732! 


Length, in ‘Section, in in 


inches? in inches in inches 


(941,720 «341, 720000 


4 455, 625, 


— — 
Dar cue | Aor cos 

| 0, | 0, | {-s090.8| 179.9 | 1930.6 

7 
— 

Te Ga | Tse Ser | Too Se [31879 | o | o | o | | | 
© | | | | © |-28180| 0 | © o | | 
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ANALYSIS BY THE REACTION METHOD 


‘The method of reaction distribution is a convenient approach in the lution = 
of grid systems. Conventional methods of handling these systems usually ~ 
‘sult in the solution of involved equations subject to numerical errors. The : 
reaction distribution method is a convergence process that is based on the 7 

_ stiffness characteristics of the members in the system. Reactions, moments, , 
and deflections at each joint are obtainable by this method for all types of | 
; _ loading conditions. The numerical solution to various grid systems is given. 
To facilitate the presentation, a table of stiffness and carry-over relation- 

ships incorporating the scope of the spans used in the illustrative examples _ 

_ The important relationship of the converging series which greatly simpli- 
po the distribution process, is established in these problems. _ Lo i 


GENERAL ANALYSIS 


a Continuous beams are readily analyzed by the various methods | available, 7 
is assumed that the methods of consistent deformation, three moment 
theorem, slope- deflection and moment distribution, and so on, are familiar . 0 


e “accomplished in the same manner by equating the joint sania 


be 
The solution of these equations may be quite tedious, especially when several 7 


a - Note.—Discussion open until January 1, 1962. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, , ASCE. This paper is part 
of the copyrighted Journal of the Structural Division, Proceedings of the American So-_ 


ciety of Civil Engineers, Vol. 87, No. ST 6, August, 1961. 


1 Structural Engr., New York Naval Shipyard, Brooklyn, N. 4 


rg 
| 
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loading conditions enter considerations. The reaction distribution 
method offers a solution that is rather easy in application. 7 


_ The basic approach is to initially introduce a rigid support at each joint - 

- that translates. The resulting continuous ‘Spans can th thus be analyzed by using , 

any of the available methods. ‘The purpose in doing this is to obtain the re- 

= at these temporary supports. Due to joint translation, these reaction 
forces must be redistributed throughout the structure to maintain static 
equilibrium, Reaction distribution as the name is a method for accom 


SECTION == 2 (c) ) TRANSVERSE ‘SECTION 


Fig. 1 (a) is a representation of a grid consisting of one longitudinal es r- 
‘two transverse members that interact at joints B and C, ' 


There is no super- 
—— loading on the system. _A deflection (A) is induced at joint B and 7 
temporary support is placed at joint C . Fig. 1 (6) isolates the a 
: span, Fig. 1 (c) isolates the deflecting transverse span, The forces at joints 


iB and C are indicated: Rp, is the longitudinal reaction at joint B; Rpy is 
_ the transverse reaction at joint B; Roz! is the longitud 


erse linal reaction at joint C. 4 
_ Sign Convention. —The sign ‘convention used in this presentation denotes all 


‘upward forces as positive. - Moments causing tension on the bottom fiber ar fiber are 
_ considered positive, downward deflection is positive. 


_ Definitions .—Absolute translation | stiffness, T, is defined as the force re- 7 
quired to induce a deflection at t the joint in question wit with translation prevented -- 


— 


q 
A 

| 


REACTION DISTRIBUTION 

- at all the other joints. , Relative translation stiffness is defined as the force — 

required to induce a unit deflection at the joint in n question n with translation 

_ prevented at the other joints. Thus, from Fig. 1 (b) the longitudinal translation ; 

stiffness at joint B [Tai] = = Rpz (a negative force) and from Fig. 1 ©) ‘the 

transverse translation stiffness at joint [Tar] = RpT. 
_ Carry-over is defined as the forces induced at the unyielding joints s when * 


a deflection is. induced at the a in question, Therefore, in n Fig. 1 (b) RCL! is 


c carry- over force. The ratio - ——is the carry- -over factor from 
‘Table 1 nas been prepared giving stiffness and carry-over relationships 
— for beams of constant EI and spans of equal length (Fig. 2). Tew : 


a i _ Distribution factor (D.F.) is the percentage of force that is distributed in 
the or transverse direction at the in Therefore, 


= 


a in . which T is the translational eaten: 4 denotes the deflection of a joint ane 
— to translation and D. R. refers to the distributed reaction. _ For Cases 2A and 


_and the transverse distr 

: For Cased 


1961 


1; 
ET 


a2 


Cases 5A: and 


am -19.2 EIA 


92EI 


(9a) 


|) "18.85714285 EI E 


Effect of Torsional .—Ina system of interacting members, rota- 
tion at a joint causes bending on one member and torsion in the member at 


ig reasonable to aan: that the torsional resistence is negligible in com- 
=_— to the bending resistence. Therefore, torsional restraint is not con- 


_ continuous spans are analyzed by any convenient method, “thus determining 
the rigid reactions. The next step is to remove the temporary supports “fl 
<n translation. The supports are removed from one joint at a time. When 

i the support is removed a force which is equal in magnitude but opposite in i» 
; direction to the initial rigid reaction is acting on the system. This force is | 

- distributed in the longitudinal and transverse directions in accordance with 
the stiffness properties at the joint in question. The joint translation induces 7 
reactions at the other joints in the system, the magnitude and sign of which is 
determined by the carry-over factor. After distribution is accomplished at a a 

the right support thereby preventing any further 


=z 
— 
r Cases 6A an 
— 
Le 
ae Mies Reaction Distribution Process.— The first step in the process is to intro; 
| 
a us cOmpleting One cycle Of distribution, [he forces that were carried over | a 7 
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_ and is then replaced. ‘The distribution and carry- -over cycles are cepented 

until there are no more forces to distribute. The system will then be in equili- 
-brium . The final reaction at each joint is — to the summation of all the > 
“rigid reactions and distributed reactions, 

It is important to note that after each time a — reaction is distributed, 

‘the temporary support is. immediately replaced at the joint. Therefore, only 
the distributed forces cause translation. The summation of the distributed 7 

forces (or distributed reactions) is the force that induced the joint translation | 
and by definition is the absolute translation stiffness. The deflection at a 

joint be easily ‘obtained by equating the expression for stiffness (from 
. ‘Table 1) equal to the summation of the distributed reactions at the ones in 
question. 

The complete process will be easily understood with the use of several 


EXAMPLE 

ial Solution (C (Constant El). Initially place a rigid support at joint B. Because 
the 100 kip load is directly at joint B, it may be assumed to be in either -. 
longitudinal or transverse span. Therefore assume 100 kipon the longitudinal © 1 

Isolating each | an in Wi. 3 with its superimposed loading yields the 

"results indicated in Fig. 4. 
From Table 1 Case a and 14) 


04 
Longitudinal D. F. 
D.F 1.000 


= 
328084 


_ To obtain distribution factors, it is more convenient to work with the re- 
eee rather than absolute ‘stiffness factors. The stiffness expression is ob- 
tained from Table 1. For relative stiffness values, A is unity. The modulus : 
of elasticity may be omitted if it is constant. It is usually more convenient to 


multiply these stiffness by a value such that the relative 


stiffness factor is unity. 
‘The distribution factors and rigid support reactions are shown in Table + 


‘The summation of the rigid reactions at joint B is the total reaction to the 7 
unyielding ‘support (120+ 13.75 = 133. 75). When the is re- 


ST 6 

each joint are now ri 
| 

} 

| 

iil 


¥ 


d |Distribu- arry | Balanced 
Reaction| ed Reaction 
| (R.R.) | Reaction 


‘ips 


Kips 


kips 


15 kip-ft 
|  kip-ft 


23.75 kips 


FIG. 4. 


— 
Joint B Factor 

9 = a 

3.75 kips kips 8.125 kips- (13,75 kips 1.875 kips 


REAC \CTION DISTRIBUTION 


OS 


2 moved the _ will ‘translate downward, The force inducing this deflection is ’ 
= therefore (- 133.75). _ This force is distributed inthe longitudinal and transverse © 
- spans in accordance with the distribution factors. Therefore, the force (or dis- _ 
tributed reaction) in the longitudinal direction - 133. 75 (0.671916) 
-89.8688; ; the distributed reaction in the transverse. direction is - 133.75 
©. 328084) = - 43.8812. The joint is now balanced for the first cycle. Because 


all the other support points are turly rigid, no distribution takes place, con- n 


is necessary ry for this problem. The balanced ‘reaction is ‘equal to the summa- 
* a tion of the rigid and distributed reactions. The deflection at Joint B is obtained | 
wt by using | the translation stiffness relationships given for Fig. 2 (Eqs. 3 
At Joint B (longitudinal) 


= 


Joint B B (transverse) 


‘Therefore, at t Joint B “longitudinal 
— 3744, 58 ki kip f ft" 


7 Since 2 D.R. at Joint B is known, the reaction components at each joint 


obtained by combining the initial value (the value obtained with the 


temporary rigid support) with the change in the component in question due to 


This is shown in the 
‘Final. Joint Reaction Component = = R. Component (Fig. 3 3) 
(Z D.R.) (carry-over factor; from Table (13) 
Ma = = +12.5 + 0.25 25 -89. 8688) (10) ~212. 172 kip ft 
Pa 


10). = +199 
(10 + 199. 672 aps ft 


5 25 (10) - 287. 172 kip | ft 


-0. = +68. 6844 kip 
| 


=< 


| 
x 
| 
i 
| 
(1a) 
— 
— — 
— 
| 
4 


12 3.8812)= + 30.0656 kip 
= + 20.0656 kip 
(8) = +108. 525 5 kip 


‘kips kips 


ae 


1844 kips 30. 1312 hips 68.6844 kips 


« 


+160. 525 ft 


EXAMPLE TWO 
assume constant El, with the longitudinal span 10 ft 
long and the transverse spans 8 ft long (Fig. 6. 
_ For a grid with two or more deflecting joints, it is usually convenient to —_ 
> balance the system for a load at one joint at a time. In Fig . 6 a 100 kip load 
is” placed at joint B. The system will be solved for this condition due to 
symmetry. The same solution will be obtained for a 100 kip load placed at 
s C. For any superimposed loading condition the rigid reactions at these 
joints” are determined and the system is balanced by taking ratios © of ' these 
reactions to the solution for 100 kip at a joint. ee pe, 
From Table 1 Case 1 | and 2 
Relative (long.) = x - = 1.000 


Relative T T (trans.) = 


final force diagram is shown in Fig. 5. q 
— 
— 

— 1.2207 
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Longitudinal D. F. (D. F. = 0. 45031 

D. ¥.(D. D.F..p) = 0.5496 


"Assume F R.R. at Jt. B (long. ) 


a: The reaction distributions for the jolats tt that translate (B and C) are are given ae 
: in Table 3. It is not necessary to list all the other joints for there will not be © 
any carry) over from - them to joints B and C, There will be, however, carry 
over from joints B and to all the other joints. This is best accomplished 


after the joints are balanced, for the effect at any joint is determined by the — 
4 ‘summation on of the fo forces that induce the change 


C.O. from Bto ce from Cte Cto B= -0. 


100 kips 


1 

rigid — B longitudinal R. = listed, When the 
support at B is removed, the joint willtranslate downward; the force inducing 
this deflection is (- 100). This force is distributed at joint B in. accordance _ 


- with the distribution factors. . Therefore, the distributed reaction at B (trans-— 


-100 (0. 45031) 031 


Because there is no rigid reaction at joint C there is nothing to distribute for 


the first cycle, and azerocarry- -over. The is 


— 
q 
qj 
1 
q 
— 
— 


August, 1961, 
action pao at joint C equal to the distributed reaction at joint B —_ 
plied by the carry- over factor. ‘Therefore, R.R.¢c CL for the second cycle is 


The distribution process is continued | until, nothing | remains to be carried — 
Pete summation of the rigid reactions and distributed reactions is listed 


TABLE 3, DISTRIBUTION FOR 20 


| Transverse | Longitudinal Longitudinal ‘Transverse 
= 


7430 
+ 15, 525 
_-6. 


-2 .755 
086 


| 

— —— 


RR, |. | mame | 

Balanced Reaction) -65.072_ +65,072, | 


an 


“separately 3) and to give ‘the balance balanced. reaction. is neces- 

‘Sary to have 2D.R. & in order to ~obtain | the joint deflections. The 


“4 

he 522 

|g 
ap 

ore 
| 


at B transverse 3. 
- 65 65.072)(8) +5552. 8 kip ft 


3 
-21.006) (10)" 1 lip 


25. (-25.639)(8)3 = +2187.9 kip ft 


Ailsebdiiias anima at all the joints are obtained by using the carry 
factors from Table 1 and 2D. The computations are as follows: 
A=- -.375 (-53. 31) +0. (- 21. .006) = +14, 740 


Wi 
53. 31) “0.875 (- 21. 1.006) = 5.450 


ip ft 
0) = +147.40 kip 
Mc (long. )=+ 25, ¢ 58. 31) san 0. 375 c 21. py (10) = = -54.50 ki 


Mp -0.5 (-65. 072) (8) = = +260. .288 kip ft 


Rg = -0.5 (-25.639) = 2. 8195 kip= Ry ay 
Mc (trans. ) = -0.5 a ni -25. .639) (8) = +102. 556 | 556 kip f ft 


=at+ar+ar +ar 


j 

— ose me balancing cycle performed 

the first solution shows that it is a converging geometric series of the 

— a. | 

(ong. point to C][-D.F.. (long.) at point | 


= 


78 | 


1-.15525 


R.R. at int B (long.) = 
 SD.R. at B Be 
‘=DR. at JT. B ( 


> DR. R. at * 
R. at j at Cc (trans. d= = - 46.644 54969) = =e -25, 
Balanced Reaction B 65. 5.071 Balanced Cc ) = -25.640 
Balanced Reaclton B (long. ) = +65. 071 Balanced Reaction C (long. )= +25. 640 
EI at joint: B (long.) = + 5552. 8 AEIat Cc )=2 2187. 9 
A EI vat (trans.) = 2187. 9 


j A EI at joint B (trans.) = 5552.7 
series form, The converging geometric series can be used to great advantage 
= 


= 5552.7 A 
‘The procedure of distributing reactions is s greatly simplified by using this 


as will be shown in subsequent problems. 


kips pe 


—10 


6 hips 


7.643 


394 kips 


aes, 
19.737 kips- 
kip-ft? 


6.760 kips q 32.123 kips 39.474 kips | 


TABLE 4,—TABULATION OF BALANCED UNIT RIGID REACTIONS _ 


Balanced Reaction anced Reaction 
wy Components for Components for 


Ra, in kips $14,740 


Rp (long.), inkips 465,072 425.639 
Rg (trans.), in 5.072 


(trans.),inkip |  -25639 | -65.072 


Rp, inkip 5.45000 
Mp (long.), inkipft = | +147.40 | 


Mc (long.), inkipft -54.50 = | 4147.40 


Re; Rp, inkip 8195 

Mg (trans.), inkipft +260,.288 = | +4102,556 


Mc (trans.),inkipft  +102.556 =| +260,288 mall 
El,inkipft® | 455528 = | #21879 


¥ 
TABLE 5, SOLUTION OF IL LUSTR ATIVE PROBLEM (FIG, 7) 


Loading Belence Fo: For Balance For Final Reaction 
Supports = = 30.143 = 48.786 Component 
46,760 


123, 


+31.746 | 


Tap 
Ac EI | 3368.49 


— 
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— 

| 

4 

| 

| 
Rp | +2095 | |  +7.191 | 


Illustrative Example.— I ior 00 k 
has been computed. Due to the ‘symmetry of the ‘system, the same 


will hold for a R.R. of + 100 kip at joint C. The results are given in Table 4, — e 
i With this table it is possible to solve the system for any loading condition. _ 
The illustrative ‘problem given in ‘Fig. loading on the longitudinal span 
. - shown. There is no loading on the transverse spans. With temporary rigid 
supports at joints B and C the resultant continuous spans can be analyzed by | 
any of the conventional methods. _ The reaction components for the superim- 7 
posed senting condition are Ra = +4.976kip, Rp = + 30.143 kip, Rc = + 48.786 © *o 
kip, Rp = + 2.095 kip, Mp = - 30.24 kip ft, and Mc = - 39.05 kip ft. The rigid -, 
_ reactions (R. R, BL = +30. 143 kip and R.R.cL= =+ 48.786 kip) must be distributed be 
and the system solved. This is performed by applying the ratios of these rigid _ 
ses to Table 4, The solution of ” problem is given in in Table 5 5 with final 
force diagrams owe in n Fig. 8. > 


a The grid in Fig. 9 will be solved for a 100 kip load at joint B. Due to 
solutions will be obtained for 100 kip loads at the other 


“4 R. at joint * “B 3” (Long. yes +100 


‘Rel. T (trans. ) 


joint. Thus, 


"Longitudinal D.F. = - 0, 5 


factor = 


1e carry over plus the distributed reaction at each joint equals zero, a : 
fore, the balanced reaction after cycle 2 is the final joint reaction. Final 
eactions at rigid supports and joint moments can easily be obtained by simple 
tatics. Deflections at each joint can be determined by any of the conventional = 
methods. However, if solution is desired by obtaining the summation of the aa I 
distributed reactions at each joint, the distribution process must continue, 
he convergency factor is rather large so that many distribution cycles are a 
cessary to complete the balancing. Convergence may be accomplished _— 


. ty by us using the geometric series used ir in example 2. a 


. 
ry 
4 
can be solved in a manner similar to that shown in example 2, 
— 
— 
Rel. T (long.) = 


REACTION ‘DISTRIBUTION 
= 


100 Mall kip-ft 


50 kips 


— 


| 
& 
— 
— is, 
) 21.875 kips ‘ba 


Longi- Transverse Transverse | 
- 
< 
|: #19,140625] +19.140625} | O- = 
-19,140625| -19.140625| 0 -19,140625° 
| + /[+16.74805 | +16,74805 
—-=16,74805 1-16, 74805 74806 | 


Solution se Sentenm —Use is made of the limit expressed in in Eq. 16 


“After cycle 2 the dis distribution takes the 1e form of a converging series. ‘Thus be 


= D.R. (cycle 1) +D.R. (cycle 2)+LimS, ..... (18). 


r=(- 0.875)2 = (0.765625) 


At Joint B [long. and trans. ]. Ps 


494 
19.140625  _ 
-50 + 0 - 50 - 81. 6667 = -131, 6667, 


6667 13,715,2812 “iP 
ET 


Joi nt C [long. and trans.]. =e 


16 
21.875 - 16:74805 


93.3333 3333 (10° ) _ 9722.219°'P_ 


‘At Joint [long. and 
=0+0- 19. 140625 _ 
"0.234375 


— @ Joints B and F are repeated in the chart to facilitate carryover. 

(from the b 

| 
Thus 

— D. = - 93.3333 
and 


REACTION DISTRIBUTION 


Transverse 


C.0. 


Cycle 


140625 
 -19,140625 


#16,74805 +16,74805, 
-16,74805 -16,74805 


21,875 «| -50 Bal. React. 


kip 
67(10)9 -g506.948 KiP 


‘ 


at it Rigid Supports an and Joint Moments.— 


= 25(- 131. 16667) - 0. 375 (-93. 3333) = = +2, +2,0833" 
Mp = =+ 260.417 | kip ft (long. and tra trans. ) ns.) 
20. 833 kip ft - M(trans.) 


= -0.375 (-93. 3333). + 81. 6667) 5833 kip = 

4 


Ry = +0.25 .25 (-93. 3333) | - 81. =+7, kip 


M ) = +145.833 ft 
q = 
MG = = +72.917 kip (long. and t 
Final force diagrams are in ‘Figs. and 10(b). ‘Member IF 


similar to member ABCD of Fig. 10(a) for which 


= 9722. 219 kip ft? 


Longitudinal] Longitudinal | Transverse | Longitudinal 
— 
| +19,140625 | +19,140625 cycle 3 
-19,140625 | -19,140625 | D.R. 
— 
BFJ is 
— 
and — 
= 


EXAMPLE FOUR 
La Example four (Fig. 11) involves a grid with constant EI and all spans ale 
to 10 ft. Part A of this problem involves solving the grid in Fig. 11 for “7 
100 kip load at joint B and a similar analysis for a 100 kip load at joint D. 
Part B involves solving the grid for a 100 kip load at joint C. For different 


conditions of superimposed can be actives in a manner 


similar to that shown in example 2. ‘wate : 
Joints B and D, from Table 1 case 4, 


‘Thus, longitudinal (D.F.) = 0.44 and transve 


C = - 0.636363; C.O. from B to D and D to B 
aaa At Joint C, from | Table 1 Case 4 and 6B a 


Thus, longitudinal (D.F.) = 0. 38461538 | and transverse @. F.) = 0. 61538462; 


C.0.fromCtoBandCtoD=-0.8 
_ Example 4 illustrates the distribution method when carry-over to several 
‘translating joints occurs in one direction. This problem differs from the prev- 


ious ones in that carry-over to more than one joint must be considered in the 
_ balancing cycle (the same procedure applies). 


| The distributed longitudinal reaction at Joint B for a ‘BL = 


the first cycle of Table 7 (- 44). The carry-over to C Sa is - 440 
(- 0.636363) = +28.0 and to D (longitudinal) is - 44 (0.272727) = - 12. 0. Since | 


in the first cycle 1 there are no other rigid reactions, this carry-over becomes» ae. 


_ the R.R. for the second cycle. The second cycle distribution indicated in the 

_ chart yields - 10.769 at C (longitudinal) and + 5.280 at D (longitudinal). The 

; : carry over to joint B (longitudinal) is the carry over from C (longitudinal) © 
7. plus the carry-over from D (longitudinal); [ - 10.769 (- 0. 8) + 5.280 (+ 0.272727) 
8.615+1.44 = +10.055.] The distributed reaction is - 10.055 (0. 440) =-4. 424. 
The remaining x cycles are obtained it in a similar manner, each time carrying 
over to all the translating joints. wet: on 
It is noted that only the longitudinal span span is indicated in the balancing 

chart. Because there is no transverse to the d deflecting joints, 


Member KCGL is similar to member EFGH of Fig. 10(b) for which 

| 
and 
— T(trans.) = 
| ] 


— 


kips 


REACTION Di DISTRIBUTION 


_~ 


) 
61. 634 kips 2.272 272 


— 15341 kip- 45 kip-ft 


FIG. 12. 


Ss 
kip-ft 
| 
a ‘= 


TABLE 7.—REACTION DISTRIBUTION FOR EXAMPLE FOUR (PART A)® 


0 “| _ oo [e 


= "0.828 | 
+ + ,520 371 

0.402 40,262 


167) 4 oa 


40,128 20,173 


- .010 - d (026-010 


+ 014+ 


19,008” 
.006 - .002 

001 + 003 - .003] 003+ 001 001 + 

5 4,058 RR, 


-48, 24 D.R. (long) 


séReaction 


- Tran. D, F,] = 

3p. R, ) = Bal. Trans 


AET (long.), in kip ft? 


Assume rigid reaction at joint B (long.) =+ 100 kip, 
b Joint B is repeated in balancing chart to facilitate carry-over ‘process. 


Baong | | Caongy | | Dilong.) Badong) | | 
0.440 | | 0.38461 | 0,440, 
-44.0 

40 +8.615|C.0, 

| 0.910 |D.R 
59 +1,608]C.C 
826 I|D.R — 
+ 071+ 3220 | 
- 047 - + 110+ 110} - .047 - .047 - .102/C.0 
40,066 | | “40.00 | 40,066 IDR 
+ 018+ 068 |- .042- .042 + 068+ 018 + ,018 + ,068/C.0 
.010 - .026 C.0. 
| 
— 

| 

61,364 | | -18,180 | 

4 
— 


; REACTION DISTRIBUTION 


TABLE 8,—REACTION ‘FOR EXAMPLE FO FOUR (PART 


430.770} +0 +30.770 +0 
692 +0| +2 (8, 616) 
+5,302 -2(1,033) 75.502 + 
2,528 +0,795 
- .689 - 7 +2 (1. 609) 


+ .159 + (0.371) | ,159 + .990]C, 


161 20,248 +0161 | +0,161 
+ ,198 (0.102) | + + | 044 + ,198|C, 
20,106 078 106 | 0,106 ID, 
029 - .062 +2 (0,067) - 062- 029) 
011 + | -2 (0. 011 + 


006 - 
= 


0,005 | +02 0,005 

+ (0,001) + + 002 +0 


mS 
Q 


DB. 

: 

.O. 

R, 

oO. 
oO. 

R, 

0. 


Reaction 


185 | x [Tr. D.FJ= a 
_(Tr.) = Bal, Trans, Reac, 


_+3030.6 [| _+757.75 AET inkip 

See +3030.5 | +757.70 | AEI (trans.) 


Pas 


Assum e rigid reaction at joint C (long.) = + = : +100 


4 ST 6 99 
| 
-+30.770 
"$1,624 
1443 45,302 | 
| 
— | 
| 
| 
| 
4 
3 
— 
— 
| 


— 


these transverse joints need not be listed. Its omission leads to a more con-_ 
venient chart arrangement, -DR.R. and =D.R. in the longitudinal direction are 
listed separately after which the balanced longitudinal reaction is obtained. 
The balanced transverse reaction is obtained by multiplying the negative of 
> R.R. at the joint in question by its transverse distribution factor. Thisis 
more convenient than distributing after each cycle. Because there is no carry- 
over to the translating joints in transverse direction =D.R. transverse = 


balanced transverse reaction. The value of 4 El is then obtained by equating © 


the values of =>D.R. = T. 


3) 
SEIB (long. ~ * 2556. 8 kip ft" 


3 
= (-61.364) (10)3 61.364) (10) +2556. 8 kip 
11.362) (10)3_ 


24000 


+ 1.785) 


ft 
SEID kip ft 


Ped at Rigid Supports and Joint Moments (Part A). ——_ reaction 


Be Table 1) by the forces that induced the change (=D. R, adel Force diagrams 
are given in Figs. 12 (a) and 12 eee ee ee 
giv gs (a) and 12 (b). 


Ra 5681818 48. 214] +0.3 11. 0. 0681818 [+1. 785] = + 23.863 kip 


Ma = +. 2954545 5 (10) +.022727 = 130. 680 kip ft 


0454545 15 (10), = +107. 71.958 kip 


= .1590909 (10) 778 kip ft 
.272727(10) -5. 678 78 kip 
+ 2954545 (10) +5. 678 k wip 


+ 30. 682 kip = Ry 


+ 153, 410 kip ft 


+900 p Ry 
= +9.00 Rip Ry 


45 kip ft ft 


— 
‘| 
| 
if “a 
— EIB ' 
0009(10) |  -0.8(10) 
- 0.5 (-6 he { 
1 M pe 0.25 (10) = 7 = 
Mc = = 7 


tEACTION DISTRIBUTION 


2) = = 1.136 kip = Rk 


= 45, 68 k kip ft 


45.46 -ft 46 46 ft ons ) kip-ft 181.825 kip- “ft 181.82 825 kip-ft 


26: 365 5 kips 


45. 453 kip- 913 45, 453 kip- ft 0.005 kip-ft 


72.73 18. 185 Kos 4,545 kips 


‘in Table 8 and final force diagrams in Figs. 13 (a) and 13 sb). The final force 
diagram for member HDK is similar to that for FBI of Fig. 12 (a). = 


-0. 5681818 8 [- 14. 289] +0.3[-45. il 0.0681818 [- 14.289] = - 4.545 
1.0 


0681818 ABE 05681818 4.545 k 
s 4.545 kip 

(0) -.005 kip ft 
272727 (10) | .2 (10) | | -.04 0454545 (10), =. = - 45. kip ft 


+ 1590909 (10) +90. 913 kip tt 
272727 (10) = - 45.453 kip ft 
=+ 022727 7 (10) (10) } =-.005kip 
Transverse Direction.— 
5 (- 185) = = +9.0925 Kip= Ry 
= 45, 46 kip ft = My 


— 


ngcycleisgiven 

9.0925 kips 9.0925 kips_ 

0.005 kip-ft 

4.545 kips 18.185 kips 

4 

— 

Ma = +, 
— 
= | M =-, 

4 
> 

| a 


= 0.5 ( = + 36. 365 


5 (10) = = -181. ft = 


CONCLUSIONS 


& The examples contained herein demonst rate the u use eof the method. The greatest , 
advantage of this method is that once understood it becomes a mechanical © 
_ operation. The solution of grid systems is not limited to the types of spans 
given in Table 1. For different types of systems expressions for stiffness @ 
~ and carry-over would have to be developed. Grids with more interacting joints 7 ae 
than those used in the illustrative problems can be solved. It can readily be 
seen, however, that the « carry-over process becomes more ‘complex with a 
_increased number of translating joints. These systems can be subdivided into 
grids similar to those solved in this presentation and then combined for final 


balancing. It is felt that the — of the more complex grid systems an 


— 


The notation u used in this paper are presented here for convenience of refer- 


and for the use of ¢ discussers: 
ar carryover; 


distribution on factor; 


+ distribution 
transverse distribution factor; 

= distributed reaction; 


= modulus a of elasticity; 

= ‘moment of of inertia; 


reactions; 


= transverse 


= rigid reaction; 


= longitudinal rigid 


ST6 
@ 
: 
Ag 
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_ REACTION DIST RIBUTION 


= longitudinal translation ‘stiffnesses; 
= transverse translation stiffnesses; and 


= deflection c of a joint due to translation. 


tose. 


. 


i> 


<5 


—_— | 
| 


ses 
=f) 


— 
— 
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STRUCTURAL DIVISION: 
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A structural model that closely resembles the plate is anal- 
-yzed for the purpose of solving problems of bending of rectangular plates nu-— 
_ merically. The plate model consists of two groups of beams intersecting at i 
‘right angles. Beams of one group are considered to be continuous over the 7 
- beams of the other group, free to rotate at points of support and capable of re- | 
_ ‘sisting bending moments only. The twisting moments in the plate are carried 
. pairs of coil springs, connected parallel to the beams at points halfway be- 
tween beam | intersections. Loads are assumed to be concentrated at beam in- 
: tersections, In In comparison with solutions by first order fir finite differences or 
by other plate | models, the procedure -e described herein has two main n advan- 
_ tages: (1) the elastic surface, moments, etc., are obtained with a higher de- 
gree of and id (2) only a small number of simultaneous « equations 


INT RODUCTION 


snatheie of analysis in the field of structural engineering became a neces- 
sity in order to avoid the complex mathematics involvedinthe solution of many 
_ problems by the formal analytical methods, Much research has been success- a 
fully: conducted in ‘this line by engineers, definite numerical 


ad 


a "Note. —Discussion open until January 1, 1962. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
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; have been established for the solution of problems related to beams, columns, 
plates, andother elastic bodies, 
ala The problem of bending of rectangular plates 1 received great attention, and 
solutions are now obtainable by the finite difference method,2 and by tu use e of 
structural models.3,4 In both approaches, it is required in general to solve a 
-number of simultaneous linear equations. The accuracy of the results increases 
with the fineness of the division of the plate surface. Nevertheless, a rapidly 
increasing number of simultaneous options wa then have to be solved, and 
the solution becomes too complicated. The object of this investigation is to 
overcome this difficulty, that is, to reducethe amount of computations as" much © 


as possible, and at the same time obtain results with increased accuracy. ‘This 


plate. model and ‘its are described herein, and solutions are ob- 
tained for plates 1 with simply supported, built-in, free, a and elastically supported 
edges, as well as continuous plates. Plates with point supports can be also 
handled by this method without difficulty. The affect of Poisson’ s ratio onmo-_ 
ments, , deflections, etc., is considered at interior points, . and at points on sup- = 

portedandfree edges, 

_ ‘The analysis of the plate model is simple and yields results with a high de- 

q gree of accuracy. In each case, the solution is obtained through analysis of 
models, and the results are extrapolated. This procedure yields answers 
which are approximately 1% in error. 

_ Generalized Scope.—The treatment by the presented structural model is 

general, _ It covers a broad scope of applicability, and it is not limited to bend-— 

| 7 ing of plates. By analogies, its application can be extended to boundary value 

problems in elasticity which involve the solution of Laplace’s or Poisson’s 
equations, as well as plane stress —” strain problems which are areusually — 


to 


PLATE 
_ Deflections and Moments of the Plate.—The elastic deformation yn of a later- 


aly is ‘governed | by the 


ax! By” 


in wee deflection of the middle of the P denotes ‘the 
Et8/[12 - u?)] re represents the eflexural 


. 


the thickness of the plate, and denotes Poisson’s ratio. __ 


2 “Numerical Computation of Buckling Loads by Finite Differences, ” by 
"Salvadori, Transactions, ASCE, Vol. 116, 1951, pp. 590- 624. 
_ 3 “Deflections in Gridworks and Slabs,” by Walter W. Ewell, ‘Shigeo Okubo, and Joel 
_L Abrams, Transactions, ASCE, Vol. 117, 1952, pp. 869-890. = = 
_* “Numerical Methods of Analysis in Engineering,” Symposium presented at Illinois 7 . 
Inst. of Technology, at Chicago, in 1949; The MacMillan Co., New York, 1949, p . 142: 
“Numerical Methods of Analysis of Bars, Plates, and Elastic Bodies, ” by N. M. New. 7 
5 “Theory of Plates and Shells, ” by S. Timoshenko, McGraw- Hill Book Co., Inc. j New 


~~ 
— 
— 
&g 
— 
— 
an 


relation§ betweenthe bending moments My and My, the twisting moment 


Description of the Model.—To ‘determinet the he structural m model, ‘consider the 

rectangular plate ABCD shown in Fig. .1(a), and imagine it to be divided into 
a number of strips by lines parallel to the sides of the plate. Lines like 1-1 : 
and 2-2 divide the plate into strips extending in the x- direction, and lines like 
3'-3" divide it into strips parallel to the y-direction, The inner 
‘strips: are all equal to A in width, while the edge strips are only half as wide. i. 

These strips are represented in Fig. 1(b) as beams. Coil springs of length A 


«* are connected at at ‘er ends to these / beams at points ha half-way between t _ 


Division 


lines Coil springs Inner beams 


intersections. The structure in Fig. 1(b) is the modelof the sete in Fig. ‘A(a). 
Beams extending in one direction are considered to be continuous and freely 
~ supported by beams of the other direction, The load on a square area like 
5 -6-7-8 in common between the two intersecting beams ab ana cd is considered | 
as concentrated at their intersection point e. Similarly the load on rectangle 
2 -3'-5-6 is replaced by a concentrated load at point c on the edge. 
_ Elastic Properties and Deformations of the Model.—In order to’ describe 
the: elastic behavior of the aes “mn distinction is made between its two com- 
namel 
(A) the continuous beams venient in the x- and y- -directions, which are 
considered to have only flexur: 


> 
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which provide. the tortional resistance of the plate 
= 

The load p, at an intersection point, is the sum of three portions p,, Py, and 
Pxy resisted, respectively, by bending of the x- and y- beams, and by torsion 
of coil springs. The “differential- -difference” equation of ‘the elastic surface 


Pxy 


in which N is the per unit width of the | Poisson’ 


“The two brackets 1 in ‘Eqs. 3 4 indicate that the inclosed terms are in differ 


wx 

2 Moy “,,)* - 2 - - 2M’ 


"moment in the coil springs of panel wn west- north of point 


= 
fot] 
in which Mx and My are respectively the 
a ve: A O LO na nal weet enuth | ij 
le Th A- Al nt in thea cnil sorinos of the 
Mo ic the twicting mome 


q 


in which the numbers: are coefficients of the deflections a the intersection 


points. Eq. 7 also gives values of Pxy at o and the surrounding points for aunit: 
deflection at o and zero elsewhere, In this case the twisting moments 


Between Plate and Model. _The us model isa representation of 
é he plate by which it is possible to Gctermine approximate deflections of the =U 

_ Comparisonof Eqs. 1 and3 shows the analogy between the differential equa-_ 
tions for the elastic surface of the plate and its model, with exception of the — 
second term of Eq. 3 which is in differenceform, Nevertheless, with increas- 

S = number of intersection points, this term approaches the corresponding one 

of Eq. 1. The elastic surface of the model istherefore an approximation of the 
_ elastic surface of the plate itself, in case Poisson’s ratio is zero. . To deter- *% 


mine the | plant deflections when Poisson’s ratio is different from zero, it is 
only necessary to substitute the proper value of N in the obtained deflections | 


a: _ Comparison of Eqs. 2 and 4 shows that the moments of the plate model rep- 
resent the moments in the plate when Poisson’s ratio is zero. The plate mo- 

_ ments, shears, and reactive forces for any value of Poisson’ s ratio different 


- from zero are obtained from model moments by use of the following relations _ 


The plate moments 


| 
7 


1961. 


- ‘The shearing forces at sections of the plate located halfway | between the beam 
"intersection points are 


‘These s shearing forces are independent of and hence are ‘equal to the corve- 


sponding shearing forces of the plate model. The reactive forces at an inter- 
_ section point on a supported edge of the plate aa are 7 >. 


aM 


in which P isthe menanat loadat the intersection point of the lini which» 


must not be overlookedin the computation of edge reactions, The pt sign | 


tive _ Sign, for on the lower and right edges. The plate corner 


OR. 2M, 2 - owe M. 


a v 
= The derivatives and differences of Eqs. 12 can be easily evaluated as will be 
shown later in the solved examples, 
At a boundary, the beams of the plate model are made built-in, simply sup- 
ported, continuous, _ elastically supported, or on point support where the plate 
is likewise su supported. vy At free and elastically supported edges, the boundary © 


conditions are more complicated, and special consideration is 


of plates with free and elastically edges are presented subse- 
ANALYSIS OF THE MODEL 
Model Unit. .—Consider a square plate, length c of si k built- in 
= four four edges. . Fig. 2 shows the plate model model with the coil springs r removed ~ 
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convenience. Let the intersection point undergo a downward unit deflec- 


are registered. The bending moments are “written above the x-beams, and on on 
the west side of the y-beams. Thetwisting moments are  scoredat the centers — 
of the panels. _ The numbers are all coefficients of N/A2, The necessary hold- 
ing force at o is found by Eq. 5to be 56 N/A2, 4 
_ Method of Analysis .—Consider the plate model shown in Fig. 3, whichis sup-_ 

portedon its four edges. Let i i represent points 1 to 9, and n 1 be one of of these 


— 


mes 
. 
points, for example, 5. The analysis of the model n may be summarized 


(1) Allow a unit deflection (wp = = 1) at only one intersection point n) 
& beams being locked against both rotation and deflection at all other inter inter- - 
(2) Determine the moments in the basic model unit (1-3-7-9) which has 
the deflected point at the center. The moments are as givenin Fig. 2, 
(3) Allow the beams to rotate. Distribute the unbalanced fixed end moments 
moment distribution procedure, The twisting moments are not altered dur- 


(4) _ Determine the self- of holding forces wget)” using 


Eq. 5. Holding forces are re only needed at intersection points on the X- dy- 


beams through the deflected point, and at the four ey leet points 


(5) — steps 1 to 4 for all intersection im, 


_PLATES _ 111, 
tudes +6 N/d2 at o, and -6 N/d2 at ends w, n, e, ands. Nomoments are ine 
. _ duced in the edge beams as a result of this displacement. Thetwisting mo- [% 
=| 
; 
4 
— 


ra 
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(6) 


‘Write an equation ~— every intersection point, stating that | the s sum ot 
holding fo forces s (R, )at the the point, with true 


"must be equal to the applied load (Pq), thus 
ware 
The holding force at point (n) due toa —_— deflection at point bath is equal 


=R _ 
i,Wpy=1 


on” (Eq. 13(a)), becomes _ 


(7) The moments, shears, and reactions of the plate are then ee 
the elastic surface of the model. The model moments are determined by su- 


Perposition of the moments in step. 3, ‘thus at point | (n) 


i‘ — result isa number of ‘simultaneous equations equal to the number of un- - 
known deflections. In view of the reciprocal relationship (Eq. 13(b) ), the equa- 
tions are of the normal type and can be conveniently solved by ee | 


and in any panel wn, the twisting moment | 


distribution ee Consider the continuous beam AB, with equal spans A, 
shown in Fig. 4(a). An unbalanced moment at support n can be distributed in 7 
the ratio of the stiffnesses Kg and Kp of the two continuous beams nA and nB. 
By the stiffness K,,, for example, is meant the moment at n required to produce 
a unit rotation of the beam nB at end n (Fig. 4(b)). The distributed moments 
~ 4 are ‘carried over to all all the supports by use of carry- over factors inthe differ- 
ent: spans. carry- -over factor from support n to ‘support 2" is Mo'/M,; ; 
from 2" to 1', My'/Mg'; and so on (Fig. 4(c)). Inthis way an unbalanced mo- 


ment is distributed and ae over to the — of the beam a only one 7 
4 
‘distribution, 


i 
cedure.—Many computations for the distribution 
| @ if 


> 


4 
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Stiftnesses:, 


468 627 3385 
(a) FAR END A HINGED 


| 
7 
a 
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(b) FAR END A BUILT-IN aM. 


Fig. 5s shows the and carry-over factors ready for use ise inthe 
distribution _ procedure, in cases of beams with far ends hinged and built- in, 
_ For illustration of the suggested simplification, consider the beam AB — 
in Fig. 6(a). The stiffness of beam 2-A is found in line (6) of Fig. 5 under ~ 
‘support 2, to be 7/2, and the stiffness of beam 2' -B is given in line (3) under 
a 2 as 24/7 . Similarly the stiffness” ‘of beams 1-A and 3' -B, 3-A and 
-1'-B can be obtained. The rest of Fig. 6(a) is self explaining. The bending 
ome diagram is drawn in Fig. 6(b). The holding forces and the elastic line 


of the beam are shown in Fig. — 


Distribution at 
carry- -over: 


©) Distribution atl 
and carry-over: 


Moments: 


q 
be the holding forces R, and with bending of beams EF 
and GH, produced bya a ‘unit deflection at forces Rxy ac- 


— 

: 

Common 

FIG. 

i 
50 
ane added in line (4). In Fig. | — 
[ and ana are 
Mi, and M! y, registered in the same order 7 
ribed previously under the heading “Basic de 
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a of intensity p, the load- deflection equation (Eq. 
14) can be at 5, thus” 


- 1534w, + +194 Ws - 1534 w, + 
Noticing that in this case there is symmetry with oneans to the diagonal line - 
_AD, the elastic surface is also — 


97 w. - 1534 w, 
= 


131 + 97 16c 
Five other equations at | points 1, 2, 3, 6, and 9 canbe written and solved to de-— 


_ termine the elastic surface of the In cases of ‘it is feasible 


) HOLDING F FORCES 2 


Holding Forces and Moments .—It was seen the he holding 
and moments of the model (Fig. 7) were obtained with great ease ina rela- — 
tively short time, after the holding forces and moments of the beams | (Figs. 

6(b) and 6(c)) have been determined. In orderto secure the same 

in other cases, Table 1 has been prepared, It contains the results of analyses _ 
wt of several beams by the | procedure previously described, The beams have con- _ 

tinuous spans in number vom two to six with combinations of simply 


— 
j 
1 
— 
| 
— 
= W,» Wy = W,, and we = We { 
— 
(3) Ry | |= 338 776 |- 388 
(2) Ry |= 930 |. 0 gh 
rie) } to allow unit deflections at symmetrical points simultaneously, thus making it 7 
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AT INTERSECTION POINTS 
Type? of Unit} ding Forces And Bendine Moments 
Support at Def- >, 
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TINUED 


TABLE 1, -CONTINUE 


Holding Forces And Bending 


re b = 
y At Support 


51 -36 9 


-51 138 -132 


~ 


198-182 


426 - 270 a as 0 


696 -1038 432 


-540 498 728872 

1824-1146 432-72 
“144-288 


712 -252 354 
432 7930 960 


a2 42 - 
48 <-72 252 -354 156 


918 -582 
936 36 - -252 


168 -48 
3936 -2490 1008 


s. 


4 


me 
» 


| 


624 930-588 


a 0 762 -540 144 -36 N/209A2 


-162 2064 -1986 864-216 36. 


7 


9-540 906 -576 


4 


D 


795 -504 4 135-36 9 
1299 1938 810 ~216 


-1008 930 -540 144 a 


270 -540 804 -504 12% OF 
-2154 2652 -1938 756 7126 


-1938 2544 -1722 468 


54-216 -1722 1788 


36 -126 468 -660 291 


571 -362 97 -2% 7 -2 
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— 
Support'at_| Def com 
Case | > | Le ey Factor 
gna | at Fa 
3) | | ©] | N/a 
els | b 
A 
N/181 
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“Holding Forces s And 


668 -388 104 -28 


-1392 2448 -1548 624 7168 

194-388 - 98 

“1548 1 1908 -1404 “168 


4 


156 624 -1404_ -1404 


-209 
7474336 


4m 
ATA 1284 (71236 540 “144 


9-336 5 564 
336 1824 - 1380 576-144 


-90 540 -1380 540 


4 1 


14466 6048 =1620 32 -108 


6942 -4032 1080 72 
432 


2016 -4032 6006 -3780 1008 -252 
6048 16086 19824 - “14574 6048 “1512, 


1080 -3780 5934 -3744 936 
6480-14574 19392 14358 5616 


| 144 -288 1008 -3744 5862 -3492 oj. 
< 432 -1728 6048 14358 18960-12846 3492 


72 -252 936 -3492 4926 
(7108 432 -1512 561 5616 
12 42 -156 582 -2172 
-72 252 -936 “4926 2172 


Ls 
~ 


a. 


ao io 


‘denotes simple s support; b, built- -in support. a 
b Bending moments are written in upper lines; holding f forces, in lower lines, 


anes and built- in ends. Table 1 is useful in the e analysis of plate models" 
_ when L/A is not more than 6, or 12 in cases of symmetry. - The table has been 
_ prepared for inner beams of the model which have widths equal to r, so that in 
case of edge beams ‘with widths /2, the holding forces are only half the values. 
Increased Accuracy by Extrapolation.—In general, accuracy of the results 
_ depends on the fineness of the chosen nets. _ Nevertheless it is feasible to use ~ 


srs 
— 

ia: "Support at | Def — 
6 |b] b 2 

— 6 | s 

— 
— 


PLATES 


less finely divided nets, ‘s, and. extrapolate the results, although it becomes” 
necessary to repeat the solution at least withtwo different nets, The extrapo- 
lation is made herein by the procedure suggested by L. BP. Richardson. 7, 8 In 
general the approximations kj, kj corresponding to meshes with nj and nj divi- 
_ sions are assumed to approach the e true value k monotonically. The extrapo- 
latedvalueisthen 


- PLATES SUPPORTED ALONG EDGES 
‘Three e plates : supported along all the edges are now oneuinie analyzed to 
“illustrate the details and arrangement of the solution by the method described — 
herein, _ The results are compared with those available by the formal analyti- 
cal methods, in order to examine the accuracy of the procedure. ee 
_ Uniforml y Loaded Square 2 Plate Simpl y Supported Along | Four Edges Fig. 
8shows two models of the plate. In model athe lengthof side L= - 2 A; in model © 
b L = 4A. In model a (Fig. 8(a)) the deflection at the center point 1 has to 
be determined. Figs. 9(a) and 9(b) show, respectively, the holding forces and — 
moments ‘consistent with a unit deflection at point 1 as obtained from Table 1 
"case 1d. The holding force at point 1 must be equal to the applied load, hence _ 


~~ which the center deflection of the athe is 


(18) 


sa “The Arithmatical Solution by Finite Differences of Physical Prob- 
lems Involving Differential Equations, with an Application to the Stresses ina Masonry > 
Dam,” by L. F. Richardson, I Transactions, Royal Soc. of London, Series — 
8 “Numerical Computation of Buckling by Finite Differences, ” by Mario | G. 
Salvadori, Transactions, ASCE, Vol. 116, 1951, Eq. 41, p. 608. 


a 

j 

7 

andthe center moments are 

tin panel 1A, for example, 

 &§ M,, = = 0.0125 p L? (18d) 


= 0.3 30, the plate deflection and moments a 
0.003125 x 12 - 0.09) P = 0, 03412 P 


i. 

= 0.0375. (1 + +0. p pL? =¢ 00. 


= - Mp (1 


[ (0.08: 0375 - 0) -0, 0.0185 

-0 0125) | pL? x L* x L = 0.1 5 


oo at G, the effect of concentrated load 0. us p nthe at G 
= 0.375 Pp he 


and the reaction Vy at Gis determined from. Eq. 12(d), hence 
= | (0.0375 - 0) (2-0,30)(-0,0125 
0125) + 0, 125 | pL? x 0.41 p 


From Eq. (12(e), the concentrated reaction at corner is 
Re x 0,00875 p L2 = 0.0175 p 


In ‘model b (Fig. 8(b)), the deflections at only three pc points musi must be deter- 


mined, in view of symmetry. The holding forces and moments consistent with 
—_ deflections at the intersection points are obtained from Table 1, cases 3e 
and 3f. Case 3e however ‘requires | a modification in order to account for, sym- 

: metry. . Unit deflections at the two quarter points are allowed simultaneously. 
; The holding forces and moments can be readily obtained by superposition. De- - 
: tails of analysis of models b are omitted for briefness, but the results of anal-— 
-yses by plate model a and b are summarized in Table 2, extrapolated and com- 
pared with the corresponding values given by S. Timoshenko. In models a and 

ng=2 and and and the extrapolated value in in Eq. 17 becomes 


— 
gg 
oA: 
— 
The modeT anc 
section of beam I-G1 
‘and to find the sheaz 
is added, hence 
— 


results given in Table 2 (2) (3) are in column» 
. The exact obtained by Timoshenko are given? in column n (5). Com-— 


q 


parison of columns (4) 2 and (5) shows the close agreement between the values 


obtained by analysis of the suggested plate model and the exact values . For 
example, the center deflection and center moment are, respectively, 1.6% and 


York, , 1940, p.133, Table5. 


_ 9% “Theory of Plates and Shells,” by S. Timoshenko, McGraw-Hill Book ok Co., » Inc. a 


— 


= 4% of the exe exact values. The small of equations which need to be 


Pe Uniformly Loaded Square Plate Built-in Along wee Edges —In general, 


more divisions are needed in the analysis of plates with built- in —_ than i 


2.—UNIFORMLY LOADED PLATE ‘SIMPLY y SUPPORTED 


‘Results Results by ‘tin Results by 
Model a Model Results Timoshenko 


Center Deflection, 


Center Moment, | 0.04875 


at ‘Mid-— 
edge, Vm 


40.0125 


0125 


(a) MODEL ‘MOMENTS 


case > of plates with simply supported edges, if results with the same degree of 


accuracy are desired. For this reason the two plate models a and b shown in | 


‘Fig. 11, with sides L 7 4 and L = 6), respectively, are considered, fe 


i 
— 
Common: 
— (b) PLATE MOMENTS (w=0.30) 
— 
— 


PLATES 


values of — 
TABLE 3.- 3,—UNIFORMLY equant PLATE BUILT- 


_ ALONG FOUR EDGES, (# = 0,30) 


Results by Results by Extrapola- Results | ‘Come 


(6) 


Results 


Center Deflection, 9 01393 


Mid-edge ‘Moment, 


negative 


Uniformly Loaded ‘Continuous Plate, Simply Supported Along All Edges “yur 
‘ton Intermediate Supports. —The procedure is also applicable to continuous 
‘pilates. . For examining, the degree of accuracy of the results the plate shown — 
in Fig. 12(a) is analyzed, and only the results of the analysis are presented _ 
‘herein, The plate is simply along the edges andalso along interme- 


7 


In model a it is necessary to determine the dellections at points 1, 
ws C 33 in model b, at points 1 to 6. Details of the analysis of models a and b are 1 7 
>: ne omitted for briefness, but the results are given in Table 3, in columns (2) and | 
7 These results are extranolatedin column (4) and compared with the 

pL4/(Et>) 

0.02275, pL2 4 

--0.0513 =| 

4 

; Comparison of Table 3 columns (4) and (5) shows the outstanding agreement | : 

— 

= 


TABLE 4, UNIFORMLY LOADED CONTINUOUS ‘PLATE, SIMPLY 
SUPPORTED ALONG ALL EDGES AND ALONG 
'TERMEDIATE SUPPORTS, (“= 0. 


Results by| Results by xtrapo a- | Co: 

osition tion of 


_ Center of Middle | Deflection, w) | 0.001525 | 0.0018 | 0.001433 
Moment, Mj, 0.0302 0,0248 0.023 


Middle of Inter- | Moment, Mgx | -0.0426 | -0,0672 |-0,0754 
(point3) T Moment, Mg 70,01278 | -0,02 -0,02262 
enter of Si eflection, w5 | 0. 85 

f Sid  Deflecti 5 | 0.002325 | 0.00282} 0.0029 


Span 5) | Moment, Ms, 0.0417 0.0399 


Moment, M5y 0.0361 | 0.035 


= (a) PLATE 


‘MODEL a 


— ae > 


| 


© ) MODEL 


— — 
actor 
| 
tm 
ee 


= 


modified, by adding | the 
My 
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Two models a and b in ‘which L = 2, and L= 4 shown in Figs. 12(b) 


and 12(c) are analyzed. In model a the deflections at point 1 and 5 arethe only | 
_ two unknowns. In model b, the deflections at ten points must be determined. | 
The deflection at points” 3 and 9 on support ss | is zero. _ The load-deflection © 
_ equations in both cases can be readily determined by use of Table 1 as — e 
| The moments and deflections at the center of the middle span (point 1), my seed 
sdiate support (point 3), and at the center of the side 
span (point 5) as obtained by analysis of models a and b are summarized and ~ 
extrapolated in Table 4. The moment at the middle of the intermediate support 
is -0. 0754 p L2, magnitude of the same moment as givenby Timoshenkol0 
is -0,0763 p 12, The two values are in good 


PLATES WITH FREE EDGE ES 
Boundary Conditions.—It was shown in the general solution of 
-portedon all edges that the effect of Poisson’s ratio canbe disregarded in de- : 
— the elastic surface of the plate model, This is, however, not possible 
in case of plates with free edges, because the boundary conditions are some- 
what more involved, For example at a free in x- ‘-direction, 


the boundary conditions are 
: 
and 
in which M_ and Vy are as given by Ses. 10 and the ge general ‘procedure. 


is the moment M' is z zero a and the plate moment M _ will not vanish at 

the free edge, resulting in residual But if the procedure is 


= 


4 


at the free edge, then boundary condition anne is satisfied. * this case ~~ 


at 


‘Similarly, if the effect of Poisson’ Ss ratio. is until. the 


surface is determined, residual reactive forces of magnitude y 


arise at the free edge of the plate. Hence reactive forces of setae 


10 “Theory of Plates and Ss. Book Co., Inc., 


| 
— 
— 
an 
— 
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xy 
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“Fig. 5. "These moments, in turn, require additional R' at the 


intersection forces V,, also ‘require | additional holding | 


‘Snir conclusions: are by consideration ota a free to the 


y 


are — at the free edge, and the additional holding forces Ry and Ry are 
obtained as previously described, 


a The holding forces Rx’ and Ry can be conveniently determined by use of © 

the “multi- operator” ‘shown in Fig. 13 in which the coefficients of the twisting a 
moments Mxy are given, The operator in thisform is applicable for all points 
on free edges extending in the direction n-s or w-e, , as well as corner points - | 
ontwo free edges. The operator is placedonthe model with point o on the edge © 


and of the operator lying outside the model are disregarded 


> L free, the other edges eupported. And let the ratio a= 2, 
‘Fig. 14 shows tw two models a and b of the plate, in which L = 2 A and L = 
‘respectively. Fig. 15 shows the moments and the holding forces in pad “a 
produced by a unit deflection at intersection point 1. In Fig. 15(a) the moments — 
in the edge beam AB are given in line (1) as obtained from Table 1, case 1d. 
moments in| y-beam 1-E are written in column (1). In column 
the additional couple M) = = -0.30 —> = -0,90 at point 1 and zero at point 


7 E is given. Thetotal moments are given in column (3). In Fig. 15(b) the hold- 
ing force at point 1 is computed as usual in lines (1), (2), aa (3). In line 
(4), ‘the | additional holding fox force 3@ Ry is given. hh In line (5), the additional hold- i, 


must be added at the free edge in order t residu 
ome _In the analysis of the model, the couples produce addit: 
q 
— 
| 
hs 
4 


‘ 


MODEL a 


| 


| 
Fo. 


ing force Ry by use use of the ope: operator a Fig. 13 is is 5 presented. In line (6) the © 7 
sum of these forces is determined. re - the load- deflection equation (Eq. . 14) 


and the bending moment 


i In nate b the deflection at four points only must be determined in view of 
ane. ' The analysis of the model is omitted. In Table 5 the maximum de- 


‘TABLE 5, ,—UNIFORMLY LOADED RECTANGULAR PLATE SIMPLY SUPPORTED 
“ALONG THREE EDGES, THE FREE (L/L) = 


Results by Results by Results by Common 
Model tion of |Timoshenko | Factor 


Maximum Deflec- 
tion, an 0. 0737 0.0775 


Maximum a 0. 0603 ( 0597 | 0.060 


flection and bending moment (at point 1) as obtained d by analysis of models a 


and b are given, extrapolated and compared with the results given by Timo- 
shenko,11 _ Comparison of the values given in columns (4) and (5) of the table | 


of the plate models and by Timoshenko, ~ 


of ‘its edges 2 are supported by elastic beams, and assume that these beams re-_ 


an 11 “Theory of Plates and Shells,” 
New York, 1940, 219, 


7 
— 
— — — 
7 from which the deflection of the — 
M _ = 3 w, = 0.06818 p L2 

my 

rk 
pL? 
shows the Outstanding agreement between the resultS aS ODtained Dy analysis 
by S. Timoshenko, McGraw-Hill Book Co., 
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PLATES 


The boundary conditions at an edge to the 


in which H is the flexural rigidity of the supporting beam. 


“MODEL 


‘The ‘conditions given by 
free edge of the plate Eqs. 21. “As explained previously, here again additional a 
couples 3 M’ (Eq. 22) and additional | reactive forces (Eq. 2 (24) must be be applied = 


the e ‘edge of the model, in to avoid residuals at the of 
"plate. To find the moments and holding forces consistent with a unit deflection | ? 
at any intersection point of the model, first the moments and holding forces 7 y, 
_ are obtainedas if the edge were free, To the holding forces thus obtained, the 
‘forces Rpx necessary to hold the supporting beam in its deflected form must — 
be added. These latter forces canbe obtained from Table 1 if N inthe common > 
factor is replaced by Np = H/A. And to find the supporting beam moments, = 


Produced! by unit deflections, from the table, N is replaced nas H int the common 

Uniform! y Loaded Square Plate Simply Supported Along Two Opposite Edges 
and the Other wo Elastically Supported. ‘square 


which Py, isthe reactive force transmitted from 
i beam which undergoes the same deflection as the a 
* 


e en flexural rigidity H suchthat H/(NL) = 2, ane z= 0, 30. The plate is uni- | 
_ formly loaded. Figs. 16(b) and 16(c) show two plate r models a and b, in asad 


L=2 Aand L = 4 A, respectively, 


L re 


« 


(2)Ry 


AT POINT 5 


3 ‘In model a the deflections at the two. points 1 and 5 must be determined. 
Fig. , shows the moments produced in the model by unit deflections at these 
"points, as obtained from Table 1, cases 1c and 1d, The moments are registered | 


inthe same order described previously. Fig. 18 shows the holding forces cor-_ 


— 

| 
— FIG, 18 


PLATES 


responding to (anne unit deflections. The forces given in lines (1) to (5) 1 have | 
_ the same significance as explained previously. In line (6) the forces Rpx which 
_ are obtained from Table 1, case 1d, and Np = = 2 H/L = 4 N are written. 
The holding forces in lines (1) to (6) are added in line (7), a 
Th loa load-deflection equations: (Bags. 13 13 (a)) at 1 5 


pad 

wee 


- The solution of these equations give the plate deflections 


94 PX’ = 0.0041 41434 134 P Pe 


from whic 


= 0. 03133 P 


___ The plate moments are oman trom the model ‘moments. by us use of Eqs. 
= 0. 05912 


od inn’, of model b is made inthe same way, and deflections an points 1 to 6. 
are determined. The analysis however is omitted for briefness, but the results — 
are given n herein. . Table 6 contains the results as obtained by analysis of mod- 
els a and b and extrapolation of these results, as well as the results computed 
by the writer from the general solution given by Timoshenko.12 The agree- _ 


12 “Theory of Plates and Shells,” by S. Timoshenko, McGraw-Hill Book Co., Inc., “a 


New York, , 1940, Article 43, Eqs. | (g), on p. 221; ; Article 29, Eqs. of and 


= 
&g 
4 
| 
= 0.0662 PE (stay 
w, = 0.028839 ~—— = 0.0018024 = 0.01968 oo 
the plate model moments are obtained, using Eqs. 15, hence at point 
=(3 w, - 2.55 w 
Bay 
— 


TABLE | 6.—UNIFORMLY Y LOADED SQUARE PLATE SIMPLY Y SUPPORTED ALONG 
_ TWO OPPOSITE EDGES AND ELASTICALLYS U PPORTED ALONG 


Results by Extrapola- Results by Common 
del b tion of Timoshenko Factor 


Results 
| 


9.05523 | 0.05856 | 0.05893 | pl 
¢ 
0.04514 


‘ment between the results in columns (4) and (5) of Table 6 shows the mane de- 
gree of accuracy of the results as obtained by the presented procedure, | ns 


é complete analysis by the first order finite difference procedure was made 
in the case of three of the previously ee plates, namely: a 


(a) ‘Uniformly | loaded square re plate, simply ‘supported along four edges; 
(b) _ Uniformly loaded square plate, built-in along four edges; and ea 
_(c) Uniformly loaded rectangular plate, simply supported along three edg 


the fourth edge free. 


> results as obtained by both procedures. It is noted from column (9), of the table 

_ that, with exception of the center deflection of the simply supported square 
_ plate, the deflections and moments have a higher degree of accuracy when ob- 
7 _ tained by analysis of the plate model. ‘ In this exceptional case, however, the 

error is only 1.58%. In case of the square plate built-in along the four edges, 


the number of equations involved in the solution by 


deflection and moments are more accurate by | analysis of the 


' Comparison is also made herein between the results as obtained by anal- 


_ ysis of the plate model and by grid analogy!3 presented by Walter W. Ewell, : 
_ F, ASCE, Shigeo Okubo, and Joel I. Abrams, A. M. ASCE. In case of s—eare 
_ plate built-in along the four edges and uniformly loaded, the deflection of a 
: grid joint nearest to the center of the plate in a grid with 9 by 9 divisions is 


approximately 2. 71% in: error. The analysis involves the solution of ten simul- 


__ 13 “Deflections in Gridworks and Slabs,” by Walter W. Ewell, » Shigeo Ol Okubo, and Joel 
Abrams, Transactions, ASCE, Vol. (117, 1952, Table 3, p. 887. 
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PLATES 
TABLE 7, .—COMPARISON OF ACCURACY OF RESULTS OBTAINED BY ANALYSIS 
OF PLATE MODEL AND BY FIRST ORDER FINITE “comic 


PROCEDURE 


Item and 
Procedure 


5) 


Center Deflec- 

Plate Model J 2 0, 
9. 0443 | 0.0443 


Center Moment, 
“4 Plate Model 
Finite Dif- 


0477 
0.0472 |0.0479 


Plate Model | — .02652 | 0, | 0.02301 
Finite Dif- 10 0.02329 


“Mid -edge 
Plate Model 


Finite Dif- 


Uniformly Loaded Rectangular Plate Simply Supported Along 
‘Three Edges, the Fourth Edge Free (u= 


Max. Deflec- 


“Plate Model 0 0737, — | — | 0.0776 |0.0775 
Finite | 0.0739 - - 0. 0775 
Max. Moment, 
Plate Dit 0.062 | 0, osos | 


Finite 0.0567. 


taneous ote equations obtained after a lengthy process of moment and torque — 
of the the error in extrapolated value 


tion “Exact” , 
il 
— 
0,04875| 0.04793] — | —_ 
0.04063] 0.04557 — | — } 
Uniformly Loaded Square Plate Built-in Along Four Edges (= 0.30) 
18 
— 
37 
| |_0,04455|-0,0482 | — | -0,0512} -0.0513 0.2 
-0.0386 | — 0.0475 -0.0505} -0,0513 1.56 
Et3 
— |0.5 
“ps 


set of three and another set of six simultaneous linear 


of the center deflection ‘is s only 1.24%. “The : the 


CONCLUSIONS 


The results of analysis of the structural plate model a ey for 


the ‘purpose of solving problems of laterally loaded a plates justify — 


(1) The procedure is ‘It is applicable cases plates with 
supported, built-in, fgee, and elastically supported edges, as well as con-_ 
"tinuous plates. Plates with point supports, though not iil the investi- 
(2) ‘procedure is simple. ‘Itis only necessary to replace the plate by 
a structural model with a small number of intersection points. For example, 
it is sufficient to divide a uniformly loaded square plate, with edge conditions 7 
_ other than built-in, into 2 by 2 and 4 by « 4 divisions and extrapolate the results. is 
In case the edges are built-in, it is recommended to choose more divisions, — 
for example, models with 4 by 4 and 6 by 6 divisions, in order to describe the 
elastic surface of the plate with enoughaccuracy. The number of simultaneous ~ 


bs: equations involved inthe solution are thus reduced to a minimum, and answers 
are Obtainedinashorttime, 
_ (3) The procedure yields results with a high degree of accuracy. In the 
problems solved herein, errors as low as 0.13%and 0. 2%are obtained, and the 
highest error is 1.58%. Results of analysis by use of the structural model are | 
more accurate than those obtained by the finite difference method. Investiga- 
_ tion of the built-in square plate justifies the same conclusion, althoughthe num- — 
ber of divisions in the analysis of the structural model is smaller than in the 
solution byfinite differences, 
(4) ‘The procedure c: can be extended for use by means of f analogies in many 


& 


gy Laplace’s or Poisson’s equations, as well as plane stress and plane strain ; 
probleme usually sol solved ed by introduction of astress function, 
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—~s a FURTHER STUDIES OF THE STRENGTH OF BEAM- -COLUMNS 


proximate interaction curve expressions, this paper r presents the results of a 
continuing investigation of the ultimate carrying capacity of beam-columns, 
_ Two basically different loading conditions have been examined, and the results 
have been presented in interaction curve form, 
_ ‘Tt The first of the assumed loading conditions is that of a column subjected — 
to an axial thrust and end moments My and B Mo. The values of f that have 
been considered are +1.0, +0.5, 0, -0.5 and -1.0. The maximum strengths of — 
columns subjected to lateral loads at their  center- line sections h have 
In all cases it has been assumed that the member is an “as-rolled” wide- 
flange shape that is subjected to bending moments in the plane of its web. It 
., further been assumed that the member is constrained to deform only in 


the plane of the web. Failure is therefore characterized by ‘ “excessive bend-_ 


ing” he of the applied moments. 


type of. deformation ‘the possibility “of -““unwinding” and thereby failing in 


Note. —Discussion open until | To extend the closing date one month, 

a written request must be filed with the Executive Secretary, ASCE. This paper is part 

: of the copyrighted Journal of the Structural Division, Proceedings of the American So- a 


ciety of Civil Engineers, Vol. 87, No.ST 6, August, 1961. = = | oS 
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INTRODUCTION 


A structural member subjected to an axial thrust plus an increasing ‘magni- — 
tude of bending moment, due either to end- moments or lateral loads, will <i 


"will be reached. Thereafter, the deformational response of the member in- 
creases atan ever-increasing rate. _ The maximum carrying capacity is reached 
when the ‘member continues to deform no increase in applied loading. 
Analytical solutions to the of determining the strength of such 

-members have, in the strictest sense, followed one of _ approaches. The 

= 


= 


Lae 


FIG. 1.—CONDITIONS OF LOADING 


first of these is based on the reasoning that a lower bound to the actual col- 


= lapse load can be realized by using the elastic ‘limit value. Since the mathe- 

- matics of obtaining such a solution are straightforward, there is much to 
recommend this procedure. However, certain fundamental objections to such 
an approach exist. First, the exact amount of the reserve in 1 strength above wy 
the predicted proportional limit load is a | variable, ‘unknown quantity. Second, 7 
‘most “as- rolled’ ’ and all “as- welded” members contain residual stresses. 7 


of an abnormally, low value for the range of proportionality would therefore 
be required. For that case in which the maximum “locked- in” stress equals 
the yield stress (a condition that can exist in certain built- -Up, welded mem- q 


since the. theoretically “failed” prior to load ‘application, 


single curvature mode ex this situation + 
discussed and the pre Spossibility. 
if 
— 
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5 7 The second method of solution re requires the determination of the actual load- 
» ‘deformational response of the member to animposed, increasing loading sys- 
tem, Curvature values, which in the inelastic range depend on both the magni-— 
tudes | of the axial thrust and the bending moment at the section in question, 
are graphically or numerically integrated to define the true deflected onli . 
‘The ‘bending moments, however, depend o on these deflections. A trial- and- — 


an assumed deflection is therefore usually vena. 3 This general approach, fe 
_ originally proposed by T. von Karman, Hon. M. ASCE, has much in its favor, — 
since the true problem is considered. It should be noted, however, that solu- © 
tions by this method are time consuming. g. 


2. RESIDUAL STRESS PATTERN 
To facilitate the the and specification writer, considerable 
attention decade has been focused on the development of ap- 


_ 2 “Plastic Deformation of Wide-Flange Beam- Columns,” by R. L. Ketter, E. L. Ka- 

- ‘minsky, and L. S. Beedle, Transactions, ASCE, Vol. 120, 1955, p. 1028. Bate: 

aa > “Columns Under Combined Bending and Thrust,” by T. V. Galambos and R. L.. 

€ Ketter, ASCE, Vol. 85, No. April, 1959. 
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proximate “interaction equations.”4,5,6,7,8 In general these have been ob- 
_ tained by fitting “acceptable types of curves” to available test results. Com-— 
parisons were also made against the few “exact solutions” that were available. — 
purpose “of this paper is to present the results of additional computa- 
_ tions of “exact solutions” of ultimate carrying capacity and thereby afford . 
more complete spectrum against which proposed, approximate, oe al 
equations can be compared. Two basically different types of loading condition — 7 
_ have been studied and these are shown diagrammatically in Fig. 1. In Fig. l(a), 
= a column of length L is presumed to be subjected to an axial thrust, P, and be 
| end moments, My at the upper end and B M, at the lower end. The values of 
B that have been considered are +1, .0, +0.5, °o, -0.5 and -1.0 in which 8 pro- ; 
: ducing a single curvature type of deformation is considered positive. oo 7 
interaction curves for the cases 8 = +1.0 and 0 are reproduced from another _ 
publication.3) The maximum strengths of columns subjected to lateral loads» 
applied | at their center- line sections (Fig. have also been determined, 


s 
: study. “It is assumed that bending moments are applied about the major | en 
of this cross-section and that deformations occur only in the plane of the web— 
bs of the member. Failure, as considered herein, is therefore ) characterized by 
“excessive bending” in the plane of the applied moments. = | =. 
The stress- strain properties of the material are presumed to be ideally 
-elastic- ‘plastic; that is, there is initially a linear range where stress is pro-- 


portional to strain, followed by a region of constant yield stress. The ¢ cross-_ 


been chosen equal to 0. in which oy is the yield stress of the material, 
This is consistent with — of cooling residual stresses as — 
in the histogram9 of Fig. 2(b). 
i Having assumed the magnitude and distribution of the residual stress aati 
the stress- strain properties of the material, it isnow possible to define for a 
given | cross-sectional shape the basic relationships that exist in the inelastic 
range between thrust, bending moment, and curvature. This was done else- 
; where2 for the 8 WF 31 shape and those relationships will be assumed in this 


Enda It should be noted that the Le M- @ relationships a assume that the full full: 


4 “Strength Analysis of Eccentrically- Loaded Columns, ans,” ” by F F.R. . Shanley, Report 
54-57, -57, Dept. of Engrg., Univ. of California, May, 1954. 
5 $5 “Lateral Buckling of Eccentrically Loaded I- and H- Section Columns,” by H. 
“Hill and J. W. Clark, Proceedings, First Natl. Congress of Applied Mechanics, ane, 
= a: “Recherches sur le Flambement de colonnes en acier A 37, a profil en double Te, © 
sollicitees obliquement,” by by F F. Campus and Massonnet, Rendus de Recher- 
ches, IRSIA, April, 1956. 
ee The Stanchion Problem in Framed Structures Designed According to Ultimeate _ 
"Carrying Capacity,” by M. R. Horne, Proceedings, Institution of Civ. Engrs., Part II, 
Vol. 5, No. 1, April, 1956, p. 105. 
_ 8 “Guide to Design Criteria for Metal Compression Members,” ” Column Research 
Council of Engrg. Foundation, Ann Arbor, Mich., 1960. 
_ 9 ©The Influence of Residual Stress on the Strength of Structural Members,” by R. L. 
—_— Proceedings, Seventh Tech. Session, , Column Research Counci Council, (1957. 957,00 == 
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magnitude of the axial thrust is 5 tiret applied and while it is held constant, © 
bending moments are increased to their maximum valves. A definite loading 
Elastic Limit Solution. —For completeness and for future comparisons, it 
2 desirable to include at this stage, before discussion of the development of — 
the ultimate strength curves , the elastic limit solutions for the first of the two _ 
loading conditions shown in Fig. In Appendix I, it is shown that for the case 
where no residual stresses are present, the following | algebraic interaction 


in which 


a “only of the ratios . and B. “The corresponding curves are given in Fig. 3. ; 
To obtain the curves shown in ‘Fig. 3, it is first necessary to determine — 

for given: values of and the location of the section of maximum 
ment (and therefore maximum stress). This is obtained from Eq. 3. Knowing» 

this location; ; that is, knowing x/L from Eq. 3, the modifying or “amplifying” 
factor, y ‘is obtained from Eq. 2. If Eq. 3 has no solution in the range 0 < x 
_<L, the maximum stress occurs at the end of the member and the amplifica- 

tion factor equals one, 

if Mo » = 0, two possible solutions exist. . The first of these follows Greely 


“second possible solution corresponds to that case in which the amplification 
factor, approaches infinity. Here, the maximum axial thrust is the Euler 
buckling load for the pin-ended member. 
_ Ultimate Carrying Capacity Solution.—The method of solution used in the 
investigation of the ultimate | carrying capacity of a given member and loading 
a" the same as that used in another paper. 3 It is essentially as follows: The 
‘member is sub-divided into eight equal segments as shown in Fig. 4(a). A 
constant value of axial thrust is assumed. In addition, a value of the end mo- 
ment, Mo, is also selected for investigation, Deflection values, based on a 
7 _ consideration of the elastic response of the member to the assumed loading 
or on information obtained from previous inelastic calculations, are assumed 
at each of the vi various ‘division points along the member. Using these, mo- - 
‘ments at the several points” are computed, and from the M- -P- relationship 
‘domme previously, corresponding curvature values are established. If 


these curvatures are numerically integrated twice, new deflection values are 
ty These will not necessarily be equal to | those originally assumed and 


7 ‘The previously described process would be repeated for successively | 
— values of me end moment. ‘This wale make possible the definition of 
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FIG. 3.— ~ELASTIC ‘LIMIT INTERACTION 
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5. —DOUBLE CURT LOADING 


— 


the load-deformational response of the einaia e.g., the relationship between © 
_ the end-moment and end-slope could be plotted as shown in Fig. 4(b). From 
_ this diagram, the ultimate moment (for that particular chosen length of mem-_ 
7 ber and assumed axial thrust) would be determined. 
_ The total process described in the preceding paragraph would be , repeated a x 
for other assumed values of constant axial thrust. For each of these there r 
_ would be corresponding critical values of the applied end moment. A curve of 7 
_ these solutions could be plotted in which axial thrust would be the ordinate 
ane the ultimate end-bending moment would be the abscissa. Such a curve will - os 
be referred to as an “ultimate strength interaction c curve” and wih te 
form for presentation of the results of this study. 
_ With the exception of that particular case for which B = -1.0 (the double 
curvature case) the numerical process described previously for solving the 
‘ultimate . strength problem is stable for all values of end bending moment less 
than ultimate. For that particular loading however, the possi- 


FIG. 7. _ULTIMATE | STRENGTH WINTER 
ACTION CURVES ( = +0.5) 
bility exists that at certain lesser values of axial thrust or end bending mo-— 
= or both, the deflection apes mand _— into one of single curva- 


by examining, at any stage of applied the response of the 
7 ‘ber to an infinitesimally small first mode deformational disturbance. For an 
— assumed ideally elastic, pin-ended member, a shift in the location of the 
= of zero moment and the beginning of “unwinding” | will occur when the a 
can ‘magnitude of the applied axial thrust reaches the “first mode- Euler buckling | 


“2 
~The “real member” in question will not “fail” in the elastic range. For the 


actual situation which the member is partially yielded, it is necessary to 
examine the response of the member to a single curvature type of virtual 
disturbance at each increasing value of | M, o- (This process is described in 


e142 August, 1961__ 
by 
| 
| 

a 


> 


 BEAM- COLUMNS 


detail in Appendix Ve where sample calculations are given. 1.) In certain ranges, 
the member will unwind prior to reaching the “double curvature mode— 

- excessive bending—ultimate carrying capacity.” This is illustrated in Fig. .-' 
_. The ultimate strength interaction curves for the cases in which the member 
is subjected to thrust the combinations of end bending moments 
n in Figs. * 4 to 11. For that situation where the 
‘member is eubdorted to a centrally applied lateral load in addition to the axial 
thrust, the solutions are shown in Fig. 
In the strictest sense these curves. apply only to 8 WF 31 members that 
contain the residual stress pattern shown in Fig. 2(b). Furthermore, the yield | 
- point of the material was assumed to be 33 kips per sq in. With regard to the | i : 
first of these restrictive conditions, it may be said that the 8 WF 31 has one 7 7 
7 of the lowest of the shape factors. It would therefore be expected that strength - 
_ predictions on the basis of these solutions should be conservative. For mem- 7 
bers having material ‘properties: essentially the same as those of the ideally 


9. -— ULTIMATE STRENGTH INTER- 


7 elastic- plastic ty type, ‘but having « a L yield : strength other than 33 kips pe per sq + sq in., 
the slenderness ratio will need to be modified to allow the use of these * charts, 
: A modification according to Eq. 5 is considered to be adequate. a 


COMMENT 


It is ‘not the intent” of this paper to > justify the adequacy of the mumerical 
integration procedure, Therefore, comparisons with test will be 


— 
| 
“5° 
= 
| 


it should noted, however, that such studies have b 
the cases B = +1, 0 and 0. 
Another publication!2 discusses the various approximate interaction ex- 


pressions that have been and are now considered. For the the case where 
B = +1,0, two general forms are listed:5,6,8 


Myit/ 


te 
elsewhere, 10 of these adequately represents the 


- 11. —ULTIMATE STRENGTH INTER- ad FIG. 12. _ULTIMATE STRENGTH INTER- 


ACTION CURVES (B = -1.0) ACTION CURVES ( CONCEN- 


For B 9 1. 0, it has been suggested that an ‘ “equivalent uniform moment © 


ct or” be defined. Based on considerations, ‘Massonnet, 


10 “Stability Considerations in the Design a6 Steel by C. Pro- 


ceedin; , ASCE, Vol. 85, No.ST 7, September, 1959. 
Bi mw... SCE- WRC Commentary on Plastic Design - - Compression Members, ” Proceed-— 
s, ASCE, Vol. 86, No. EM 1, January, 1960. ee ee 
pe 2 “Guide to Design Criteria for Metal Compression Members,” ’ Column Research | 


Council of Engrg. Foundation, Ann Arbor, Mich., Chapter 5. vi 


Design in Steel,” A.LS. C., New York, N. Y. 


— 
— 
— 
FER 
— «CLC 
ig 


oe 0. 29 and 0.42 
"place of 0. 3 and 0.4, respectively, in Eq. 8. ‘Horne, using a somewhat different — 
approach has obtained? results similar to Massonnet’ 
Forthe values considered herein, Table 1 gives!4the modification factors 
tt that correspond to: each of the studies referred to previously. ey ee 
- For the case in | which the maximum moment occurs at the end of the )mem= 
ber, has the following ‘second, ‘interaction 


y ‘should be noted that a double standard 1 has been proposed, 


3 


which the equation predicting the least M)/Mp i orrect for that particular — 


TABLE 1. _MODIFICATION ‘FACTORS. 


ream bending moment factor for the “as-rolled” 8 WF 31 section -. 
_ strained to bend in the plane of its web. It should be observed, however, that 
3 such comparisons | contain four dimensionless parameters instead of the two 


possible to combine the axial slenderness terms and 


thereby “Obtain; a form more— easily compared w with the elastic 


Three values of P/P, have been selected for illustration: ome: = 0.2, (0. 4 -" 
6. In Fig. (13 B is plotted as the abscissa, MEqu/Mo as the ordinate, and 


__ 14 “Guide to Design Criteria for Metal Compression Members,” ‘Column Research 


Couneil of Engrg. Foundation, Ann Arbor, Mich., Table 5.3, 80. 
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FIG. 14. o EQUIVALENT MO} MOMENT FACTORS WITH MASSONNET’S. 
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ae = standard is required for the full range of definition of the equivalent 7 
moment, then the ratio P/Pr is a major factor in determining the equivalent re. 
moment expression. It is as significant as the factor 8. Of somewhat lesser 
importance is the influence of P/Py. It should be noted, however, that 
P/P, becomes larger, its influence is ‘more pronounced, 
double standard is used, the region of influence of ‘the ratio can 
reduced. _A comparison of _Massonnet’s equivalent ‘moment expression 
(Eq. 8) with the values - given’ in Fig. 13 is shown in Fig. 14. Only that range . 
wen requires the use of the equivalent moment has been included. That is, 
_ for smaller values of P/PE and f, the strength of a given member as pre- 
dicted by the Massonnet procedure would be governed by Eq. 9 rather than by Any 
6 and 8, ‘This region has been excluded in the 


APPENDIX, DETERMINATION OF ELASTIC LIMIT INTERACTION 
EQUATION (NO STRESS) 


o ‘deflection at t any section “x” along the length of the member is 
sae = Mo [Sin k(L- 


kL={=n 
Pr 

The moment at the corresponding section is then 


Since the maximum stress will be reachedat the section of maximum moment, a 
it is first neces essary to define the location wad the critical section. That is 


= -2 + ( 


a “Ts “Theory of Elastic Stability,” by S. Timoshenko, New York, N. Y., McGraw-Hill _ 
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GENERAL LOADING C¢ CONDITION on 


APPENDIX - —INSTABILITY DUE TO UNWINDING” 


= ‘The seine and loading condition in question is that shown | in Fig. 16 (b). : 
a The problem is to determine that particular value of (M,/M ) for which the 


To solve this problem it will be assumed that at the stage of loading in 
question the member is subjected to a virtual, deformational disturbance —- 
the single curvature - type. The tangent to the moment- -curvature relationship _ 
for any particular value of moment represents the instantaneous resistance 


of the member to an increasing moment, Where moments decrease, the - 


information, Newmark’s procedure for determining the buckling load can be 
_ directly applied. 16 “Unwinding” will occur when the critical axial thrust thus us 


computed equals that applied to the member. 


member, ‘it. is therefore possible to ¢ define the effective stiffness. With 


a 


_ 16 “Numerical Procedures for Computing Deflections, Moments, and Buckling Loads,” o 
ASCE, Vol. 108, 1943, p. 1161. 


> a Substituting Eq. 16 in Eq. 14 give 
|. | 
a ~ The interaction relationship is tht 

— 
¥ 

| q allow a in the po ero moment (initia med to be at the 
q 
: 
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Fig. 1 12400 indication mn of the degree of accuracy obtained. ‘The member 
in eo: is the same as that considered in Fig. 16. The end meamend is ead 


6. —DETERMINATION OF THE CRITICAL LOADING ASSOCIATED WITH : 
UNWINDING 


ai 


umed Deflection (inches) 


FIG. 17.—INDICATION OF THE DEGREE OF ACCURACY USED IN THE 7 


st stability. Fig. 18 summarizes the numerical 


— 

rp 
0 0005 = 
| 

0.257 0.408 0469 0459 0.392 0.263 0.148 0 Computed Deflection (inches) 
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«See ratio, the eight division spacing 1 Oy = = 33 kips per sq in, and E 
30.000 kips per sq in. 


Galeulated Deflection ) (20) 
we. 


sr considered: 


0.280 0672 0 952 1.000 1.000 I. 000 
om = 0.520 0.832 983 1.000 0.889 0.665 _ Assumed Deflection 


(EL) 
(From M-P-@ Relation) 


Concentrated Angle Change 
_ (Parabolic Assumption } 


Calculated Deflection 


__Assumed Deflection 


Calculated Deflection — 


‘FIG. (18. —SUMMARY OF THE ‘COMPUTATIONS REQUIRED IN THE DETERMINATION > 


- ‘OF THE CRITICAL THRUST ASSOCIATED WITH UNWINDING 


The Silitieass symbols have been adopted for use in this paper: 


rea of cross- -section; 


— 


= Young’s Modulus of Elasticity; 


at end of member; 


Mp fully plastic moment under pure moment; 


initial yield moment under pure moment; 
MEqu ~ = equivalent uniform moment in a beam- 


7 
a 
— 
Qa 
~ Therefore 
— 
— 
7 
19. 5 48 il. 10. 7894 4215 
19808 15.268 12.548 11.920 10555 7694 4250 
49.044 78280 92.248 93668 83.168 62.113 33.164 
9.01060 001062 0.01065 0.01067 001068 001070 001070 
if 
4 
7 
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= Euler buckling load 


= - radius of gyration; 


= non- ieee parameters; 

ratio of applied end moments (Fig. 1); 


— 


— 
— 
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STRENGTH © OF PL PLATE GIRDERS IN BENDING 
Konrad Basler,! A.M. ASCE, and Bruno Thiirlimann,2 M.ASCE 


- oan first the bending strength (this paper), then the shear strength (in al 
the October 1961 Journal of the Structural Division), and finally the interaction 


_ between bending and shear (in the October 1961 Journal of the Structural Divi- 


of the Journal of the Structural Division, Proceedings 0 of the American 
ciety of Civil Engineers, Vol, 87, No, ST 6, August, 1961, Cee ae — ne 
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| 
2 investigation of welded plate girders was conducted at Lehigh University 
oe mt cl during the years 1957 to 1960. The objective of this project was to determine | 7 
aay the static carrying capacity of transversely stiffened plate girders, The study " 
7 _ was grouped into an experimental and a theoretical phase, andthe results of [am 
7 > the experiments were published as a Welding Research Council Bulletin, i 
| Note.—Discussion open until January 1, 1962, Separate discussions should be sub-_ 
mitted for the individual papers in this symposium, To extend the closingdate one month, 
| 1 


1961 


adi 


n analysis is presented of ‘the static capacity of plate 
eo subjected to bending. "First, the phenomenon of “web buckling” is 
4: discussed to explain the discrepancy between actual behavior (1)3 and theory. — 
‘Pious this is an analysis of the bending strength ot plate girders. Finally, 
an attempt is made to interpret the findings as a design concept. 


as possible in the extreme fibers. By keeping the web area as small as possi- — 
ble, the lever arm of the internal forces is maximized and with it the carrying | 


A asa the Gotten: of plate girders the tendency is to arrange as much : material } g 
capacity. In the past web buckling has been set as a limitation to this tendency * . 


towards an optimum utilization of the material. Consequently, considerable > @ 

effort has been spent in establishing web buckling values (2),(3),(4). 
conventional plate-buckling theory predicts the load intensity under 

which a plane plate subjected to edge stresses deflects on of its plane. The | 


the same e word “buckling” is used to ‘describe the phenomenon in a plate. 7 
_ Because the computed column-buckling value gives an adequate measure of — q 
the strength of a column , it was natural to consider that a web-buckling value —_ 
is of equal significance. Such an assumption is not true. The strength of a 
plate can exceed its buckling limit, and this additional 1 margin of strength is 
termed “postbuckling strength” 


— 

6 

® Numerals in parentheses refer to cor 


= 


one: of safety was prerne to web buckling than that used for primary col- 7 
umn buckling. What should be the value of this factor of safety? In order to 
Clarify this problem the plate- girder investigation a at Lehigh University 
(Bethlehem, Pa.) was started, 


— Cmpersanenee (1) demonstrated that the concept of expressing the post- 
a girder certain percentage of the web- 


where they first in the illustrations or in the cont, and are 


alphabetically, for convenience of reference, in Appendix Il. ‘cael 


ty 
WEB BUCKLING © 
the purpose of this section to to explain the b behavior c of a plate: when 
_ Phenomenon of Web Buckling.—A rectangular plate subjected to edge hia 7 
is shown in Fig. 1 (a). For a given edge-displacement pattern the strain €g -_ 
expresses the applied deformation and Y is a nondimensionalized measure of ‘ 
_ the plate deflection out of its plane ata particular location. Assuming com- 
pletely elastic behavior of the material, the result of the linear buckling © 
‘theory that is commonly used can thus be illustrated by Fig. The plate 
remains plane’ for all positive values of the | applied strain éa (tension), te 
also for a limited range of negative strains (compression). However, = 


q 


this X-Y diagram bifurcates. The plate ‘must then bulge out of its plane. This” 

_ phenomenon is referred to as plate-buckling and the corresponding load is 

termed t the plate-buckling load or critical load. ee 
The actual behavior differs radically from that predicted by the linear = 

buckling theory because initial plate imperfections are present, With the help _ 


of Fig. this contradiction will be explained. In Fig. 1(c the relation be- 


Y, is again yeotted, assuming that the plate has a small initial deflection of 
the magnitude Yy as shown for X = 0. Discussing first the upper half of Fig. = 
1 (ec), , it is clearly seen that only the strain-deflection path of a perfectly plane 
plate exhibits a bifurcation point of equilibrium, and that the sudden rate of _ 
increase of plate deflection disappears with increasing initial deflections . 
_ These initial imperfections need not be very large to completely obscure the 

- buckling behavior. The majority of the tests conducted at Fritz Engineering 


Laboratory on shop fabricated, welded gy girders furnished web-deflection 


. 


have, in general, smaller initial Gutestess and may lead to a path r more e like 
_ the curve for Y; = 0.2. In aluminum girders the relative initial imperfections 
are even smaller due to the fact that the lower modulus of elasticity requires 
a sturdier “web to furnish the same > critical stress as an equivalent steel 


sie ss The existence of postbuckling strength has been pointed Out ever since [i 
— | 
= considers the influence of the flanges and transverse stiffeners on the overall 
capacity. The strength of girders with r t tobending isthe 
i; 
: 
— 


girder. Therefore, a more pronounced rate of nveiee of web deflections can 
be expected such as indicated (7) by the curve for Yj = 0.05. Curves which > 
come consistently closer to the singular case Yj = 0 seem to be obtainable 
only in carefully conducted laboratory tests on isolated aluminum plates (8). 
The conclusion is that, in general, a web-buckling load cannot be observed 
on a full scale, shop fabricated plate girder built of steel, because the trans- q 
2 ition from the pre-buckling into the post- buckling range is not accompanied - 
_ by a sudden increase in deflection. Thus, concerning the behavior of the web, — 
: there appears no reason why the computed web-buckling strain should not be 
exceeded in a plate girder under working load, 
~~ Fig. 1(c) also indicates that, incases where postbuckling strains are toler- _ 


ated, an effort to cold straighten the web of a slender girder in the shop — 


7 


several have derived similar | deflection curves for special 


in which X is €a/l€cr|, measure the applied deformation, Y; 4 
_ the initial plate deflection, and Yis the resulting plate deflection. The plate de- 
_flections Yj and Y are nondimensionalized in terms of the plate. thickness and 4 
a constant which depends on the shape of the deflection. Eq. 1 was derived (14) © 

_ by assuming that the deflected configuration grows under load without changing — 


its shape and by neglecting the influence of the membrane shear ‘stresses. qt s 


rx 

a 

— 

— = 
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PLATE GIRDERS 

, — be further pointed out that the branches resulting in negative values of na 
Y (lower half of Fig. 1 (c)) are also equilibrium positions for the plate under the é. 
stated assumptions. If the plate could be brought | to a position on the dashed ~ 
_ branches which are unstable equilibrium positions, , it could then snap toeither , 
side of the reference an and stabilize there. This is the so-called canning 

Web Participation under Bending .— 

ofa plate when strained beyond the buckling limit, an explanation is now given — - 
as to how the deflected web participates in carrying part of the moment applied — 
_ to the girder, Of all the stresses in the web only the membrane stresses in 

4 the longitudinal direction of the beam give rise to a bending moment about the | 

‘ girder’s neutral axis. Therefore, reference is 3 only made to these stresses or or 


equivalent strains in the subsequent material 
As an example, , measured web deflections and membrane stresses are 
shown in Fig. 2, plotted against the outline of the cross section and the ele- — 
- vation of the particular test girder from which the data were observed (1). ' 
The membrane strains were . obtained by averaging the strain readings on both 
7 sides of the web. The strains or stresses—depending on the scale used—are > 
_ plotted fo for three different bending moments, These moments are expressed © 


- yielding naneedine to ordinary ‘beam ‘theory. The test data corresponding to a 
particular applied moment are connected by straight lines. Thus, the stress_ y 
distribution over the entire girder depth can be visualized and compared with | - 

_ the predictions of beam theory, predictions which are also plotted in thin lines —_ | a 
each of the three moment 


The initial distortions in this cross section and the additional web deflec- 
4 tions 3 under the » three | ‘moment values appear on the left-hand side of Fig. 2. | 
4 How the stresses are related to the web deflections becomes apparent when 2 
longitudinal web strip, which extends over the panel length, is considered. It 
is readily seen that an axial shortening canbe “digested” not only by straining 
_ alone, but also by lateral deflection of the web. Thus, strips in the compres- 
; sion range (upper portion of Fig. 2) avoid carrying their full share of axial 
stress, as shown quite clearly in Fig. 2. 
‘This modified stress distribution could be estimated using the method on 
which the derivation of Eq. 1 is based. However, because it would depend on 
the initial distortions which are quite random, ‘an assumption as to the result- 
. ing stress distribution will be made rather thanto the shape and magnitude of 
the initial web distortions. It will be assumed that the contribution of the | 
= web portion may be disregarded except | for an effective strip 
along the compression flange (Fig. 3(c)). This leads ‘to the effective-width — 
as already used in existing specifications (15),(16). 
There exists little information on the effective compressive strip i in plates 
a aa to edge displacements similar to those of webs of bending girders. q | 
a 5 But it seems safe to assume that at ultimate load at least the same effective 
— strip width can be expected as for a thin plate under uniform — compres- — 
sion, namely about 30 t for mild steel. 
seen from th the previous material that the attainment of th web- 
buckling load does not constitut ute failure of a of a plate girder. When \ wil a a, 
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girder subjected d to bending fail? Obviously, not ‘before one of th — framing 
_ members around a web panel fails, because only then will the deflected web 

be unable to burden the framing members with that portion of the bending 
7 moment which it alone cannot resist. In a symmetrically proportioned plate © 
5 girder subjected to bending, the static carrying capacity is reached when the — 


flange fails—barring brittle fracture. Therefore, this 


be concerned with the bracing and proportioning of the compression — ] . 
Ih order to classify the different buckling types, the compression flange 
“may be considered as an isolated column. With regard to buckling, such a 
column has three degrees of freedom: it can buckle laterally, torsionally, or _ 
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: .- MEASURED STRESSES AND WEB DEFLECTIONS 


“show the three directions of ‘motion. In ‘Fig. to used later 
Vertical Buckling of the Flange. —If the compression flange possessed 

7 rigidity in all directions, it could be the lone carrier of the compressive force. 
Together with the tension-flange force this would provide the resisting bending» 
moment. This concept is realized in a truss. A plate girder of I-shaped cross © 

- section could act the same way provided the web were able to brace continu-_ 

ously | the compression flange which otherwise lacks rigidity in the web 

Ss | the required bracing stiffness is small, the danger of 
compression-flange failure in the vertical direction is limited to high web 

ae ratios. By setting an upper limit to the web slenderness para- 
meter f = b/t, it should be possible to arrive at a design concept which is un- 
affected by this type of failure. To establish such a limit will iia cial a num- 


— 
7 
— 
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the process of bending, curvature gives rise to transverse flange- 
force components which cause a uniform compressive stress oy on the upper — 
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PLATE | GIRDERS 
- "and lower edges of the | web. . This is indicated in Fig. 4 aus the reaction 
ee On needed to keep the flanges in equilibrium are shown. If the plate 
i were subjected to stresses of this kind only, failure of the web would occur 
similar to that of an Euler column a 


Actually, the plate is also “subjected to 1 longitudinal stresses as pictured in 
Fig. 2. Assuming that the adverse influence of the compressive stresses in 
_ the longitudinal direction is offset by the tensionin this direction, and that the 
influence of yield penetration into the web at high flange strains i is offset by 
zg restraint offered to the thin web by the flange, the value given above can 
be considered as an estimate of the web’s resistance against vertical buckling | 


FIG. 3. _ NOTATION, (b) BUCK LING MODES OF THE COM- 
PRESSION FLANGE, AND (c) EFFECTIVE GIRDER 
CROSS SECTION =— 


Over the length the transverse component of the force amounts 


to Ag of dx, where o¢ is the flange stress and €¢ the strain. 


vertical buckling it is required that the .e applied force be 


¥ is re uired that ev ever flan je fiber ‘Teaches ield stress | before failure, ‘then a 
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a 7 This slenderness limit depends somewhat on the ratio of web areatoflange § 
area, Ay/Ag = p. In general, this ratio of area is not below 0.5. In order to 
prescribe an “upper limit for b/t applicable within the range of practical 


1961 
flange force Ag @. @y is obtained. To fulfil this condition the flange strain mu: must | ‘ Ti 
_be somewhat greske r than the yield strain €y= dy/E because an elimination of 
_ residual stresses 0; . requires a strain of e¢= Oy + or)/E, as sketched in Fig. Be 
Under these. assumptions and by numerical values, Eq.3 


(4) 


and amounts to b/t < 360 for mild steel where oy = 33 ksi and dy is 
_ Although this derivation for web- slenderness limitation appears somewhat _ 
it agrees fairly well with the test observations, points out the p para- 
ane } involved, and establishes a relationship between them. Eq. 4 indicates 


= 


that a ‘siete built of high-strength steel will require a lower b/t limit than 
_ one made of mild steel. This is due to the fact that the flange force and curva-— 


4 ture can increase in proportion to the higher yield stress, thus giving rise to 
a larger vertical force component. For a girder built of ow-alloy high-strength 


steel with a minimum yield point of 50 ksi and 0, again assumed to be 16.5 
the web-slenderness ratio should not exceed 250. 
_ This vertical failure of the compression flange was observed (1) in a test 
a girder whose > web-slenderness ratio was - = 388, web-to-flange area. 
ti Aw/Ag = 0.68, and whose flange strain €¢ was slightly more than the 
yield strain experiments also included a girder with a tubular com- 
pression flange, a cross section of which appears in Fig. 2. It differed from 
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mentioned previously only in the the shape of the compression 
Regardless of the amount of strain imposed, this girder could not be made i 
to fail: in this manner. In this case the flange was a self-supporting element 
_ _ with respect to vertical buckling and needed no bracing by the web. There- | ; 
_ fore, it can be concluded that the stated limit for the web-slenderness para- 
_ meter is conservative for all those girders whose flanges provide a certain , 
degree of vertical rigidity, _ such as in riveted girders with flange angles. _ 
_ It must be pointed out that girders which are loaded not only at points of 
_ transverse stiffening, but at intermediate points as well, require webs s sturdier 7 
‘ than stipulated by Eq . 4. This will also be the case in curved compression. 
flanges, such as those sometimes used over interior supports of centimous 


a. Lateral Buckling.—Lateral buckling of beams has been analyzed thorough- 
y (16),(17). The topic is taken up again only in order to derive and justify a 


"design concept especially applicable to plate girders. 


In an effort to arrive at a design formula the paper by K. de Vries 


de Vries, which depends only on the parameter of (span 
x girder depth) - (flange width x flange thickness). This formula, 

appears in the upper left corner of Fig. 5. 

It is the objective of this section to show that the ; second formula of the 


column whose effective cross s section is composed of the ‘compression ae 
and one sixth of the web. ‘This concept is pictured in the upper right corner | 
ae Relation Between St . Venant Torsion and Warping Torsion.—Since the p pa- 
_ pers by Timoshenko (2) were published inthe early 1900’s, it has been. recog- 
nt nized that the resistance of an I-beam against lateral buckling consists of two 
parts: St. -Venant torsion and warping torsion, often referred to as pure | tor- 
sion and flange bending torsion, The St. Venant part is due to twisting of each 
component plate where the angle of twist gives rise to a shear stress flow as 
7 pictured in the left of Fig. 5. The sum of the torques due to this stress flow 
is the St. Venant torsional resistance considered in the accepted analysis of 
lateral-torsional buckling of beams. The warping contribution, shown on — 
right side of Fig. 5, is due to lateral a of the flange plates and will be 
‘discussed in detail subsequently. 
_ In order to reflect these two parts, the expression for the critical bending 
- ‘moment is rewritten between the two sketches of Fig. 5 as expression o 
The assumptions which it ean be transformed into expression (b) of 
Fig. 5 are shown in Appendix I, where the two approximate formulas appear - 
as the equivalents of these two terms. The first term takes only St. Venant 
torsion into account, the second only warping; hence the symbols Sor (v) al 


three | the that between the sides of a right tri- 

angle. The length of the hypotenuse represents the correct critical stress, 

"whereas the two sides indicate the prediction by each formula alone. ‘There- 
fore, a conservative estimate: of the lateral buckling stress is is obtained if 
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7 4 - it was pointed out, especially by G. Winter, that there should be not one buta [im 
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either 
 Oer(v) or Sor te) predominant, the other can be neglected, since the 


of the hypothenuse is only slightly more (in no case more than 41%) than that — 
. It has been (18) convincingly demonstrated that for hot-rolled beam sec- 
tions the o,y(v) expression applies. However, in plate girders the aseaon 
torsion is generally the governing factor and the design should be based on 
_ the aforementioned column concept. This is seen from Fig. | 6, where the 
critical buckling stress of a plate girder is plotted against the slenderness 
&§ ratio of a column whose effective cross section is composed of the compres- 

- sion flange and one-sixth of the web. By assuming that the flange thickness is 
ee related to the web thickness, the stress resulting from the formula Sor (v) and 
ie _ the exact critical stress ocy(v, w) can be plotted in the same coordinate sys- - 

This is done for a flange to web thickness ratio of 3 and various web 


tem 
ratios. As seen, the _ buckling stress 


St Venant Torsion 
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5. LATERAL BUCKLING 
Cerlv, w), considering both warping and St. Venant torsion, exceeds the simple 


prediction | only for a web ratio B= 200 and even 


>= 


“Design Considerations: Test Results” the 
- program was confined to a low range ¢ of column slenderness ratios, from Vr 
15 to 40, that is, A=0.16t00.42,0 
Physical Interpretation of Warping Torsion. — Analyzing the mode of buck- 
— ling from the lateral- torsional buckling equations, it can be shown that, for — 


the case e of St. torsion, ‘rotation of the cross section 


mode of a beam under pure bending and with supported ends: 
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a 
which | Wo | is on maximum deflection of the top flange at centerline of the 


x, y) = 


span and normal dics web plane. The lateral force component of the — 


pression flange ii is ‘equal to 0 


. This overturning force is held in iti 


CKLING STRESSES FOR GIRDERS WITH SLENDER 


Inserting w(x, dad b/2) in this differential equation yields 


r=¥ Thisis is simply the lateral- buckling stress of the ‘compression 
flange plate. The t tension flange neither a stabilizing nor a detrimental 
: effect since it remains s undeflected in the lateral direction and, as stipulated, 
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St. Venant torsion. is neglected. ‘However, the stresses in the web also create 


lateral force component which amounts to 
per unit depth and is pictured in Fig. 7, Lacking rigidity in the lateral direc-— 


tion, the - certainly burdens the two flanges w with the resultant lateral forces : 


FIG. ihe LATERAL FORCE COMPONENTS OF WE! 
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£2.24 


4 
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deflection. To the lateral (Eq. 9a) and the equilibrium condition 
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vey 


mate 
“which is the” expression derived in Appendix I. should t be noted that the 
am critical stress represents the stress at the centroid of the compressed ele- 
ment furnishing the lateral r rigidity, thus” it is not the extreme fiber stress. * 
_ From the previous derivation it can readily be seen what modification of 


Eqs. 10 would be needed ifthe stress distribution over the girder depth =s 


from the assumed linear one. For example, a stress distribution similar to 
' ¢ that of the full plastic moment necessitates taking into account one quarter of 
the web area instead of one sixth. ._ Fora stress distribution as pictured by the 
e plotted points of Fig. Be a value smaller than Ay/6 would result. ‘Since the 
- influence of the webon lateral buckling is small compared to that of the flange, a 
_ it is justified to use the average value of one sixth, corresponding to the unesr 
stress distribution, also for the other extreme cases. 
_ If the girder cross section differs somewhat from the doubly symmetric 
7 I-shape, the neutral axis may shift from the mid-depth of the web. In ‘ae 
- girders this effect is generally limited in extent, and thus it is not necessary 


to change the equivalent web portion of b t/6 = “Ay/6 to a more precise value. 
Of much greater importance for profiles without double symmetry is the fact 
. that compression and tension flanges have unequal lateral rigidities, an ex- 
ample of which is shown in Fig. 7(b). From the physical interpretation it be- 
comes clear that in this case of negligible St. Venant torsion only the lateral 
: rigidity of the compression flange is of significance, Therefore, Sor (w) as 
_ given in Eqs. 10b and c should apply fairly well to all plate girders with sym- 
metry about the vertical axis. To check this, reference can be made to an_ 
excellent survey on simplified lateral buckling formulas recently assembled © 
(17). In Figs. 5 to 8 of that paper the lateral-buckling stresses as obtained by 
various approximate formulas applied to o unsymmetrical I-sections are com-— 
pared with the rigorous solution (Eq. 1 of that paper). This comparison indi- | 
cates that the formula Cer does not provide a good approximation to the 
exact value for the case of 1 = = 300 in. According to the previously discussed 
relation between Scr (v) and p(w), it should be expected that o,;(w) would 
ee more adequately. That fhis is the case is shown in Appendix I, where | 
| “Figs. 7 and 8 of the previously cited paper (17) are reproduced as Fig. 13. 
‘Because the phenomenon of lateral buckling if St. Venant torsion is neg- 
lected is simply one of lateral buckling of the compression flange, the buckling 


urves used in the inelastic range must be those of weak axis buckling of 
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_ flange type columns, It can be concluded (20) that, depending whether t the 


girder is welded or riveted, the critical stress in the inelastic range must — 

- significantly different. In order to avoid too great a differentiation, itis — _ 
4 recommended (16) that the basic column curve be used, which represents an = 
average of the extreme transition curves ina stress versus slenderness dia- 


FIG. 8.— LATERAL BUCKLING CURVE, 


column curve be used for later: lateral — of plate girders ernie with the 7 


€y 


2 I 
‘The curve of Eqs. 11 is plotted in Fig. 8. First the standard slenderness 7 
_ ratio ¢m , which makes the plot independent of yield stress, _is used for the 
abscissa scale, Next the scale for the parameter 1/r, valid only for éy =9 y/E 7 


= 33/30, 0.001 1, is shown. Finally, the old established parameter “puck- 
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cable only if the shape of the compression flange is a “ae. i eye 
_ Preservation of Shape.—The derivation of the general lateral buckling ex- 
‘pression (Eq. 23 in Appendix I) is based on the assumption that the shape of | 
the girder cross section is preserved at the instant of buckling. The validity 
this assumption | becomes u uncertain for plate girders with high web- 
“slenderness ratios and in cases where the transverse stiffeners are not fitted 
to the tension flange. However, this uncertainty concerns only the St . Venant — 
torsion | part since it is dependent on an undeformed cross section whose com- 
«ponent plates are assumed to be forced through the same angle of twist. If, 
for example, the joint between the tension flange and web were pinned, the 
torsional constant K in Eq. 23 would reduce to about one-half the value for a_ 
rigid joint, (exactly one half if the web’s contribution in 5 b t3/3 is neglected), = 
a reduction in “critical” stress, of almost 30% would result. The 
derivation. of Scr (w) in Eqs. 10, however, is not affected at all. Since —— - 
torsion in a a plate girder is the predominant contribution to lateral stability, a 
deformation of ‘Profile shape has but little effect on the resulting buckling 
and an on the basis of — alone isa conservative one 


‘direction is needed (Eq. ger having shown that a preservation of the right 7 
angle between tension flange and web is not essential, a theoretical confirma- 
tion is obtained of the test evidence that transverse stiffeners may be cut 


short for a limited distance at the side of the tension flange (21), > _ 


If all restraint on the flange from the web is neglected, the case reduces 
* buckling of a long, hinged plate under pure edge compression at its ends. ; 
- Hence, the only parameter on which the flange plate buckling stress depends” 
is the ratio of outstanding width to plate thickness, c/d. In Fig. 9, the critical _ 
flange stress is plotted as a function of this parameter and i. The quantity 
A is a dimensionless parameter defined subsequently. In the inelastic range 
the curve is obtained by assuming (22) that the onset of strain-hardening of 
flange plate is at A = 0. 45, that compressive residual stress of foy/2 (A= v2). a 
exists, and that the transition curve is tangent: to the curves at these two 


‘for 0.45 <x < V2, 


a 


for » > v2, with 


— 

7 4 < 

° 
_ of failure is simplified in that it does not recognize the dependence of the [i 
— 7 different failure modes on each other. It was correct to single out vertical — : _ 

buckling, because it occurs in a direction of symmetry. However, when tor- 

— 

| 

ig 

points, hen the analytical expressions Of the DUCKIINE Curve are ‘ 

om’ 


and | a plate buckling coefficient (3) = 0. 425, Assuming a yield strain ¢ € 
= 33/30,000 = 0.0011, the critical etress can be expressed in terms of the 
plate slenderness c/d as shown on the second abscissa scale in Fig. 9. 
' Above the computed buckling limit the strength of any plate is, in general, 
made up of two parts. First, a favorable rearrangement of forces where ele- 
ments adjacent to the plate participate, and second, the plate’s own post- 
buckling strength. The first part is not significant for the flange plate of a _ 
girder with a thin web. For the second, Winter concludes (23) that only for aa 
an outstanding flange of more than about thirty times the flange thickness can 
a sizable effect of post-buckling strength be expected. In such a range the 
ultimate stresses are less than one half of the yield stress, thus it is of little 
interest in the design of plate girders. Therefore, it is ‘recommended that the 


_ post-buckling strength be ignored completely (as is the established practice _ 
steel structures (15). 


the design of cold- ‘formed light-gage 
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“FIG. 9.- _TORSIONAL BUCKLING OF THE PLATE 


_ AS described previously (1), a test girder was built with an ileal wide 


flange (e/a = 24), Its” ultimate load exceeded the computed critical flange 
_ buckling stress by about 10%. The flange distortions of such a stricture ap-— 
- proach magnitudes under which serviceability is impaired prior to reaching 


_ In designing plate ¢ girders, the compression flange should be ‘made as wide — 


: as possible to increase its lateral rigidity and its lateral buckling strength. 
- But if this is done in excess, torsional buckling of the flange — will al 


place lateral buckling at a lower ultimate stress. 
In order to eliminate torsional buckling as a primary cause — failure, the ae 
e critical stress of the flange plate given ‘by Eqs. 13, should exceed that of 
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lateral buckling ¢ given by 11, in 


avy 
¢ for range, 26, ond ag 


(that is, a ‘flange width- not exceeding 12, the ratio 


_lateral- length to flange would exclude the of 


= 


0 
- AND TORSIONAL BUCKLING OF 
FLANGE 

nts buckling as a primary cause of failure for girder sections under uni- is. 
_ form bending. Should it be desirable to exceed this limit for the flange width- 

: to-thickness ratio, Eq. 15c could be used to find an equivalent column length 
with which the critical stress is found from the basic column curve, Fig. 8. 7 
This | procedure of eliminating plate buckling by specifying ‘the plate- 
slenderness ratio as a function of the column slenderness ratio has been 

2 recommended (16). Inserting the appropriate values into Eq. 3. 3b of f Ref. 16 
results in exactly the expression given for the elastic range, Eq. 15a. 


The curves in Fig. 10 for the inelastic range are —— on the assump 4 


4 _ for inelastic range, = < 26. This co a 

Fig. 10. As indicated, in the conditi 
— | 
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‘ c/a = 8. 5 for oy ty = 33 ksi). This is plotted as — A in Fig. 10 _ If the ci column 
: curve used would take strain hardening into account, the curves of Fig. 10, 
poet form the boundary | between torsional and lateral buckling, would neces- 


Ultimate Bending —Although | the prediction ofa agirder 
segment subjected to pure bending « appears to be > simple because only three 
ae types of compression-flange buckling have to be considered, there 
still difficulties in specifying admissible compressive- 

stresses. This is due to the presence of four r independent parameters which 
influence the result, namely 1/r, c/d, Ay,/Ag, and b/t. Whereas the first two 

__ essentially control lateral and torsional buckling of the compression flange, 

two latter parameters also influence the ultimate moment to 


certain degree. 


72, 


A slender “web burdens the flanges with the stresses which the web cannot 


as a 


_ sionally in terms of the "aes tee My. At very high web ‘uedawenes : 
~ ratios, vertical buckling of the singed would take place manne the extreme 


suming. that a ‘width of the compressed portion of the web equal | to only 30t 


— 
at 

— 

Assuming first that the flange is prev Fig. 11. 
oment M,, is indicated in Fig. 11 
_torsionally, the predicted ultimate bending m u * 


stress pattern would eventually approach point B, corresponding to the - 
aaa According to Fig. 8 (22), strain hardening could take place - 
_ When 8B < 53, because the web could endure an edge strain of more than 12 =: 
without the occurrence of inelastic instability. _A moment equal to the full 
value can be developed when €f ¢<12 €y and strain-hardening can com- 
- mence, hence My > Mp when €¢ > 12 ey. Thus, the ultimate moment will ex- 
ceed or at least equal the plastic moment at 8 < 53. 54. 
Rather than to speculate on the course of this transition from the estab- © ; 


— lished intersections A and B in Fig. . 11, a straight line may be considered to 


—_ with sufficient “eye ' the strength predictions in the range of 


plastic moment is reached at B= = 53 and an effective v width of 30 t is 5 available 


at B = 360. Of greater importance than the exact ma of the transition curve 


.12.— 
INFLUENCED BY THE _PARAME- 
is: the fact that the ultimate bending moment My depends: 5 also on a “shape : 
factor,” p= A w/ Ag, the > ratio of web to flange area, as is born out in Fig. 12. 
_ The range ‘of heeihe in which most plate girders are built is cross- hatched. _ 
"Details of computation for this figure are shown in Appendix I. epoeennae 
«af the information in Fig. 12 is | to be expressed analytically it can be. done 
‘conveniently in the form — 


“in which the coefficient C is ‘the slope of a line in Fig. “a and ie’ the wens 
secting regen of My/M = 1 andthe same line. The constant C could be ex- 
pressed as in 17a, since most are built in the 
Aw/A < 2, the simpler Ex , tm may be a 
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eee that the s stress computation would be made with the — iis. > 
modulus concept, the nominal ultimate bending stress Oy would equal M,/S, 


where S is the section modulus. Since My/S = dy, My/My = Oy/oy. ‘Thus, Eq. 


“The expression as assumes that t instability of the compression 4 


- flange does not influence the carrying capacity. Itis convenient to incorporate _ 
_ the influence of lateral or torsional buckling simply | by r replacing oe with th the 


= oer [1 - 0,005 (a - 


But, unlike the case in the plastic range, there cannot be : an increase of the 
nominal ultimate bending stress when elastic instability occurs because a 
_ stress distribution similar to that for the plastic moment will not take place. ; 
In order to keep the concept simple it is proposed to use Eq. , 18b only forthe _ 

_ range of web-slenderness ratios where B > Bo, namely in the post- buckling : 

_ range where a linear stress distribution according to Navier-Bernoulli would © 

- Defining the post- buckling ra range as one where. Rennideitoen computed bend 

“ing stress o is greater than the critical web stress s (3), the limiting v value al 


value of b/t (Bo) can be as follows: 


wie 

=. 6 E 


(1 - 2 


which = 170 

In this way, the influence of t the two Aw /Ag anc and b/t would be 
restricted in the ordinary plate-girder design to a range of high web- 

a where a reduction of the allowable compressive —_ | 


‘stress is indicated. The reduction, in percentage, of the nominally computed > 
flange stress o can generally be expressed as given subsequently in Eq. 21a, 
the limits of which are obtained from Eqs. 4 and 20. For A7 steel Eq. 21b 
possible reduction, while for a low-alloy high- strength 
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‘Stress in ksi Reduction, in 


-1 170) 


Correlation with Test .—The e which have presented 
for computing ultimate load can be substantiated by the results of nine tests © 
carried out on five different plate girders. The girders and the tests have 
= described (1). In Table 1 some ofthe girder properties are summarized. r ; 


TABLE 1 1. .—SUMMARY OF OF BENDING TESTS WEL WELDED ‘PLATE GIRDERS 


_w 


in inches, 02 


yield 
e 2) 2 


FAILURE | MODES Q 


Below | the cross section the section modulus is given as the ratio of the mo- - 
ment of inertia, I, to the distance from the neutral axis to the extreme com-_ 
* _ pression flange fiber, @q. The effective lateral- -buckling lengths, lk, are chosen 
as 100 in. in tests T1 and 50 in. in tests T2, although the actual bracing dis- 
- tances were 75 and 375 in., respectively. The determination of these values 
was explained in Sec. 2.5 of Ref. 1, where the failure modes are described. 7 
The radius of gyration is computed as r = VI/A where A is the sum of the = 
4 compression flange area and one sixth of the web area and I is the uae ~ 
inertia of this area about the vertical axis of symmetry. 


a The yield stresses of the compression flanges, listed to the right of = 
respective cross sections in Table 1, are somewhat above the value of 33 ksi 
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structural steel with a vi 50 ksi Eq. 21c would apply: 
7 i 
» 
— 
— 
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“ess, for the slenderness ratios involved the ¢ error r in using these figures, 8 
_Girder Gl. —With a slenderness ratio of the compression flange as low — 
19 and a flange -width- -to-thickness ratio as high as 48, torsional buckling — 
"must be the governing failure cause. Referring to Fig. 9 when c/d = 24, Scr 
= 0.56 Oy = 0.56 x 35.4 = 19.8 ksi. This stress need not be further reduced in _ 
accordarice with Eq. 21a since this expression applies only to slenderness 
ratios B = b/t in excess of By = 5.7 VE/o = 5.7 x V30,000/19.8 = 222. Hence 
ou = Ocr and My = dyl/eg = 19.8 x 14,380/25.9 = 11,000 kip-in. For a moment 
arm of 150 in., P,th = - M,,/150 in. = 73 kips, As shown in Table 1, PyX = 81 
giving a ratio ‘of experimental to ultimate load of 1. 11. 
_ Girder G2, Test T1.—Substituting the value of c, d, and], it is found 


‘that Eq. 15¢ holds. Thus torsional failure should not be the buckling mode of © 7 


the flange. Fig. 8 either at l/r 32, or 1/2 c = 8 with = 1. 
«Seine = 0. 97 oy = 0.97 x 38.6 = 37.5 ksi. With Bo = 5.7 VE/o 
= 5.7 v30,0 NOUS. = 161, “the reduction in accordance with Eq. 21a becomes © 
0 05 x 1. 4(185- 161) = = 2%, Hence, My = 0.98 x 37.5 x 14,940/25. 9 = 21,200 kip- 

‘ d ‘in, and P,th = 21,200/150 in. = 141 kips. This would be compared with the 
_ Girder G2, Test T2. —Eq. 15c shows that the mode of failure can be either 
lateral: or torsional buckling. The curves for lateral as well as torsional 
ae. (Figs. 8 and 9) furnish der = dy. . With the 2% reduction according to 
Eq. 21a, it results in My = 0.98 x 38.6 x 14,940/25. 9 = 21 ,800 kip-in., P,,th 


_ Girder G3, Test T1. —With an 1/r = 39, Fig. 8 furnishes Sor = 0.96 Oy 


od 96 x 35. dee 34.1 ksi. Reduced 1% according | to Eq. 21a, « Cu = - 33.7 ksi. AS 
7 ‘pointed out | under the heading “Ultimate Bending Moment: Lateral Buckling” 
é' _ the- stress should be computed at the centroid of the compressive element > 
which provides the lateral ‘rigidity. This was not done in the cases 


of the ‘flange controid, But for girders G3 and G5 eg is taken as the distance - 
from the neutral axis to the center of the tubular compression flange: My 
= 33.7 x 16,200/28. 9 = = 18,900 kip-in. Pat th = 18,900/150 in. = 126 kips. Of 
course, a tubular flange has a very significant torsional stiffness. But with a > 
lateral buckling stress equal to 96% of the yield stress, considering the tor- 
sional stiffness would at most increase the ultimate moment by 4% _ 
,. Girder G3, ‘Test: T2. nine * = 0.99 oy = 0.99 x 35.5 = 35.1 ksi. With 1% re- 
- duction, Ou = 34.8 ksi, My a 8 x 16,220/28.9 = 19,500 kip-in., Py th = 
Girder G4, Test: failed by lateral buckling. According to 
Fig. 8, ocr = 0.98 oy = 0.98 x 37.6 = 36.8 ksi, B. = 163. Reduction: 0. 05 
x 0. 70(388< 163) = 8%. Hence, Cy = 0.92 x 36.8 8 = 33.8 ksi, a = = 33.1 8 x 8 x 13, 420 420+ 
25.7 = 17,700 kip-in., P,th = =11, 700/150 in. = 118 kips. 
Girder G4, Test T2,—My =0. 0.92 ay I/eg= 0.92 x 37.6 17.6 x 13,420/28, 7 = 18, 100 
‘kip-in., Pyth = 18,100/150in.= 121 kips.§ 
a Girder G5. — Computation similar to that for girder G3 gives: 
7 Test T1: My = 0.97 oy a- 0.09) )1/eg = 34.4x 0.91 x 14,710/28.9 = 15 900 aie 
in., Pyth = 15,900/150 = 106 kips, 
"Test T2: My = 0.99 ay (1-0.09)I/e, = 35.1x 0.91 x 14,710/28.9 = 16,300 


16, 


in., Pyth = '16,300/150 in. = 100 kipe. 


| 
as 
f 
- 
- 
4 


a _ The last test result, G5- T2, , exceeds the | predicted considerably because 
the actual I/r ratio was so small that the flange could strain-harden. Thisis _ 
‘also borne out with the load versus centerline deflection diagram shown _ } 
_ Fig. 2.9 in Ref. 1. If the analysis were based on a column curve which takes © 
strain- -hardening into account (Ref. 20, Fig. 27) this discrepancy would be re-- 
All of thes mene design considerations have been concerned with the ultimate 
compressive flange stresses, From these stresses the ultimate bending mo- | 
mo was computed. Except for the possibility of a brittle fracture, vertical, — 
4 torsional, or lateral buckling of compression flange will always be the cause 


_ of failure “under a statically loaded and symmetrically proportioned plate 


girder. Two facts contribute to this. First, it is the web portion adjacent to 
the compression flange which, through its insufficient participation of mem-_ 
stresses, causes an. overstress of the compression ‘flange. ‘Second, 
two equal members are subjected to axial forces of equal magnitude but | 
_ opposite signs, it is the compression member which fails first because of 
instability. Only in the case of girders with smaller tension flange area = 
compression flange area is it possible that the ultimate moment, predicted 
2 one or more of the three different types of compression flange ao 
discussed, , could not be reached. In such case, it is is suggested that the 
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APPENDIX 

Derivation and For the subsequent derivations it is 
sumed that the lever arm between the flange centroids and the distance be- 


- tween the extreme fibers can be put equal to the web depth b. The notation read 
_ defined both in Fig. 3(a) and in Appendix III. The expression wl the — 


of a beam 1 with simply ‘supported end ends is 


Ss mome 
mo 


= girder's ment of inertia kaxis = 
K = torsional constant = = 2 c d?’ + 


The fraction | B in Eq. 24 can be put equal to — This is the case 


thickness d is about 2 ‘the web web and in 


= 

hy. 

& 

4 

— 
A; (2 d2 + —~ t2} 

— 
cand 

— 
— 
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the web area is with the flange area Aw/ Af 
Then the remainder A in the first term is the square of oor (v) (Fig. 5), 
= the numerator equals (0.65 E)2 using Poisson’s ratio v = 0.3, 
‘The radius of gyration r of the equivalent column composed of the com- 


My/My.—With the. notation defined in Fig. 3, the 
=bt ,p= = Ay/Ag, Na = ya/b, np = Yp/b, and the thet the flange 
thickness i is infinitely the computation is as shown below 


b/t = 360,— 


#3 


Section modulus; 


p 


Itimate moment: 


4 

ay a The expression D in Eq. 24 can thus be interpreted as r*, and the product of . 

(29) 

— 

7 

M,, = Np) 12 (53 nb) (34) 

Mu (=f = “shape factor”) ... (36) 


. es shugiitsing assumption of infinitely sr small flange thickness affects the 
result My/My only slightly, since both numerator and denominator are af- 


fected in the Same way and thus cancel possible errors to a large extent. 
y poss g 


EQS. OF REF 17 | 
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ro) 


14 
= L FLANGE Lance FLANGE LARGE FLANGE 
IN COMPRESSION | IN COMPRESSION COMPRESSION IN. COMPRESSION 


A, oa oF FLANGE 


13. —LATERAL BUCKLING STRESS C COMPARED WITH 
PRESSIONS OF REF. 17 
: “same remark applies to cross sections w which differ somewhat ‘from the dou- 
ble symmetric I shape. _ Furthermore, the foregoing ratios My/ My change 
little in the ranges of high web slenderness ratios f. _ Therefore, a value 6 | 


_ = 360 appears to be appropriate also for other ratios of Aw/Ag, even though 7 _ 
* limit of 360 is derivedfor Ay/Ag=O.50nly, 


Lateral Buckling of Unsymmetrical I- I-Sections.—The radius of 


the flange (2c x and done sixth of the web (b t/6) 


— 
— 2 0444 O82 
: @ee 
— 
— 
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PLATE (GIRDERS 
asacolumn 


9.86 x 30, 3,000 2 
Using the k= + that is is, 


Case d = 1 in., 


Ww 

critical stresses apply for 0<K< 1/2, “while for 1/2<K«< 1, Per is 
is equal to the ‘value for x = 1/2 as plotted in Fig. 13, 
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c = = half of flange 


VS 


“a - thickness of flange; 3 
E = modulus of elasticity (30,000 ) ksi); r : 


= distance from neutral axis to extreme in 


= moment of inertia _ corresponding to r; 
= = buckling coefficient; 


1 buckling length ofa 
M= bending moment; é 
- Pe = load on test girders (defined in Fig. 1.1 of Ref. 1); 


and one sixth of the as column cross section); 


radius ote gyration (in 1 lateral ‘direction, by considering Compression flange 


— 


9 


= deflection perpendicular to the ; plane of the web; a 
§ = = ratio of of applied strain to critical ‘strain; 


= nondimensionalized plate deflection; my 

ratio of web depth t to we web thi thickness = b/t; 


= normalized slenderness 1 ratio; a 
‘Poisson’s ratio 0. 3); 
’ pe ratio between web and flange area = Ay 


cripts or Superscripts: 


— flange; 


th = = theoretical; 
u us ultimate; 


= web; and 


= 


= 


othe 
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= 
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_ VARIOUS INSTABILITY MODES OF THE FIXED BASE 


DONALD DA, SAWYER,® A. M., “ASCE, — Larson has out ut the writer’ 
inaccuracy in terminology by correctly distinguishing between a “rocker” and a 
a *roller.” His derivation of the fundamental buckling load for the true roller _ 
case has extended the scope of the paper. It may be observed that an eccentric — 
“frictionless” roller v would not be permissible because any tilt of the top of the m4 
column would cause a tangential force on the roller, Fortunately, friction al- 
_ ways exists to some degree and eccentric rollers could be tolerated in prac- i 
tice. It would be interesting to study the significance of frictionfor this case 

7 with p particular attention paid to the minimum requirements for friction under id 
magnitudes of eccentricity, 

_ With regard to the applicability of Eq. 34 to lift slab conditions, the writer — 
‘de totake issue with Larson, _ Larson states that Eq. 34 describes a struc- 7 
ture in which there is no lateral movement of the slab during buckling and 4 
which there is clearance between the slab collar and the column, The writer — 

- maintains that if a condition of zero sway of the slab will result ina smaller 7 

buckling load, then the system will accommodate itself to this lower buckling 7 

- mode, _ With respect to the clearance requirement, practical reasons dictate 

that a small clearance be provided, as Larson states later in his discussion, 

_ The initial buckling would occur in the small gap. iF After this initial buckling © 

had closed the gap, however, a different condition would exist for which Lar- 

_ son’s point would apply. Actually, this is an academic point, because the most 7 
a critical stage of the lift occurs either at the start of the lift or at the conclu- 

; sion depending on how adequately the slab is restrained from sidesway. Eq. 
“4 (34 would apply to an intermediate stage of greater strength, 


ling equations can be made, This approach is aptly demonstrated by his 


examples that verify the answers obtained by the writer. This is 


of structures loaded by links. a. . Hartmann, F. ASCE, was kind enough to 
urnish an additional reference? as did R. L. Ketter, M. ASCE, 8 R. Dabrowski 


July, 1960, by Donald A. Sawyer (Proc. Paper 2544). 
mi 6 Asst. Prof., Univ. of Florida, Gainesville, 


4 Closure by Donald A, Sawyer 
ae 
| 
i 
= 
| 
| 
| 
— 
fh 7 w Loads,” by H. L. Cox and | 
wh. Stahibau, Bin Mandbucn tur oludium und Fraxis, vol. 1, Verlags-GMBH, Koln, 1956. 


out that his lectures at the ‘neue of Florida gave r reference to a series of 
- letters and other data published in Engineering News Record.! 10- 14 The general 
~ nature of these additional references reinforces the writer’s original statement 
that the link-loaded column obscurely in the the literature. 


m 


do 


9 “Knicksicherheit des Portalrahmens bei Richtungsanderung it im La Last- 
_ angriff,” by R. Dabrowski, Der Bauingenieur, Vol. 35, No. 5, 1960, p.178. 
10 “Bridge Collapse Traced to ” Anonymous, Engineering Record, 
_ oil “Learning from Accidents,” Editorial, Engineering News- Record, July 3, 1952, ma 
12 “How Careless!,” by A. K. Wetzig, E gineering News- -Record, July 17,1952, p. 10, | 
13 “Controversial Rocker,” by T. P. Young, Engineering News-Record, July 24, 1952, 
“Was Viaduct Support Stable ?,” ” by R. Boynton, Engineering News-Record he 
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DYNAMIC RESPONSE OF ELASTO- PLASTIC 


Closure by Joseph Penzien 


_ JOSEPH PENZIEN,® M, ASCE.—As pointed out by Ereminthe elasto- 

characteristics of secondary components suchas walls, partitions, stair wells, 

_and so forth, as wellas their stiffness and strength characteristics provide the 

- apparent ability of many existing structuresto withstand strong motion — 
quakes. This important potas was, of course, a general conclusion reached 
| Eremin has expressed interest in a “repeated stress-strain loop that ~<a 
sulted from a dynamic vibrating force” and has raised the question regarding 


recovery time of steel when strained plastically. Iti is interpreted bythe writer — 

_ that he is concerned withtime-effects onthe stress- -strain properties of steel, . 
These effects as reported inthe literature can be important if strain rates are 
_ sufficiently high. The time between repeated stress-strain loops can be esti- | 


= mated approximately fromthe elastic response of a low damped single degree _ 


system to a random type of input such as used in this investigation, _ Such re- 4 
sponse is an oscillation near its own natural frequency that has an amplitude _ 
envelope and phasing that change in a random fashion. When plastic deforma- » 

“tion is allowed inthe system, thetime interval between repeated stress- strain 

loops will become somewhat larger depending onthe degree of plastic deforma- _ 

_ tion allowed. If one is studying the dynamic response of a specific struc- _ 

ture or element that has sufficiently high f frequency response characteristics, ’ 


then time effects on stress-strain properties of materials cannot always be 
ignored, In the general investigation reported by the writer no attempt was 
_made to allow for such effects because the primary purpose of the investiga-_ 
ton was to show the gross effect that inelastic ‘deformations have on | dynamic: 
response. 
-Eremin’s comment regarding the effect of damping on dynamic | stresses in 
‘his comparison of the results in Figs. 5 to 12 with the results in Fig. 13 ap- 
“pears perhaps to be a misinterpretation of these results. First, it must be 
understood that the forcing function used in obtaining the results of Figs. 5 to 
7 * was the El ‘Centro ground acceleration w whereas the forcing function | used in 


paren exactly with the natural frequency of the system. _ Secondly, the ab- > 
scissa in Figs. 5 to 12 represents whereas in Fig. 13 it t represents 
a= strength ratio 6. In each case (in which = 0), it is observed that the _ 
maximum response must reach a sufficiently high level over the elastic a ‘ 
sponse So > that the energy fed into the system can be absorbed, ‘This observa- 7 


7 
July, 1960, by Joseph Penzien (Proc. Paper 2545), 
6 Assoc. Prof. of Civ. Engrg. Univ. of Berkeley, Calif. 
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mame has discussed the importance of separating the amount of attenation 
received from viscous damping andinelastic deformation, This separation can 
be done in all cases even for those cases involving damping at 0.05, 0.10, and 


: ~—60.15 of critical, because Curves No, 2 in Figs. 8, 10, and 12 also | represent 
: [ - those cases of an infinite yield point with damping of 0.05, 0,10, and 0.15 of a 
respectively. Therefore, iti is to be observed that these curves show 
the attenuation due to viscous damping | only and that Curves 2 through 4 in Figs. , 
5 _ 5 through 7 showthe attenuation due to inelastic deformations only. The atten- 
uation shown by most of the remaining curves includes that due to both viscous _ 


damping a1 and inelastic deformation. 
— The writer is in complete agreement with Blume’s excellent appraisal of the the or 
needs for future research, both analytical and jenpertmental, in the general field 


“gwers are still needed before the stresteral designer can properly consider 7 
of inelastic deformations in earthquake resistant struc- 


=» 


‘Errata.- —On page 85, in the second part of Eq. 


ae be changed to Ugo/ oF . On page 94, in line 1, Housner should be changed to 7 


FREEZING Al AND ‘THAWING | EFFECTS ON ON PRESTRESSED SSED CONCRETE® 


Gg 
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Closure by M. J. and F, E, Musleh 
: 
GUTZWILLER, 13. M. ASCE, and F, E, MUSLEH,14 A, M. ASCE.—It 
should be pointed out that all of the mixes were designed to have an air content 
of 4%+ 2%. The main purpose of the tests was not to compare the durability 
of the different strengths of concrete but rather to examine the effect of the pre- 
stressing on the durability of a particular strength of concrete, namely 5,000 
psi. - The results of the 5,000 psi plain concrete versus the 5,000 psi prestressed 
concrete, as tabulated in Table 3 and Table 4, clearly show that the ‘Garability = : 


The results shown in Table 3 do show, as pointed out by Kunze, that the 3 ,000 

aw plain concrete (if specimen D-2 was omitted) has a higher durability factor 
thanthe 5,000 psi plain concrete, or 5,000 psi prestressed concrete. The writ-— 

a ers feel that the data obtained in this study are too limited to drawa general _ 
conclusion on comparison of concrete durability of different mixes. Within the 
_ range of variables inthis experiment, it canbe said that the richer mix showed 
a higher rate of of deterioration the leaner mix, 


- 4 October, 1960, by M. J. Gutzwiller and F. E. Musleh (Proc. Paper 2¢ 2633). 7 
Prof. of Structural Engrg., Purdue Univ., Lafayette, Ind. 


14 Design er. , Pierce, Gruber and Beam, Inc., Indianapolis, ind. 
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comments. A review of electrical engineering literature shows 
practical applications ¢ of the ‘proposed in that field; it is hoped that 
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— & October, 1960, by Kurt H. Gerstle (Proc. Paper 2634). a 
8 Assoc. Prof. of Civ. Engrg., Univ. of Colorado, Boulder, om 
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Discussion by George B. Begg, Jr. 


32 JR., _ ASCE.— _This paper has 
insights into the problems in structural safety confronting the 
_ present-day designer. Brown strongly develops the need of integrating broad = 
economic viewpoints withthetechnical approach to structural-design assump- 
tions. To admit and welcome a dialogue with participants fromdisciplines other = 
than the design office, university, and testing laboratory is a long overdue in- pe 
« ~ novation . However, this writer doubts that the question is qereenily as “open” 
- the author infers in his paragraphs on “Structural Life” = = a 
; - Brown suggests that the design ofa static structure may be approached on 
‘either the “definite life span” or “monumental” basis. Limiting comment to 
‘ the design of buildings and hydraulic structures, the writer has seldom observed 5 
the definite life span as operative ROF does there seem to be any gathering 


vertical extensions in buildings in which no such provisions were contemplated J 

in the original design, ‘The unprovided for vertical extension is deliberately 

" chosen because of the great ‘demands it makes on ‘concepts of structural safety. fe 
Ina more subtle, but just as real, way it offers considerable challenge tothe _ 


studies on apr tee vertical extensions to existing buildings. Such studies re- 
quire the combined appraisal of the architect, mechanical oS 

and structural engineer, . From a structural viewpoint the request cooetany 
inquires, | “Can we put one, two, or three stories on this structure for such and | 
such live loads?” ‘Usually the structural engineer is told by those concerned a 


ne 


— A 


tionally the only desiredarea and if we cannot have it there we probably cannot "es 


¢ Confining structural comments to the context of this paper, several features t 


reoccur in these feasibility studies: ta | iy 


hee of Structure. —The - buildings under cc consideration ‘most st often f: fall in th in the 
20 yr to 50 yr old age | bracket. This means that they are experiencing the au- 
thor’ s “fatigue, wear, and corrosion [that] necessitates a definite life toa 


structure. Thus, whereas the vertical extension is the struc- 


December, 1960, by C. B. Brown (Proc. Paper 2678). 
_ 32 Chf., Design and Investigation Sect., Structural Engrg. Branch, Pub. Buildings 


To document this position the writer would like to discuss the —_— of ¥ 
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- _ The Design- -Functional Records of the Structure. —In this area let us con-_ 
| sider a specific design office. The Public Buildings Service (and its prede- 
cessor agencies) is over 100 yr old. As a matter of policy, and vindicated by 
% experience, the structural engineering records have always been maintained 
at a high degree of completeness - perhaps unequalled in any other structural 
ie office. Therefore, in a sng study we can utilize the original de- 
sign computations, 
and building management history in reaching our conclusions. However, in - 


these extensive the ability to to understand the he existing 


* a. Conflicting Data.—The results of a design computation may be ‘contra- 
dicted by the requirements of the contract drawing, =” 
Operational History. —The design assumptions and execution often ap- 
pear satisfactory, but a building superintendent can be continuously reporting 
but unanalyzed cracks and deflections, = = 
Absence of Original Design Team. — Computations and drawings that were 
es enna clear to a designer preparing them 30 yr ago, can confuse a later in- 
- -vestigator, particularly where a vital assumption was not explicitly recorded, 
The absence of the original design team is pertinent to this discussion because 
7 it precludes any documented position onthe proposed life span of the structure. 


» W 


_ Investigator - Client Relationship.—Two important ground rules governing 


feasibility studies ; should be recognized herein, The first is that such studies 
are invariably limited in time - with the resulting value judgments on what to 
gg investigate and how far. The second rule to remember is s that the structural 4 


“@ be. Generally, the client asks what he thinks is a simple question only to re- 
ceive a very complex answer or a series of countering questions. And the an- 
_swer or the questions rarely have a unified note. When confronted by apositive 
; answer from his architect, mechanical engineer, and estimator the client is 


_ from his structural engineer. . In those cases in which the e structural e engineer 7 
is swimming against the tide, it is usually necessary to. pursue a factual but 
_ simplified argument aimed at the common denominator, for example “the ex- 

isting foundation system will not accept even a minor part of the | new loads that 


will result from the proposed 
a The preceding review of vertical extensions in the advanced lif. life of of | a  build- 


‘ing leads to some questions, A similar group of questions would have emerged 

if remodeling instead of extension had been dev eloped. The writer is convinced — 

that the answers to the following questions are either r implicitly statedin Brown’s _ 


paper or or naturally grew out of it: 


Safety.— — How for safety i ina na structure that is s substantially 
: altered late in its operational life? Or again, what modification of variables as 
_ developed in the discussion of the ratio and probability methods is required 
- for a vertical extension that uses new materials and stresses and is executed - 
7 by a new design team but rests on an ‘old structure that was not designed for 
it? what probability of such occurrance should be assigned in a new 
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The Definite Life Span Assumption.- a definite life span 
to compete withthe monumental design in building inthefuture? Will the client 
agree (except in an unusual speculative case) to a building that is designed to 
= be replacedin a specific number of years? Why should we look for a leveling- 4 
off of the current practice of vertical and horizontal extensions plus remodel- 
ing? If the public ingeneralisto participate in judging the design assumption, ~ 
what type of preliminary agreement onterminology | and premises must the de- 
signteam agree on ? Should the architect or the engineer be the moderator in 
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‘Discussion by M. Gregory 


GREGORY. 9_ The author has presented an paper on the lat- 
- eral inelastic buckling of a tied arch, It is quite permissible to make a South- 
well Plot onthe observed lateral deflections inthe inelastic range. An upward» 
curving line is expected, as shown in Fig. 6, but how one should interpret the : 


inverse slope of this curve (and it varies considerably from one end of the | 

curve to the other) in terms of the critical load, is difficult to say. If itis 
realized that the tangent-modulus critical load is a property of an inelastic 
no initial imperfections either of form or (of loading, one 


re 
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BE in practical | cases, though it is interesting to note some ‘small measure 


when en imperfections are minimized, as the author has caretalty 


oxtiien! load furnishes a solutionto allinelastic instability problems. It is still 
_ a critical load, a property of the perfect structure. The problem of the prac-— 
tical structure is much more difficult and intractable, It is basically a prob- 
lem of strain configuration (as is indeed emphasized even inthe philosophy one 
uses indiscussing the relevance of the tangent modulus). To obtain vital fun- 
_ damental information on this problem it would be useful to measure strains at — 
_ critical points as the arch is loaded, in order to know which portion of the 
stress- curve applied at stage of the loading. > 
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2 a a havior of the structure under load only if its imperfections were extremely [im 
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‘SAFETY, RELIABILITY AND STRUCTURAL DESIGN@ 

_ Discussion by Jack R. Benjamin and Theodore c. Zsutty 


JACK BENJAMIN, 10 M, ASCE, and THEODORE c. ZSUTTY! — This 
“most interesting paper contains many general concepts of great value in struc- 


turalengineering. The paper is concerned with important aspects of the math- 


ematical theory of probability. It must be understood that decisions cannot be 
rationally made purely on the theory of probability. - Two equally important 
additional bodies of theory must be considered, _ First, the statistical nature 
of the problem must be solved, This involves the values of parameters, fit of 
‘Getz hypothoste tests, and forecasting concepts, among other items 


Probability supplies the initial mathematical framework, statistics interprets 


nature, and decision theory adds value concepts. 


_ In probability s studies it must be e emphasized that if R and Ss are random vari- 


2, ig only the first step. The variance of v and associated PF is ‘more 
critical inthe decision analysis thanthe characteristic long run value, Finally, 


real structures: involve many - modes of behavior each with its own v and Pr. 


These structures can be designed to a constant v or a constant Pf. The Pal 
designs differ significantly. It can be readily shown that constant Pp may in- 
volve widely varying load factors throughout the structure, Thus current load 
factors, based on judgment, cannot possibly yield a structure with constant 
probability of failure unless it is a simple beam under one loading condition, 
Even two _independent loads applied | to a simple beam demands two different 
load factors to maintain a constant probability « of failure, _ The differences in. 
+ variances of the loads and resistance are responsible for the phenomena, - 
; The ‘first: basic question of how to determine a realistic probability distri 7 
bution of load must be resolved, The log- ‘normal distribution for loads is od 
venient for computation purposes, but it has certain characteristics that re- 
‘strict its use in real situations. First, the distribution is limited at zero ol 
‘rather than at a constant dead load of major magnitude, . The “Weibull distri- 
bution” might be used to overcome this difficulty. Second, the tail of the i 
normal is very extensive; a consitionthat gives forecasts of very large iam 
Elements of control always exist with man-made loads and, thus, forecasting 
from sample data for these loads into regions of low probability is hazardous 
indeed. A similar difficulty applies also to extreme value distributions when 
small data samples, are used to estimate the parameters of the distribution. 


— 


_ The small data sample has basic limitations of description such that the accu- 
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rate estimation of parameters required for forecasting is not yontthie, Also, 
lus the sam- 

"pling variance of estimated value. The sampling variance increases as the 

_ magnitude of the forecasted load increases beyond the range of values in the 
data simple. Until these variances are considered and confidence limits es- : 
_ tablished, the equations, charts, and tables given in terms of theoretical dis- 
Perseromte will remain as exercises in mathematics without practical use. It 

should be noted that v and PF are forecasts, 

_ The second basic question of howto determine the distribution of structural 

_ resistance has a better chance of being answered because many test records 
of resistance of structural elements are available. In order to use probability 
of failure as abasis for design it will be necessary to have a formula that will © 
predict some central value such as the mean of the resistance of a structural 
_ element and that has a known probability distribution of the population of re-— 
sistance values ¢ about this central value, The confidence intervals of the pre-— 
“dicted central value and the » distribution parameters should also be known. it 
4 is doubtful whether any of the present formulas for resistance accurately pre- _ 
dict a central value, and if it ‘is even more doubtful if anything is accurately _ 
~ known about the distribution of actual values about the ‘predicted value. If rec- 
= of alarge number of widely varied tests of resistanceof a structuralele- — 


are available, one way of obtaining the required formula is by the 
cal methods of multiple regression analysis and dimensional analysis. 
event that empirical methods proveto be the only way of obtaining the formulas 
‘perimental design, will be required for nearly all structural elements. In con- 
junction withthis program, studies are required to determine the variation and | 
covariation of member dimensions and strengths due to workmanship errors. 
The implication that the probability distribution of resistance of a structural — 
element is the same as_ that of the material of which it is composed may be 7 
highly erroneous in practical cases. For example, the ultimate moment of an_ 
under-reinforced concrete beam is quite insensitive to large variations of its — 
concrete cylinder strength, Even if the concrete strength has an extremal dis- 
_ tribution due to poor quality control, there is no reasonto believe that the ulti- 
Be moment resistance has the same extremal distribution, Similarly there” 7 


: resistance, then a large program of testing, guided by the principles of ex-— 


fabrication stresses and points of stress concentration, would have adistribu-_ 7 
tion of plastic moment resistance the same asthatof a yieldstrength specimen > 
of the same steel, In general, each structural element may have a distribution 


of resistance that differs from that of tests of its base material or materials, 
The author is to be congratulated on his continued pioneering effort. The 
= _ writers believe insufficient data are presently available to allow ny eanpenl 
tural design using statistical concepts. _ However, it is time that the long-range 
educational process toward this goal be initiated; this paper balk an onsennng 
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